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SYNOPSIS 


Olefins and diolefins are important intermediates in 
the petrochemical industry, and the future promises further 
substantial increases in demand. Bismuth molybdate and tin 
antimonate catalysts are very versatile and are effective in 
several different processes such as oxidative dehydrogenation 
of butenes to butadiene, oxidation of propylene to acrolein, 
amm on oxidation of propylene to acrylonitrile etc. 

Several catalysts have been formulated for the 
abstraction of ally lie hydrogen from butenes or propylene. 

These catalysts, however, are inefficient in the abstraction of 
first hydrogen from a saturated hydrocarbon molecule like 
n-butane. Abstraction of first hydrogen atom from n-butane 
molecule involves the rupture of a strong a bond between the 
hydrogen and the carbon atoms. The hydrocarbon molecule thus 
has a tendency to burn off completely rather than dehydrogenate 
to yield butenes and butadiene. 
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The aim of the present work has been to develop a 
suitable catalyst for the oxidative dehydrogenation of n-butane 
to butenes and 1, 3-butadiene. Preliminary experiments were 
carried out on a molybdic oxide -aluminium phosphate catalyst to 
investigate the feasibility of oxidative dehydrogenation of 
butane. Experiments were performed at different conditions of 
temperature, contact time, and catalyst composition. The 
catalyst showed some activily for the oxidative dehydrogenation 
of n-butane. 

In order to improve the activity of the molybdic oxide- 
aluminium phosphate catalyst, bismuth was introduced in the 
catalyst, and the catalyst activity was investigated further. 
Three types of aluminium phosphate-bismuth molybdate catalysts 
were chosen, namely, 

1. Bi20^.3Mo0^ (a -phase) + AlPO^ 

2. Bi20^.2Mo0^ (P -phase or Erman phase) + AlPO^ 

3. Bi 20 ^.MoO^ (Y -phase or koechlinite) + AlPO^ 

These catalysts were characterised in terms of chemical 
compositions, surface areas, characteristic X-ray diffraction 
reflections, infrared spectra and thermal stabilities (BTA, TGA 
and DTGA). X-Eay diffraction studies and infrared spectra of 
Bi 20 ^. 3 Mo 0 ^-AlP 0 ^ catalysts indicated its unsuitability as a 
catalyst, since it had monoclmic crystal structure with too many 
edge shared Mo - 0 octahedra. The other two phases viz. 
Bi20j.Mo0^ and Bi20^.2MoO^ containing corner shared Mo - 0 
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octahedra were investiga'ted "to establish the best conditions for 
oxidative dehydrogenation. 

Optimal conditions for oxidative dehydrogenation were 
established using techniques developed in applied statistics, 
namely, exploration of response surface using orthogonal and 
rotatable experimental designs. Various factors, and their 
region of determination for finding the optimal were 

Temperature 400-5 00 °C 

Catalyst to feed ratio 0.25-2.5 

W/P , 

(0.25-2 sec contact time) 

Butane to oxygen ratio 0. 5-2.0 mol/mol 

mol bismuth molybdate 0-100 

per 100 mol aluminium 

phosphate 

Oxygen to nitrogen ratio 0.5 mol/mol (held constant) 

Bismuth molybdate p- and y-pbases 

(qualitative variable) 

Optimization parameter was chosen as the yield of butene 
and butadiene. Pactorial experiments were designed with temper- 
ature, W/P, butane-oxygen ratio and per cent bismuth molybdate 
as quantitative factors, and the phase of bismuth molybdate as 
qualitative factor. 

To start with, fractional factorial experiments (one -half 
replication) were designed, and a first order response surface was 
fitted to the data. Method of steepest ascent was used to move 
towards the optimum. In the near-optimum region, a central 
composite design was used, and experimental results were fitted to 


g catalyst 
mg mol/min 



a second order model for the response surface. The exact optimal 
conditions were located using the second order equation, and the 
second order model was subjected to canonical analysis to examine 
the exact nature of the response surface. 

Mechanistic modelling of the dehydrogenation process 
was then attempted* Redox, Rideal-Eley and Langmuir-Hinshelwood 
models were constructed and point selectivity defined as the 
ratio of rate of formation of butene to that of butadiene was used 
to Identify plausible models. Finally recommendations for future 
work have been proposed. 



CHAPTER 1 


IHTRODUCTIOH 


1.1. Introduction . 

Naturaal gas and associated gas, presently used as 
sources of energy can serve as potential feedstocks in the 
production of petrochemicals like ethylene, propylene, butenes, 
butadiene, maleic anhydride etc. Use of such a valuable 
material as a source of energy when other cheaper and abundant 
energy sources are within easy reach could bring forward a day 
when such natural resources are extinct. There has been a 
growing concern in the industrial as well as political circles 
to utilise these expensive petroleum products very judiciously, 
with a view to conserve this precious, limited energy source. 

1 .2 . Utilisation of Gaseous Hydrocarbons — Indian Scene . 

3 

About 400 million m of associated gas comprising of 
75% methane, 8% ethane and 12% propane and butane will be 
available on the shores of Gujarat and Maharashtra by 1980, 
from the off-shore oil production platform of Bombay-High [55] . 
Claims made by the state governments and the industries have 
led the Indian Government to evaluate the alternatives available 
for the best use of Bombay -High gas [36]. The methane rich gas 
from Bombay-High can supply the natural gas req.uirements for 
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three urea plants of 1500 tonnes per day capacity, and already 
there is some thinking regarding the use of Bombay-High gas as 
a feedstock for the fertiliser plants based in Bombay and 
Gujarat. 

Methane may also be used to produce methanol, 
acetylene etc., whereas ethane, propane and butane can be used 
for the production of olefins, butadiene, maleic anhydride and 
a host of other petrochemicals [61, 114]. 

The present study is aimed at the utilisation of 
n-butane. The dehydrogenation route provides a means for the 
synthesis of other hydrocarbons that yield ultimately various 
useful products. The dehydrogenation of n-butane has been 
discussed briefly in the following section. 

1.3. n-Butane Dehydrogenation Processes . 

Conventionally, n-butane is dehydrogenated to butenes 
over a chromia-alumina catalyst [34, 41, 42, 55] at temperatures 
of 500-550°C. The reaction is equilibrium limited [23] , and 
the catalysts tend to loose their activity rapidly due to 
severe coking of the catalyst pellets. 

Halogens and halogen based compounds were used in the 
reaction in order to alleviate the above difficulties in the 
dehydrogenation process [42, 61, 98, 115]. In such processes 
when iodine was the halogen used, butane was dehydrogenated to 
1 , 3-butadiene in a single step with very good selectivity and 
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yield. However these processes did not gain much commercial 
importance due to the prohibitive cost of iodine. 

Butane can also be dehydrogenated to butenes and 
butadiene over various types of solid catalyst--' [49, 61, 115] 
using mixtures of butane and air. Such a process appears very 
attractive since the oxygen present in the reaction media 
prevents the coking of the catalysts and secondly the dehydro- 
genation of butane to butadiene can be realized in a single step. 

Literature on oxidative dehydrogenation of n-butenes 
to 1, 3-butadiene [10-21, 45, 61, 104] indicates that the bismuth 
molybdates catalysts are very active and selective compared to 
other catalysts. It was generally observed that the molybdic 
oxide provided the activity in the dehydrogenation while bismuth 
oxide was responsible for the selectivity. 

The feasibility of the one -step oxidative dehydrogenation 
of n-butane over a molybdic oxide -aluminiiim phosphate catalyst 
was established by Gaspar and Pasternak [49]. The activity of the 
molybdic oxide -aluminium phosphate catalysts could probably be 
enhanced by the addition of bismuth to molybdic oxide. 

1.4. Object of the Present Investigation . 

Having looked into the feasibility of the oxidative 
dehydrogenation reaction over a molybdic oxide -aluminium 
phosphate catalyst and considering the reported influence of 
bismuth oxide on the selectivity of the reaction, the present 
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investigation was directed towards the study of oxidative 
dehydrogenation of bismuth molybdate -aluminium phosphate catalysts. 

The structure of the catalysts were determined first 
and their thermal stabilities were investigated before actually 
studying the butane dehydrogenation reaction over them. The 
results are given in Chapter 4 . 

The effects of process variables like temperature, 
contact time, butane to oxygen ratio and catalyst composition on 
the yield of butenes and butadiene were studied using response 
surface methodology in oi*der to establish the optimum process 
conditions. These results are presented in Chapter 5. 

Various mechanistic models have been proposed in 
Chapter 6. Point selectivity (defined as the ratio of rate of 
formation of butene to that of butadiene) method has been used 
to identify plausible models. Pinally recommendations for 
future work have been proposed in Chapter 7. 
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GHAPOEE 2 

SURVEY OE BELATED WORK AED SEIECTIOH 
OP CATALYST POR THE PRESENT WORK 

Several processes based on different raw materials 
have emerged for the production of olefxns and diolefins. Due 
to the relevance and importance, literature on the production 
of n-hutenes and 1 , 3 -butadiene is reviewed in the following 
sections. 

2.1. Sources of Butenes and Butadiene . 

Butenes and butadiene are produced by the following 

methods : 

1. By cracking of hydrocarbons. 

2. Manufacture from ethylene, propylene and ethyl alcohol. 

3. By thermal dehydrogenation of n-butane and n-butenes. 

4. By oxidative dehydrogenation of n-butane and n-butenes. 

Various aspects of the above processes have been 
reviewed by several authors [11, 42 , 44 , 55 > 61, 104, 114, 

115, 124]. 

2.2. Cracking of Hydrocarbons . 

In the petroleum industry, naphtha is cracked as a 
means of obtaining raw materials for the manufacture of 
organic chemicals. Cracking is directed mostly towards the 
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production of ethylene and propylene but invariably some butenes 
and butadiene form as side-products [11, 42]. Cracking is 
normally done at low pressures with steam as a diluent at a 
temperature of 700-750°C. lower hydrocarbons like butane can 
also be cracked to yield butenes, however the yields are very 
low [94, 116]. 

2.3. Manufacture from Ethylene. Propylene and Ethanol . 

Commercial processes for the synthesis of butadiene 
based on ethylene and propylene have been developed, but it is 
unlikely that these will compete economically with the processes 
based on butenes and butane. Butadiene can also be obtained 
from ethanol. 

Ethylene dimerization over a Ziegler catalyst such as 
aluminium triethyl with nickel acetylacetonate yields butenes. 
Alternatively ethylene may be oxidised to acetaldehyde, which 
may be further converted to acetaldol by aldol condensation. 
Eurther hydrogenation and dehydration yield 1, 3 lutadiene [42]. 

Propylene disproportionation over MoO^ -C oO-AlgO ^ 
catalyst at about 200°G and 30 atm pressure yields butene and 
ethylene as products. Ethyl alcohol dehydration over 
Al20^-Zn0, MoO-GoO catalyst gives 1, 3-butadiene [39]. 


2.4. Thermal Dehydrogenation of n-Butane and n-Butenes . 


Processes which involve reactions of the type. 
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C^H^q Catals^t^ ^ A H = 28.5 kcal 

C^Hg C^Hg + Hg Z1S.H = 28 kcal 

are referred to as 'Thermal Dehydrogenation processes’ . These 
reactions are endothermic and the yields of butenes and 
butadiene are constrained by thermodynamic equilibrium [23]. 

Practically all operating commercial butadiene 
processes are based on the thermal dehydrogenation of n-butane 
or n-butenes or mixtures of these. 

1 . Thermal Dehydrogenation of n-Butane . 

Dehydrogenation of butane to butenes is conveniently 
carried over a chromia-alumina catalyst [26 , 34 , 41, 48 , 55, 

56, 73, 91 ]. The chromia-alumina catalyst deactivates during 
dehydrogenation due to coke deposition, and demands frequent 
regeneration. The summary of the work reported in literature 
IS presented in Table 1 . 

The most important thermal process using n-butane as 
the chief raw material for large-scale manufacture of butadiene 
is the Houdiy process [42], 15-20?5 chromia on activated 

alumina is used as a catalyst at a temperature of 620-650°C 
and a pressure of approximately 0.2 atm without any diluent 
gas. Several reactors are used in swing operation to tackle 
the problem of coking. 
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Table 1. THEmAl DEHTDROGENATIOH OE n-BUTlUE TO n-BUTEHES 


Catalyst 

Temperature 

°C 

Pressure 

Atms. 

Contact 

Time 

Seconds 

Conversion 

% 

1 1 

00 

I > 

CM 

O 

350 

1 .0 

- 

7.5 

[26] 

530 

1 .0 

0. 5-2.0 

- 

Al 20 ^-Cr 20 ^ [ 41 ] 

569 

1 .1-2.7 

0 . 25 - 9.0 

4-36 

Alg 0^-20% CrgO^ [55, 

56] 447-550 

0 . 19 - 3.331 

i - 

0 

• 

to 

11 ^ 0 ^- 22 % Cr^Oj [ 91 ] 

501 

0.06-1 .8* 

10.0 

3-6 

Al20^-Cr20^ 

510-550 

0.1 -0.5 

- 


+ Na20, Li02 [34] 






* 


Nitrogen diluent. 
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11 . Thermal Dehydrogenation of n-Butenes . 

Several processes are in existence differing mainly 
in their choice of catalyst. The feed -stock used is a 
olefin-rich hydrocarbon distillation cut, from which isobutene 
has been removed by acid treatment. Steam is the most commonly 
used diluent, and dilution ratios of steam to hydrocarbon upto 
20:1 are used. Considerable effort has gone into the development 
of selective catalysts. Some such catalysts are shown in 
Table 2 . 

Catalyst deactivation due to coking is a severe 
problem in the dehydrogenation of butenes to butadiene. Those 
catalysts which contain potassium salts are to some extent 
self -regene rating, in particular the shell 205. Potassium 
salts have a catalytic effect on the water-gas reaction 

c + H2O ^ CO + 

and in this way coke is removed from the catalyst bed. The 
Dow catalyst is at a disadvantage since it requires air 
regeneration and has a short on-stream life. It has the 
advantage, however, of a significantly higher selectivity. 

Dehydrogenation of butenes to butadiene over 
chromia-alumina catalysts demands temperatures higher than 
that needed for dehydrogenation of butane over the same 
catalyst [55, 57, 120]. Coking problems and poor selectivity 
to butadiene impede their acceptance as a suitable catalyst 



Table 2. COMPAEISON OP n-BUTEEE DEHYDROGENATION 
(THERIIAD) CATAiYSTS 



Esso 1707 
[22, 103] 

! 

! Shell 205 
! [125] 

t 

} Dow Type B 

i [92] 

\ 

C hemic al 
composition 

18.4 EegO^, 

72.4 MgO, 

62.5 

2.2 Cr^O^, 

0aglIi{PO^)g 

+ Cr20, 


4.6 CuO, 

4.6 KgO 

53.3 KgCO^ 


Temperature °C 

625-650 

650 

625-650 

Steam: butane 
ratio 

14:1 

18:1 

20:1 

Conversion % 

28-35 

37 

28-40 

Jk 

Selectivity % 

75-70 

76 

86 


Selectivity is defined as the ratio of moles of butenes 
forming butadiene to the total moles of butenes reacted. 



OXIDATIVE DEHYDRO&ENATIOIT 


Thermal dehydrogenation of n-hutane or n-hutenes 
yields hydrogen as one of the products. Build -up of hydrogen 
concentration in the product stream imposes a thermodynamic 
limitation on the equilibrium conversion [23] . If the 
hydrogen is oxidised during the course of dehydrogenation, the 
equilibrivun of the dehydrogenation reaction can be shifted in 
a favourable direction. In addition, sufficient heat can be 
generated for the dehydrogenation reaction which otherwise is 
endothermic in nature. 

Various types of hydrogen acceptors, like — chlorine, 
bromine, iodine, sulphur and selenium may be used. The 
dehydrogenation reaction with these hydrogen acceptors is 
exothermic except for the case of iodine, for which it is 
endothermic . Suitability of various hydrogen acceptors has 
been discussed by Vadekar and Pasternak [123] . Extensive 
reviews of the oxidative dehydrogenation processes have been 
given by Voge and Adams [124] in 1967, Sittig [114] in 1968, 
Skarachenko [115] m 1969, Sachtler [104] in 1970 and Hucknall 
[61 ] in 1974. 

Thus, oxidative dehydrogenation can be effected in 
two ways. The first is the direct oxidative method, in 
which the hydrocarbon to be dehydrogenated is mixed with an 
oxygen containing gas, usually air. This mixture, to which 
steam may be added, is then passed over a suitable catalyst. 
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The second method makes use of halogens and sulphur. The 
gaseous mixture containing hydrocarbon, air and halogen is 
passed through a catalyst bed. In some process modifications, 
the halogen may be reacted in a separate stage with the 
hydrocarbon and at a later stage, the oxygen is used to 
regenerate the halogen from the resultant halide. 

These processes are discussed in detail below 
for butane and butenes. 

2.5. Oxidative Dehydrogenation of n-Butane . 

The production of diolefins by the oxidative 
dehydrogenation of alkanes is generally not successful due 
to low reactivity of these compounds. At temperatures below 
300°C, aldehydes, alcohols and oxides of carbon are formed. 

At temperatures of 550°C and above cracking of butane dominates. 
The yield of oxygen containing compounds reduces and the same 
of ethylene and propylene increases as the reaction temperature 
is increased from 300°C to 7 OOOC, whilst, the yield of butenes 
passes through a maximum [115]. 

The yield of oxygen containing compounds like 
aldehydes and alcohols decreases when the oxygen to hydrocarbon 
ratio IS reduced. Water is the sole oxygen containing compound 
at temperature of 500°C and oxygen to hydrocarbon ratio of 
0.1 [ 10 ]. 
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Levin et al. [81] used a number of oxides in their 
study and arranged their activities as follows: 

Ni> Mo> Go>W> Pe > Gu^ Mn>Gd> Sn> Pb^^Gr^^Bi 

Patent literature on oxidative dehydrogenation of 
butane is not very extensive. Varied reactor designs have 
been used for the purpose. McDonald and Ifcintyre [82] have 
used fluidised bed reactor whereas Lester [80] designed a 
reactor in which air and butane streams were separated by a 
porous concentric tube containing an oxygen transfer agent. 

The hydrogen generated during butane dehydrogenation diffused 
across the porous wall, into the air stream, thus reducing 
the hydrogen partial pressure and increasing the yield of 
butenes and butadiene. Table 3 provides a glimpse of the 
results of a few important patents and literature. 

Inorganic complexes used as catalysts allow 
dehydrogenation reaction to be carried out at lower temperatures. 
Kovaleva and Tulupov [71 ] used the complex compound 
[ (GH^ )2K'^G1’'] 21 Catalyst in a melt containing 

53% KWO^, 40% NaN 02 and 1 % BalfO^, with a catalyst concentration 
of 0.1-5%. Oxygen and butane mixture were bubbled through 
the catalyst containing melt at 170-200®G. The product was 
only butenes with undetectable side -products. 

1 . Halogen Based Dehydrogenation . 

Iodine can be used for the dehydrogenation of 
hydrocarbons, especially butane [98]. The endothermic reaction 





Ip 


^4^10 2 I 2 ^ + 4HI 

proceeds at temperatures above 400°C as studied by Raley et al. 
[98] and recently by Kuz'michev and Skarachenko [74]. 

In the presence of a metallic oxide capable of 
reacting with hydrogen iodide [ill], the reaction was found to 
go to completion with negligible formation of undesirable - 
by-products such as acetylenes. The iodine could be easily 
regenerated by oxidising the metal iodide with air or oxygen 
to give metal oxide and free iodine. The two steps can be 
combined into one by passing hydrocarbon with air and less than 
stoichiometric amount of iodine over a suitable metal oxide 
catalyst like PbO, ligO [111], and CaO [96]. In this system 
the elemental iodine is regenerated continuously to carry out 
the reaction to completion. 

Bromine and chlorine can also be employed in place 
of iodine in the oxidative dehydrogenation process. The 
reaction is exothermic in the case of chlorine, and mildly 
exothermic with bromine. These processes use mostly alkali 
metal halides for catalysing the reaction [118, 11 9 ]. 

Use of mixtuore of halogens gives better yield and 
selectivity as compared to the use of a single halogen in the 
oxidative dehydrogenation. Results of some of the experiments 
reported in literature are shown in Table 4. 

Halogen-based processes offer many advantages over 
conventional methods, notably high conversions and select ivit 10 s , 



Table 4-. HALOGEN AND SULPHUR BASED DEHYDROGENATION OP n-BUTANE AND n-BUTENES 
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and continuous operation without the need for catalyst 
regeneration. Iodine offers several advantages over other 
halogens, especially in the case of butane dehydrogenation. 
Iodine unlike chlorine and bromine does not react with butadiene 
at elevated temperatures and there are, thus, no undesirable 
side reactions. 

II . Sulphur Based Dehydrogenation . 

Sulphur vapours were used as hydrogen acceptor by 
Rasmussen et al. [100] in dehydrogenation of butane, at a 
temperature range of 650-760°C. Nascent sulphur is an 
efficient hydrogen acceptor as compared to sulphur vapours. 
Pasternak and Vadekar [93] have used hydrogen sulphide with a 
mixture of butane and oxygen to effect oxidative dehydrogenation, 
where the hydrogen sulphide combined with oxygen to produce 
nascent sulphur and water. Small quantities of halogens mixed 
with hydrogen sulphide significantly improved the conversion 
and selectivity of the dehydrogenation reaction. A molten 
salt bath of LiCl-KCl-MnCl 2 catalysed the reaction equally well 
when hydrogen sulphide, oxygen and butane were used in the 
absence of halogens [51] . Table 4 shows some of the results 
found in the literature. 

III . Other Hydrogen Acceptors . 

Ethylene and propylene have been used as hydrogen 
acceptors in dehydrogenation of butane over C hr omi a -Alumina 
catalyst [66, 77-79]. Absence of oxygen in the reaction media 
leads to coking of the catalyst. 
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2.6. Oxidative Dehydrogenation of n-Butenes . 

Numerous catalysts for oxidative dehydrogenation of 
butenes have been reported in the literature [44, 53, 69, 114, 
115, 124]. Ihe catalysts recommended are generally mixtures 
of oxides of elements from group Y (Bi, Sb, As, P, V ) and VI 
(Or, Mo, W) of the Periodic Table. 

i. Oxide Catalysts . 

Pure metal oxides themselves are not very selective 
and give very poor conversions. The addition of oxides of 
metals from neighbouring groups often greatly enhances conversion 
and selectivity of metal oxide catalysts. 

Those catalysts based on Bi-Mo, Sb-Sn and Pe-Sb are 
probably the most useful and have been studied in detail in 
literature. 

Bismuth Oxide -Molybdenum Oxide Catalysts . The oxidation of 
butenes over bismuth-molybdates has been the subject of 
considerable study. Over this catalyst, oxidation proceeds 
with a high selectivity towards butadiene, even at high 
conversions. Adams et al. [1, 2] and Bleyenberg et al. [25] 
have reported selectivit les of 90-95% at conversions ranging 
from 20-80% in oxidation of 1 -butene at 460°G . Cis-2-butene 
and trans -2 -butene are about 5 times less reactive than 
1 -butene [2, 63] , trans -2-butene being the least reactive of 
the three. 
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Over a range of binary catalyst compositions — 
Bi20^.3Mo0^ (a-phase), Bi20^.2Mo0^ (P -or Erman phase), 
Bi 20 |^.Mo 0 ^ (y - phase or koechlinite ) and the ternary catalysts 
based on bismuth molybdate (e.g. Bi 20 ^-Mo 0 j-Al 20 ^ ) , oxidative 
dehydrogenation of butenes has been reported to be first order 
in alkene and zero order with respect to oxygen over a wide 
range of temperature [2, 122]. Keizer [63] and Bakshi et.al. 
[12], however, have reported a deviation from zero order v/ith 
respect to oxygen over koeihlinite catalyst at temperatures 
in excess of 380°C and high conversions. Appearance of a 
dependence of rate expression on partial pressure of oxygen 
was explained by the blocking of active centers [l2] . 

Structure and properties of bismuth molybdates 
together with the mechanism of oxidative dehydrogenation 
reactions over these are discussed in detail in Sections 2.7 
through 2.9. 

Antimony Oxide Containing Catalysts . Many catalysts containing 
antimony oxide are selective in the oxidation of butenes to 
butadiene. However they are not as active and selective as 
bismuth molybdate catalysts. 

Of the antimony oxide -containing systems, iron (III) 
oxide -ant imoiQr oxide has been extensively studied by Boreskov 
and co-workers [27, 108, 110]. These catalyst were found to 
have a good selectivity to butadiene when they were present 
in the reduced state. 
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Catalysts based on antimony oxide-tin (IV) oxide 
system have also been fairly well investigated [106, 109, 12l]. 
These catalysts give good results when the antimony to tin 
ratio ranges from 1;16 to 4:1. Antimony oxide-tin oxide 
catalysts are not very sensitive to the olefin structure of 
the hydrocarbon being dehydrogenated and with 2 -butenes 90% 
selectivity at 63% conversion has been claimed. This was 
obtained with steam dilution, the ratio of butene : air: steam 
being 1:5^4 and a contact time of 4 sec at 440°G . 

KIIOiTICS MD BEACTIOU HBGHMISM OVER 
BISMUTH MOLYBDATE CATALYSTS 

The oxidative dehydrogenation of butenes over bismuth 
molybdate catalysts has been a subject of considerable study. 

A brief review is presented below. 

2.7 Kinetics of Dehydrogenation of n-Butene over Bismuth 
Molybdates . 

Over a range of binary catalyst compositions, 
Bi20^.3Mo0^, BigO^.SMoO^ (Erman phase ) , Bi20^.Mo0^ (koechlinite ) 
and ternary catalysts based on bismuth molybdate (e.g. 

, the oxidative dehydrogenation of butenes 
has been reported to be of the first order in butene, and 
zero order with respect to oxygen over a wide range of 
temperatures [2, 63]. These results indicate that the 
chemisorption of butenes on the catalyst surface is slow and the 
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chemisorption of oxygen is extremely fast. Detailed studies 
on the kinetics of dehydrogenation of butenes over bismuth 
molybdate catalyst have been carried out by Komarovskii et al. 

[693 . 

2.8. Mechanism of Oxidative Dehydrogenation . 

Two mechanisms were proposed simultaneously by Adams 
[l] and Sachtler and DeBoer [105] for the oxidative dehydro- 
genation of 1 -butene to 1, 3 -butadiene. The mechanism 
postulated by Adams involved initial adsorption of oxygen on 
the Catalyst surface prior to the occurrence of the reaction 
while Sachtler and DeBoer showed that the oxidation is a 
classical example of the Mars-van Krevelen mechanism [85] . 
Batist et.al. [20] later established the details of the 
Mars-van Krevelen mechanism for the catalytic oxidation of 
1 -butene to 1, 5 -butadiene over a bismuth molybdate catalyst 
as follows: 

C^Hq + Q + 0^" C^h:^ + OH" 

Mo^'^ + [Mo = (2.1) 

[Mo = + 0^" Mo^'^ + P + OH" + C^Hg (2.2) 

20H~ ^ 0^~ + P + HgO (2.3) 

^2 -> 20^" + 2Mo^'^ (fast) 


(2.4) 
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In reaction (2.1) and subsequent steps, Q denotes an anion 
r 1 5 + 

vacancy and L Mo = represents a Ti-allyl complex of 

Mo^ and the all 3 a. carbanion (C^Hy). The kinetics of the 
reaction were reported to be of the first order in 1 -butene 
and zero order in oxygen, leading to the conclusion that 
reduction steps (2.1) and (2.2) are rate determining and that 
reoxidation step (2.4), is comparatively fast. It is of 
interest that the reverse of reaction (2.1 ) leads to cis- and 
trans-2-butenesi.e . double-bond isomerisation, though the 
isomerisation reaction is slow compared to the oxidation of 
1 -butene to 1, 3 -butadiene. 

The validity of the Mars -van Krevelen mechanism was 
confirmed from the interesting experiments performed by 
Batist et al. [21 ] on the reduction and reoxidation of 
BigO^.MoO^ catalyst at 470°C. The yellow compound MoO-^ 

was reduced by 1 -butene to a blue black compound Bi^O^.MoO^ ^ 
(i.e. Bi^ + Mo^ ) which showed an X-ray diffraction pattern 
similar to that of the original catalyst, implying that the 
reduction did not affect the position of the metal atoms. The 
catalyst could not be reduced beyond a state formally described 
by Bi20^.Mo02 ^ using 1 -butene. However the reduction with 
hydrogen resulted in a compound containing bismuth metal and 
molybdenum (IV). The reduction of Bi20^.3MoO^, Bi^O^. 2MoO^ 
and Bi^O^.MoO^ with hydrogen studied by Beres et al. [24] 
has also indicated a rapid diffusion of oxygen through the 
crystal lattice. 
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Prom thexr adsorption studies of butenes, butadiene, 
water, oxygen etc. over some bismuth molybdates, Matasuura and 
Schuit [87] have established the existence of tv/o different 

surface sites. These sites were termed the A-sites and the 

2 _ 

B-sites. The "A-site" is an 0 -ion in a special position and 
IS associated with a Bismuth ion. The A-sites adsorbed 
butadiene slowly but strongly, and adsorbed water weakly and 
reversibly. Oxygen adsorption on the partially reduced catalyst 
was irreversible. The "B-site" is comprised of an anion 
vacancy on a Mo -ion sorrounded by corner-shared 0 -ions. 

The adsorption of butadiene on B-sites was found to obey a 
dual-site mechanism and the same was fast and reversible 
representing a weak adsorption. The configuration of the 
oxygen ions in the site is shown below. 

Og ( 1 ) ^ ^ layer 

0^ layer 

Og(3) layer 

Matsuura and Schuit [87] proposed the following steps for 
the formation of butadiene from butene. 
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In step (2.1), the adsorption of butene leads to the formation 
of a hydrogen atom and an allyl group (H-atom located on 
0^(1) and allyl group on 0^(2)). The dissociation of a second 
hydrogen atom and the formation of a butadiene weaicly attatched 
to 0g(2) and 0^(3) taJkies place in the second step. 
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At this point, there axe two alternate steps. In the first 
one the weaJcly adsorbed butadiene is desorbed and two hydrogen 
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The desorption of water and furan results in the reduction of 

the catalyst. The catalyst is reoxidised at an A-site either 

2 _ 

by diffusion of an 0 ion from the lattice or by the gas phase 
oxygen^ 

Oxidation of propylene to acrolein over Bi20^.3Mo0^, 

Bi20^.2Mo0^ and Bi20^.Mo0^ showed that activity of the catalyst 

was insensitive to the type of the molybdate [54] . Configuration 

of oxygen atoms around molybdenum atoms was different in the 

three molybdates and thus it was concluded that the Og sites 

should be Mo^"*" and not the oxygen associated with Mo^^. 

5 6 

Tracer studies using exchange of 0 isotope with partially 

reduced catalyst confiimied that the A-site is an oxygen ion 

associated with a bismuth ion [65, 129]. 

Attempts have been made in the past to relate the 

activity of the catalyst with the reactivity of the double 

bond type lattice oxygen Mo^"*" = 0 [8, 9? 54] . Mo^'*' = 0 bond 

strengths in bismuth molybdates are in the order Bi20^.MoO^'> 

BigO^ . 2MoO^ Bi^O^ . . However, in the butene oxidative 

6 + 

dehydrogenation, the relative bond strengths of Mo = 0 was 
of no consequencG [54]. 

2.9 Acidity and Activity of the Catalysts . 

Ai and co-workers [4-6] have studied the acldii;y 
of bismuth molybdate catalysts in relation to their activity 
for the oxidative dehydrogenation of butenes. Acidity of pure 
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BioO, was foijnd to be lower than that of MoO, or MoO^-P^O^, and 

that wxth increase in Bi^O^ content of the catalyst, the acidity 

rapidly increases at first, passes through a maximum at a Bi/Mo 

atomic ratio of 0.2, and then decreases to the value of pure 

Bi 20 ^ catalyst. Basicity of the catalyst was found to increase 

with the Bi/Mo ratio. The acidic sites were ascribed to Mo^^ 

5+ 

cation, since it was more electronegative compared to Bi . 

Anion vacancies proposed by Schuit and co-workers [I9j 87] were 

2 _ 

associated with surface lattice oxygen, 0 which were earlier 
confirmed by Wragg et al [130]. 

The activity of the bismuth molybdate catalyst for the 
oxidation of butenes to butadiene was found to be the best at a 
Bi/Mo ratio of 0. 5-3.0 (medium acidity). At lower Bi/Mo ratios 
(high acidity), acidic products like maleic anhydride were 
formed in preference to butadiene. 

2.10. Choice of Catalyst for the Present Work . 

The dehydrogenation of butane involves abstraction of 
hydrogen atoms from the parent molecule. Paraffinic hydrocarbons 
act upon the metal oxide catalysts more or less as electron 
donating agents; this tendency being very weak, in butane, its 
dehydrogenation activity will be governed by the activity of 
the surface sites. Therefore, an efficient catalyst will be 
the one which has either active Lewis acid-sites that can 
abstract hydrogen atoms from butane as hydride-ions, or 
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has suitable Bronsted basic sites which can abstract protons 
from butane molecule or both. Acidic sites consisting of metal 
cations are then necessary for adsorbing the remainder of the 
butane molecule. Thus the dehydrogenation of butane can 
proceed by stepwise abstraction of hydrogen over these centers. 

Molybdenum oxide on aluminium phosphate has been used 
in the past as a catalyst for the dehydrogenation of butane in 
the presence of oxygen [50] . Mo^”^ ions provide good active 

2 — 

centers for the adsorption of olefins. However, the oxygen, 0 
ions associated with the molybdenum tetrahedra are not mobile, 

and hence are not efficient in oxidising the hydrogen abstracted 

2 _ 

from the butane molecule, lattice oxygen, 0 associated with 
Bi has been proven to be very mobile, and can oxidise the 
abstracted hydrogen easily. 

The aim of the present work was to obtain dehydrogenation 
products like butenes and butadiene from butane using oxygen, as 
the hydrogen acceptor with as little oxygenated products 
(aldehydes, alcohols etc.) as possible. The molybdic oxide- 
aluminium phosphate is an acidic catalyst and hence may have a 
good selectivity for oxygenated products. 

Bismuth molybdate -aluminium phosphate catalyst was 
chosen for the present work from the above considerations. 

Bismuth oxide helps not only in reducing the acidity of the 
catalyst but also provides mobile lattice oxygen. This should 
make the bismuth molybdate -aluminixmi phosphate catalyst more 
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active and selective than the molybdic oxide -aluminium phosphate 
catalyst in the oxidative dehydrogenation of n-butane. 

BigO^.MoO^ (koechlinite or y -phase), 

Bi20^.2Mo0^ (Erman or P-phase), 

and BigO^ . 3MoO^ ( a -phase ) 

were chosen as the different phases of catalyst for the study. 
These three phases give catalysts with different Bi/Mo ratios. 
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CHAPTER 3 

MATERIALS , EQUIPIUT AND EXPERIMENTAL PROCEDURES 

3.1, Peed Cases , 

n -Butane . C.P, grade n-butane obtained from Mat he son Gas 
Company was used in all the experiments. Typical composition 
of C.P. grade n-butane is as follows: 

Minimum purity 99 vol % with small amounts of 

isobutane and 2-2, dimethylpropane 
as impurities. 

Sulphur content less than 0.005 wt % 

Moisture about 100 p.p.m. 

Oxygen . Industrial grade oxygen was supplied by Indian Oxygen 
Ltd., Kanpur. Oxygen gas had a minimum purity of 99 vol 

Nitrogen . Technical grade nitrogen (99.5% Pure) supplied by 
Indian Oxygen Ltd., Kanpur, was used in the present study. 

Air . An oil-free rotary compressor with teflon seal was used 
to compress air. The compressed air was passed over a bed of 
calcium chloride and activated charcoal in order to remove 
moisture, carbon dioxide and any other impurities. 



31 


PREPABATION OP CATALYSTS 


Effect of preparation of catalyst on the activity of 
bismuth molybdate catalyst has been investigated by Batist 
at al. [ 17 ] and G-rzybowska et al. [52]. Procedures suggested by 
Batist and co-workers [17, 70] have been used in the preparation 
of bismuth molybdates in the present study. Aluminium phosphate 
was prepared according to the procedure of Gaspar and Pasternak 
[ 50 ], Reagents used in the preparation of catalysts are as 
follows; 

Bismuth Ritrate . (Bi(N0^ )^ . 5 H 2 O). BDH Analar grade bismuth 
nitrate was used in most of the catalyst preparations. Bismuth 
nitrate was analysed by EBTA titration with xylenol orange 
indicator (see Appendix I) and was found to be 98.97^ pure. 

Bismuthyl Nitrato . (BiORO^ .H 2 O) . Bismuthyl nitrate was 
prepared in the laboratory as follows- 500 g of bismuth nitrate 
wrere placed in a 4-litre beaker to which 2.5 litre of distilled 
water were added and stirred vigorously for 5 minutes. The 
bismuthyl nitrate slurry v;as filtered and the unwashed cake was 
dried at 80 oQ [89] for 24- hours in an oven. Bismuthyl nitrate 
was analysed by EBTA titration with xylenol orange indicator 
(see Appendix I for details) and was found to contain 81.09% 
Bi 20 ^ by weight. 


Molvbdic Acid . (112^00^). Pure grade molybdic acid was 
obtained from Reachim, U.S.S.R. It was analysed using oxine 
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reagent and was found to contain 86.934^ MoO^ by weight. 

Details of the ano.lysis are given in Appendix I . 

Ammonium para Molybdate . ((lilH^)g MOr^O^^ . 415^0 ) . J.T. Baker 
Analysed Reagent. Analysis certificate indicated that it 
contains 82.5% MoO^ by weight. 

Aluminium Nitrat e. (A1(R0^ )^ . 9H2^ ^ * Sarabhai Chemicals G.R. 
grade and J.T. Baker Analyzed Reagent grade aluminium nitrate 
were used# 

Ortho-Bhosnhoric Acid . (H^PO^). G.R. grade o-phosphoric acid 
(85% weight), supplied by E. Merck (India) ltd. was 

used. 

Ammonia Solution . G.R. grade ammonia solution (Sarabhai 
Chemicals), containing 25 weight percent ammonia - sp. gr. 0.91 
was used in the present study. 

Nitric Acid . (HNO^). BDH Analar grade, and Pfeiaer Certified 
Reagent containing 69 wt % FNO^ were used. 

3.2. Bismuth Mo lybdates . 

1 . Y-Bismuth Molybdate (Bi oO^,^Mo0^2 . 

Procedure suggested by Konings et al. [7l] was used 
in the preparation of y-bismuth molybdate. 114.9403 g (0.4 mol) 
of bismuthyl nitrate wore added to 34.686 g (0.2 mol) of 
molybdic acid suspended in 2 1 of boiling water. A white 
slurry, with a pH of 2.5, was formed which slowly started 



55 


changing its colour after about half an hour. It was lemon 
yellow after 3 hours. The suspension was stirred vigorously 
for 20 hours under boiling conditions. The mass was then 
cooled and filtered. The filtered cake was washed with watei:* 
and then dried at 110°C for 24 hours in an oven. The filtrate 
from the unwashed cake was analysed and was found to contain 
very little bismuth molybdate in it. 

Dried bismuth molybdate cake was ground to a fine 
powder and calcined in a muffle furnace at 500 + 10°C. Chemical 
composition of the compound was verified by chemical analysis 
(see Appendix I). 

11 . 3-Bismuth Molybdate (Bi ^O ? • 2MoO ^ ) . 

0.5 mol of catalyst was prepared according to the 
method suggested by Batist et al. [1?]. 

490.075 g (1 mol) of bismuth nitrate were dissolved 

3 

in 180 cm of nitric acid and 2 litres of water. 

174.49 g of ammonium para molybdate were dissolved in 

3 

2 litres water and 50 cm of 0.91 sp. gr. ammonia solution. 

Ammonium molybdate solution was added to bismuth 
nitrate solution resulting in a white precipitate. Dilute 
ammonia solution was added to the slurry to adjust the initial 
pH of 1 to 6.5. It was then filtered and the residue was 
dried at 110°C for 24 hours to obtain a faint yellow coloured 
product. The powdered sample was calcined at 500 + 10°C for 

3 hours and then further kept at 450 + 10°C for 90 hours 
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during which time the colour of the sample turned yellow. The 
composition of the compound was checked by chemical analysis 
(see Appendix I). 

111. g-Bismuth Molybdate (Bi ^O ^ .3MoO ^ ) . 

0.5 mol of catalyst was prepared by using the procedure 

suggested by Batist et al. [17] . 

490,075 g (1 mol) of bismuth nitrate were dissolved 
3 

in 180 cm of nitric acid and 2 litres of water. 

261 .735 g ammonium para molybdate were dissolved in 

3 

2 litres of water and 50 cm of ammonia solution of specific 
gravity 0.91 . 

Ammonium molybdate was added dropwise to the bismuth 
nitrate solution and a white precipitate was obtained. The 
solution was kept well stirred and was heated to evaporate the 
water. Later it was dried in an oven at 110°C for 24 hours and 
calcined at 500 + 10°C for 4 hours. 

3 . 5 . Ortho-Aluminium Phosphate (o-AlP0_^ ). 

The procedure suggested by Gaspar and Pasternak [ 50] 
was used in the preparation of o-AlPO^. 

472.7 g (1.26 mol) of aluminium nitrate were 
dissolved in 3.8 litres of water. 145.3 g of ortho-phosphoric 
acid were then added and the contents were mixed. An ammonia 
solution of 0.91 specific gravity was added slowly to the well 
stirred aluminium nitrate -phosphoric acid solution till the pH 
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of the slurry reached 7. 0-7.1. The white coloured aluminium 
phosphate slurry was tnen filtered and the cake was washed with 
water and dried in an oven at 110°C for 48 hours. 

Aluminium phosphate was prepared in six batches. To 
avoid the batch to batch variations, aluminium phosphate 
prepared in different batches were mixed together, crushed gnd 
ground to a fine powder in a stainless steel ball mill with 
stainless steel balls. 

Aluminium phosphate thus prepared contains ammonium 
nitrate and some ammonium phosphate as impurities. Three 
samples of aluminium phosphate were ignited separately at 700°C 
in silica crucibles. The ignited residue was cooled in a 
desiccator and weighed. The percent o-AlPO^ in aluminium 
phosphate was computed. The samples contained 73.156 wt %. 
o-AlPO^ in the prepared aluminium phosphate. The aluminium 
phosphate was stored in the unignited condition and was ignited 
prior to its use. 

3.4. Molybdic Oxide -Aluminium Phosphate Catalyst . 

Molybdic oxide -aluminiiua phosphate catalysts with 
various MoO^/AlPO^ ratios were prepared according to the 
procedure described below for the preparation of a typical 
catalyst with a weight ratio of MoO^/AlPO^ of 0,1. 

41.676 g (0.25 mol) of aluminium phosphate were 
Ignited at 600 + 10°C for 5 hours in a muffle furnace. The 
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Ignited sample was ground to a fine powder and then diy mixed 
with 4.13694 g (0,025 mol) of molyhdic acid. Dry grinding was 

3 

continued for 15 minutes. About 60 cm of water were then 
added before wet grinding. Later it was dried in an oven at 
110°C for 12 hours. The dried catalyst was crushed t o 1“2 mm 
fraction and calcined at 600 + 10°C for 4 hours in a muffle 
furnace. The catalyst was cooled and stored in a desiccator. 

3-5. Bismuth Molybdate -Aluminium Phosphate Catalyst . 

Bismuth molybdate -aluminium phosphate catalysts were 

prepared for varied bismuth molybdate to aluminium phosphate 

ratios and different phases of bismuth molybdates ( a, p and y). 

Methods of preparation of these catalysts are identical. A 

typical preparation is described below for the case of a catalyst 

with BigO^ . 2Mo0^ : AlPO^ ratio of 0.05. 

41.676 g ( 0.25 mol) of aluminium phosphate were 

Ignited at 600 + 10°C for 3 hours in a muffle furnace. 

Ignited aluminium phosphate was ground to a fine powder and 

then mixed with 9.424 g (0.0125 mol) of Bi20^.2Mo0^. The 

3 

mixture was ground for 15 minutes. About 55 cm of water were 
then added and wet grinding was carried out. Later it was 
dried in an oven at 110°C for 12 hours. The dried catalyst 
was broken into 1-2 mm fraction and activated at 500 + 10°C for 
10 hours in a muffle furnace. The catalyst was cooled and 
stored in a desiccator. 
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EXPERIMENTAL SETUP MI> PROCEDURES 

3.6. Description of the Experimental Setup . 

A schematic diagram of the experimental setup is shown 
in Pig. 1. Three capillary flowmeters (5) were used to meter 
butane (1), nitrogen (2) and oxygen (3) separately. Manometers 
were provided to measure the differential pressure across the 
capillary as well as the absolute downstream pressure. Water 
was used as manometric fluid. Needle valves VI and V2 were 
provided to control the flow rates of gases. The valve VI was 
a fine control needle valve. Diying tubes containing calcium 
chloride (7)j were provided in the flow-loop to remove moisture 
if any, from the gases. A separate air line (4, 6, 7) was also 
provided for use. 

The reactor was made from a 15.1 mm I.D., 625 mm long 
l/2 inch schedule - 40 tube of type 304 stainless steel 
(Or 19%, Ni 9%, C 0.08% max:). Out of the total 625 nm length, 
280 mm long section was used for packing the reactor with 
catalyst and inerts (9, 10). The remaining 325 mm length was 
used for the preheating of the feed gases. A nitrogen-oxygon 
mixture was heated in the annular space formed by the reactor 
tube and a 12.5 mm O.D. stainless steel tube. Butane was 
heated separately inside the 12.5 mm O.D. tube, and was mixed 
with nitrogen— oxygen stream just before entering the catalyst 
bed to avoid any pre-catalyst bed reactions. Two chromel-alumel 
thermocouples, TCI and TC2 in 6 mm O.D. stainless steel sheaths. 
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Experimental setup. 



coaxial with the reactor tube, were provided for measuring the 
temperatures at the entry and middle of the catalyst bed. 

The reactor was placed in a 42 mm I.D. , vertical 
tubular furnace (8). The reactor temperature was controlled 
using "Alnor" type N-I5 temperature controller, 0-880°C range, 
supplied by Alnor Instrument Company. Thermocouple TC3 was 
used as the temperature sensing element for the controller. A 
cooler-condenser (11), fabricated from a 6 mm O.D. stainless steel 
tube was provided for rapid cooling of the reactor effluents. 

Trap (12), separated condensates if any, from the reactor 
product stream. A three-way stopcock (15), allowed the 
sampling of the product gases, which were passed over a column 
containing potassium hydroxide on asbestos for chromatographic 
analysis . 

Chromatographic Assembly . The chromatographic analysis of the 
reactor effluents, namely, ethylene, propylene, n-butane, 

1-, trans-, cis-2-butenes and butadiene required a column 
temperature of 21 °C. This necessitated a chromatograph with a 
refrigeration facility for the column oven since the average 
room temperature in the city of Kanpur is usually above the 
desired temperature of 21 °C* A chromatograph waa, therefore , 
constructed in the laboratory, the details of which are shown 
in Pig. 2. 

The chromatographic assembly consisted of an injector 
(11), a column (?), and a detector block (8) placed in a 
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water-bath (1). The injector had a rubber septum to facilitate 
the injection of gases to be analysed into the column. Along 
with the carrier gas, these gases were pasjsed into the column 
in which they were separated before entering the thermal 
conductivity detector (8). 

The chromatographic column and the detector block were 
placed in an insulated water bath, the teii5)erature of which was 
controlled within an accuracy of + 0. 1 °C by circulating chilled 
water (2, 3, 4). The chilled water circulation pump (3) was 
connected to a Jumo mercury contact thermometer (6) via a 15 
ampere contactor. 

The detector block was essentially the same as the 
one used in Varian Chromatograph Model Aerograph 90-P3. A 
matched quad of Gcw-Jilac l/7X-type 9225 tungsten-rhenium filaments 
with 12 inch long kovar leads were fitted in the detector block. 
The detector block assembly was placed in a 150 mm x 150 mm x 
220 mm aluminium container and was insulated with glass-wool. 

A Gow-Mac power supply (13), Model 9999-0 i 1 , 220 Y, 

50 Hz was used to supply current to the detector filaments. 

The ripple of the power supply was measured to be 0.7 mV. A 
Sargent recorder (14), Model SE-50, fitted with a Model 204 
Disc chart integrator, was used to record the chromatographic 
signals. High purity nitrogen (oxygen content less than 
0,001 vol %) was used as the carrier gas. 
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3.7. Process Measurements . 

Chromel-alumel thermocouples protected in. stainless 
steel sheaths were used to measure temperatures of the reactor 
and the furnace. The standardization procedure for the thermo- 
couples IS discussed in Appendix II. A Honeywell, Model 2720i 
portable calibrating potentiometer test set was used to measure 
the voltage across the thermocouple junction. The accuracy 
of the potentiometer was + 0.01 mV. 

Butane, oxygen and nitrogen flow rates were measured 
with capillary flow meters. Air flow rate was measured with a 
rotameter. The gas flow rates were calibrated at a down stream 
pressure of 50 cm water using a soap bubble flow meter of 

3 

200 cm volume . 

3.8. Experimental Procedure . 

The reactor tube was assembled and the threaded joints 
which were in contact v/ith high temperatures were sealed with a 
high-alumina/sodium silicate cement. The refractory cement was 
allowed to set for at least 6 hours. 

3 

The reactor tube was later packed with about 40 cm 

3 3 

of catalyst (1-2 mm), sandwiched between 10 cm"^ and 5 cm of 
porcelern packing (inert) on the exit and inlet sections 
respectively. Air was passed through the reactor at a flow 
rate of about 200 cm^ per minute and the reactor furnace was 
switched on. When the required temperature reached, nitrogen 
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and oxygen were supplied to the reactor and subsequently the 
air supply was stopped. 

When the catalyst bed temperature was steady, butane 
was introduced into the reactor. The flow rates of butane, 
oxygen and nitrogen were adjusted to the required values using 
the valves VI and V2 (Fig. l). When the measured temperatures 
at the inlet and the middle of the catalyst bed registered 
steady values, it was assumed that steady state has been 
achieved. Two reactor effluent samples were drawn out 
separately, one of which was taken into an Orsat apparatus to 
estimate carbon dioxide, oxygen and carbon monoxide and the 
other was used for chromatographic analysis. The latter sample 
was passed over a bed containing potassium hydroxide on 
asbestos before analysis in order to remove the carbon dioxide 

from the product stream which interfere with the analysis (see 

*5 

Section 3.9). About 10-15 cm of this gas sample was injected 
into the gas chromatograph for analysis. The rate of formation 
of the condensate, which contained small quantities of 
oxygenated compounds, was negligible. 

When the experimental conditions were altered for 
the next experiments, care was taken to change the partial 
pressures of butane and oxygen very gradually. A sudden change 
in the partial pressure of butane was found to increase the 
temperature of the catalyst bed catastrophically. 



44 


After the desired experiments were completed, air was 
again introduced into the reactor and the nitrogen and oxygen 
streams were shut-off j Butane partial pressure was tapered 
off to avoid sudden increases in the catalyst bed temperature. 

The oxide coating foimed on the reactor wall was 
frequently removed by pickling the reactor tube in an acid bath 
containing 3% HCl, Jfo H^SO^ and 94% water, at a temperature of 
about 70°C [58]. The reactor tube was thoroughly washed with 
distilled water after such a treatment and dried before placing 
in situ. 


AJAIYTIGAL METHODS 


3.9. Methods of Analysis . 

Chromatographic as well as Orsat analyses were 
employed to establish the composition of the reactor effluents. 

The Orsat analysis of the product gases was carried 
out to determine the percentage of carbon dioxide, oxygen and 
carbon monoxide in the product stream. Details of the reagents 
used in the Orsat analysis are given in Appendix III. 

The chromatograph used for the analysis of the 
hydrocarbon mixture has been described in Section 3.6. The 
chromatographic column was a coil of copper tubings 15.4 m 
length and 5.2 mm I.D., packed with 30/44 mesh Chromosorb P 
(obtained from Chromatography and Instruments Co., Baroda, 
Gujarat, India) coated with propylene carbonate (35% by weight). 
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Prosynth grade propylene carbonate was obtained from Riedel, 
Germany. Die hlorome thane (G.R. grade, E. Merck, Germany) was 
used as a solvent to deposit propylene carbonate on Ghromosorb P. 

The above column could completely separate methane, 
ethane, ethylene, propane, propylene, n-butane, 1 -butene, 
trans -2 -butene, cis-2-butene and 1, 3-butadiene. However, the 
carbon dioxide peak interfer with the propylene peak and 
methane, oxygen and carbon monoxide emerged as a single peak 
which warranted a separate method of analysis for estimating 
these components. Orsat analysis as described above was, 
therefore, used in addition to the chromatographic analysis for 
determining CO, CO 2 and O 2 . 

The separating power of the column depends upon a 
number of factors [47, 64 ]. Experiments were conducted to 
find the best conditions. The following operating conditions 
were arrived at based on minimum H.E.T.P. (height equivalent 
of a theoretical plate) for n-butane separation: 

Column temperature : 21°C 

•Z 

Carrier gas flow rate : 36 emv minute 

Pig. 3 shows a typical chromatogram. 

3.10. Calibration . 

In order to interpret the chromatogram in quantitative 
terms, the thermal conductivity detector of the chromatograph 
was calibrated using pure gas samples. The method of absolute 
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standardisation was used to calibrate the detector response 
[ 47 , 64, 83]* Ibe volume of gas handled being large, the 
errors involved in the accurate measurement of gas volumes was 
small . 

IXinng calibration, pure gas samples (O.P. grade), 
obtained from Matheson Gas Company, were injected into the 
column using 100 ^1, 250 pi and 500 pi Hamilton gas-tight 
microliter syringes. A Sargent Recorder (fitted with a disc 
integrator) with a 1 mV range was used to record the 
chromatograms* The retention times and the response factors 
were calculated by injecting various volumes of the same gas 
and averaging the peak area per unit volume of the gas injected. 
Reproducibility of the peak area was checked by injecting 
identical volumes of the gas at different times. The response 
factors thus calculated are given in Table 5* The specific 
volume data of the hydrocarbon gases, used in the calculation 
of the response factors, were taken from Matheson Gas Data 
Book [86 ]. 

The retention times of hydrocarbons changed with the 
operation of the column over a period of several months. Ho 
efforts were made either to change or repack the column during 
the investigation since the separation was always complete. 
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Table 5* 

EETENTIOir TIMES AND 
POE HYDROCAEBONS 

RESPONSE FACTORS 

Component 

Retention Time 
minutes 

Response Factor x 10 
mol/integrator count 

Ethane 

10.6 

0.1655 

Ethylene 

11.5 

0,202 

Propane 

13.6 

0,1626 

Propylene 

17.7 

0.216 

n-Butane 

20.0 

0.2297 

1 -Butene 

30.3 

0.2172 

t r ans -2 -But e ne 

36.0 

0.2424 

CIS -2 -Butene 

42.0 

0.2464 

1, 5 -Butadiene 

60.0 

0.2483 
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3<11< Process Calculations . 

Quantities like conversion^ yield etc. were calculated 
by a material balance of carbon# hydrogen and oxygen^ The 
details of this procedure is presented in Appendix IV . 
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CHAPTER 4 

CHARACTERISATION OP BISMUTH MOLYBDATE -ALUMINIUM 
PHOSPHATE CATALYSTS 


4.1 4 Introduction 4 

The activity of a catalyst is related to its solid 
state properties like surface area» crystal structure, phase 
composition and thermal stability. Such properties aid in 
understanding the catalytic behaviour of the solid and 
elucidating the reaction mechanism. 

The catalysts used in present study were characterised 
in terms of chemical compositions, characteristic X-ray 
diffraction reflections and thermal stability. The various 
properties studied and their implications are listed below, 

4.2. Chemical Gomnosition . 

Bismuth molybdates used in the present investigations, 
were prepared in three stoichiometric compositions, with varying 
bismuth to molybdenum ratios. The y-i^base (BigO^.MoO^) has 
been obtained by slurry reaction while the p-phase (BigO^ . 2MoO^ ) 
was prepared by the precipitation method. The a-phase 
(Bi 20 ^. 3 MoO^ ) has been prepared by evaporating the slurry 
obtained from the precipitation reaction between bismuth 
nitrate and ammonium molybdate and sinee the reagents involved 
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are of known purity > the product has not been analysed for 
its chemical compositioni 

The bismuth molybdate suspensions of y- and P-phases 
were separated from their mother liquor and the chemical 
composition of the same in the diy fonn confirmed that the 
compositions are essentially that of the Y- and p -phases 
(Bi20^.Mo0^ and Bi20^t2Mo0^) , This confirmation at this stage 
was necessary since any possible losses in the filtrate would 
haYe considerably affected the chemical compositions of the 
bismuth molybdates prepared. The molybdenum contents in these 
two bismuth molybdates were analysed using oxine reagent. 

BDTA was used as a masking agent to prevent the interference by 
bismuth ion [43# 97]* Experimental details for the same are 
indicated in Appendix I. The P- and Y“T bases contained. 38.35% 
MoO^ and 23 *8 9% MoO^ (by weight) respectively and these closely 
correspond to the calculated percentages. 

4.3. Bulk Density . 

Various catalysts were prepared using bismuth 

molybdates (i.e. the a-, P- and Y-phases) in different 

proportions with ortho-aluminium phosphate. On characterisation 

since the a-phase has been found to be unsuitable as a catalyst 

(on the basis of X-ray diffraction patterns and infrared spectra 

detailed in Section 4.5 and 4.6 respectively), it was not used 

in the present study. As such only the catalysts involving 

* - ■ i- > 

Abb, No. ^ 
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the P- and y-phases of the bismuth molybdate were adopted. 

Typical bulk density values for the catalysts 
containing P- and Y -phases are presented in Table 6. It was 
observed that the bulk density of the prepared ortho- aluminium 

’Z 

phosphate was around 0.27 g/cm"^, which is very low compared to 

•2 

the reported value around 2.57 g/cm^ of the crystalline ortho- 
aliiminium phosphate [127]. The aluminium phosphate in the 
present case was amorphous as evidenced by the X-ray diffraction 
data. 

The bulk densities of Bi^O^.MoO^-AlPO^ and Bi20^.2MoO^- 
AlPO^ catalysts were foiond to increase with the increase in the 
proportions of the active components. This is expected since 
the densities of the constituent bismuth molybdates were high 
(3-6 g/cm^). The bulk densities of both types of bismuth 
molybdates were comparable at equal ratios of bismuth molybdate 
to aluminium phosphate. 

4.4. Surface Area Measurements . 

The surface area determinations were made using liquid 
nitrogen adsorption isotherms (adopting a standard value of 
16.2 A^ for nitrogen molecule) with a static apparatus. The 
typical results for the catalysts involved are presented in 
Table 7. 

It has been observed that the addition of bismuth 
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Table 6^ BUIE DENSITIES OF BIStlUTH MOLYBDATE -AlUMIFIUM 
PHOSPHATE catalysts 


S.No. 

I f 

t t 

Catalyst Composition { 

*2 

Bulk Density g/cm^ 


mol/100 mol o-AlPO^ { 

1 

( 

\ t 

BigO^ »MoO^ 

BigO^ * 2Mo0j 

1 

0.0 (o-AlPO^) 

0.26 

0.26 

2 

1 .0 

- 

0.25 

3 

3.0 

- 

0.27 

4 

5.0 

- 

0.29 

5 

7.0 

0.32 

0.32 

6 

9.0 

- 

0.34 

7 

15.0 

0.41 

0.40 

8 

21 .0 

- 

0.47 

9 

23.0 

0*46 



10 


25.0 


0.49 


0.49 
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Table 7. SUKPAOE ASEA VALUES OP THE CATALYST SAMPLES 


S.No, 

1 

I 

1 Catalyst 

1 

f 

\ 

f 

1 - - - - - - - 

1 

1 

1 Specific Surface 

1 Area 

j m^/g 

f 

1 

AlPO^) calcined 

49»2 

2 

BXgO^.MoO^, calcined 

9.3 

3 

BigO^ . 2MoO^ , calcined 

3.8 

4 

0,1 MoO^ ; 1 AlPO^, calcined 

37.8 

5 

0.07(Bi20^ ,MoO^ ) ; 1 AlPO^, calcined 

21.1 

6 

0.05(Bi20^.2MoO^) ; 1 AlPO^, 
dried 1 uncalcined 

83.0 

7 

O.ObCBi^O^ . 2MoO^ ) : 1 AlPO^ calcined 

11.1 
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molybdate or molybdxc oxide to aluminium phosphate reduces the 
specific surface areas of the catalysts. The uncalcined 
Bi20^<2Mo0^-AlP0^ catalyst (with around 12% of water as indicated 
by the TG-A data) has a surface area higher than that of calcined 
Bi20^.2Mo02-AlPO^ as also the calcined AlPO^ base. 

4.5f X-Ray Diffraction Studies . 

X-ray diffraction studies enable the determination of 
the nature and the extent of the individual constituents present 
in any material^ while the chemical analyses indicate the overall 
chemical composition taking into consideration all the consti- 
tuents within the material. In the present investigations this 
technique was adopted to ascertain the crystallinity of the 
starting materials prepared as also to observe any possible 
changes in the starting materials on calcination to around 500 °0 . 

The X-ray diffraction analyses have been carried out 
on General Electric XRD-5 diffractometer using Ei-filtered Cuk^j 
radiation. The samples were scanned between 10° and 60° (20) on 
a scale of 2° (29) to an inch. The 'd* spacings were calculated 
for various reflections obtained in each case using the standard 
Bragg *s equation. 

X-Ray Diffraction Results ; 

1. Ortho-Aluminium Phosnhate (o-AlPO^ ). 

The characteristic diffraction pattern (Pig. 4) for 
this material has indicated only a very broad hump in the range 








57 


of 12° to 18° (20), No other* reflections could be noticed. 

Several peaks reported in the literature [75] for a crystalline 
phase of this material are absent here# thus indicative of its 
amorphous nature. 

11 . Bismuth Molybdates . 

In the present case, diffraction patterns of the a-, 

and Y-i)haseB were studied. Typical diffraction pattern for 

BigO^.MoO^ sample is indicated in Pig. 4 with the corresponding 

'd' spacings and the peak intensities are presented in Table 8. 

The ^d' spacings in the present case are in close agreement with 

the values reported by Prondel [49] for the bismuth molybdate 

(mineral koechlinite) with composition BigO^.MoO^ (Table 8). On 

the basis of the 'd' spacings for 200, 002 and 060 reflections, 

0 

the cell parameters have been computed to be a = b = 5.46 A and 
d 

b = 16.08 A. These are also in agreement with the reported 
values [45# 49]. 

The results of the X-ray measurements of Bi20^.3MoO^ 
are shown in Pig. 5. The corresponding *d' spaeings are given 
in Table 9. The X-ray diffraction data of the prepared Bi^O^. 
3MoO^ closely resemble those given by Erman et.al. [45] (Table 9). 
Thus the Bi20^.3Mo0^ prepared in the laboratory has the monoclinic 
stiueture as reported by Erman et.ai. 

In Pig. 6 the characteristie X-ray diffraction 
reflections for the prepared -bismuth molybdate (3120^.2100^) 
are indicated. The ‘d' spacings for the reflections (Table 10) 
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Table 8< X-RAY DIPRRACTION RATA FOR Bi«0,,MoO,. 

dp P 

*d' SPACIRG-S MB PEAK INTENSITIES (I)< 


Present Work 

Prondel t49] 


H 

* 

f 

0 1 


f~ — 

r 

f 


H 

I 

d, A ! 

I 

! hkl 

HHHIi 

H 



i__ -1 


1 

f 

» , - - . 

3.26 


3 




3.15 


100 

3.13 

100 

131 

2.87 


4 




2.73 


74 

2.73 

60 

200, 002 

2.69 


79 

2.68 

50 

060 

2.60 


11 

2.59 

10 

220, 022 

2.478 


25 

2.473 

30 

151 

2.265 


8 

2.263 

20 

240, 042 

2.006 


2 




1.938 


49 

1.936 

60 

202 

1.926 


86 

1 .918 

80 

260, 062 

1 .888 


13 

1.879 

10 

222 

1 .652 


69 

1.647 

90 

331, 133 

1 .632 


70 

1.628 

70 

191, 280 

1.575 


32 

1.570 

60 

260 


* 


Arbitrary units indicating the line intensities. 
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Table 9* X-RAY DIFFRAGTIOF DATA OF Bi20^, 
SPACIITGS AUD PEAK INTENSITIES (I) 

3MoO^ , 

» 

Present Work 

1 

f 

1 


Erman et 

.al. [46] 

jr. 

0 i 

d, A 1 

t 

r 

** 

I 

! 0 
; d, A 

I 

- - 

t 

♦ 

I 

f 

* 

I 

( 

1 

t 

1 

... Vr 

hkl 

4.87 

26 

4.83 


w. 

012 


3.59 

16 

3.60 


w* 

102 


3.41 

8 

3.42 


w» 

210, 

112 

3.18 

138 

3.18 


s » 

032 


3.05 

100 

3.05 


s. 

131 


2.87 

59 

2.84 


s * 

204, 

133 

2.481 

28 

2.477 


w* 

024 


2.243 

24 

2.248 


w. 

240, 

034 

2.001 

28 

2.002 


med* 

330, 

044, 106 

1.914 

20 

1.914 


w. 

061 , 

035, 226, 404 

1.877 

31 

1 .878 


med* 

160, 

316, 321 

1 .806 

11 

1 .805 


w* 

320 


1 .760 

9 

1.764 


w* 

400, 

045 

1 .713 

23 

1.715 


med. w. 

336 


1 .692 

26 

1.691 


med. w. 

435, 

431 

1 .592 

14 

1.589 


w. 

064 


1 .554 

14 

1.550 


V. w. 

337 


** Arbitrary units 

indicating the 

line intensities. 

* s. = strong; med. = medi\m; 

med. w. = medium weak; 

w. 

= weak; 

v.w. = very 

weak. 
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Table 10. X-RAI DII’I'RACTION DATA OD BigO^.aMoOj. 
'd' SPACniGS ADD PEAK INTENSITIES (I). 


Present Work 

f 

t 

1 

Batist et#al.^ [17] 


t 

0 t 

d, A I 

1 

f 

** 

I 

t 

t 0 

! d, A 

1 

f ^ r 

f 

1 

1 

f 

t 

I i 

f 



* 

hkl 

8.04 

18 

8.01 


20 


4.90 

33 

4.92 


70 


3.55 

78 

3.44 


5 


3.20 

362 

3.20 


600 

301 

3.15 

402 

3.15 


400 


3.05 

136 

3.06 


40 


2.99 

44 

2.98 


20 


2.88 

42 

2.87 


20 


2.80 

107 

2.81 


100 

321, 330 

2.75 

96 

2.73 


30 


2.70 

100 




002 

2.491 

42 

2.486 


12 


4.987 

44 

1.986 


65 

402* 600 

1 .941 

133 

1.945 


100 


1.924 

122 

1.928 


30 


1 .884 

42 

1.883 


20 


1.687 

80 

1.687 


100 


1.653 

98 

1.650 


30 

303 

1 .634 

107 

1.641 


95 


1.598 

42 

1 .592 


65 

602 

** Arbitrary units 

indicating the 

line 

intensities. 


* Prom the 

data of 

' Erman et.al. 

[46]. 
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are compared with the values reported by Batist et.ali fl?]* ^^he 
prepared Bi 20 ^i 2 Mo 0 ^ corresponds to p -bismuth molybdate as 
identified from the characteristic reflections at 5-95^ 4«90< 
2.80, 1.987, 1.941, 1.654 and 1.598 A (2© = 14,85°< 18*1°, 32.0°, 
45»8°, 46,85°, 56.3° and 57.7°). However the presence, to a 

restricted extent, of tho a-phase (as identified by its character- 

0 

istic line at 3.05 A at a 20 value of 29.2°) and y-phase (as 
inferred from its characteristic lines at 8.04, 1.924 and 
1,653 A at 2© values of 10.94°, 47.25° and 55.6°) has also been 
evidenced in the X-ray diffraction patterns. The coexistenco of 
these two phases in minor proportions with the p-phase is due to 
the partial decomposition of the latter at room temperature 
according to the following reaction 

BigO^ .2MoO^ ( P) — '>« BigO^.MoO^ (Y ) + Bi^O^ . 3MoO^ ( a ) 

(4.1) 

111 . Ei ^O ? .MoO^-AlPO^ Oatalyst . 

The X-ray diffraction patterns for the unealcined and 
the calcined Bi^O^.MoO^-AlPO^ catalysts are shown in Big. 4. 
Typical 'd^ spacings are presented in Table 11. As can be 
evidenced from the Big. 4, the unealcined and calcined specimens 
of the BigO^.MoO^-AlPO^ catalysts have similar 'd* spacings 
indicating that they have not undergone any noticeable 
modification in crystal structure. However, an interesting 
feature is apparent in the diffraction patterns. A comparison 



64 


Table 11* 

X-ray 

DIPPRAOTIOR DATA OP Bi^O^ .MoO^-AlPO^ 

GATADYST* 


BigO^.MoO^/AlPO^ = 

0,25 mol/molf 




SPACINGS AND PEAK 

IKTEKSITIES (l). 


T 

Uncalcined | 

. I 

Calcined 


0 — 

d, A 

" f 

t 

r 

1 .... 

5F 1 

I 1 

. ^rrr ^ ,, .,rT.,.r. 

d, A 1 

. - ^ . .. . 1 _ 

' * 

I 

8.04 


8 

8.05 

5 

3.26 


4 



3.15 


100 

3.15 

100 

2.88 


2 



2»74 


42 

2.74 

36 

2.68 


34 

2.69 

29 

2.60 


6 

2.60 

50 

2.50 


13 

2.49 

10 

2.27 


4 

2.276 

4 

1.941 


28 

1.939 

24 

1.926 


50 

1.927 

41 

1 .892 


7 

1.891 

6 

1.653 


38 

1.653 

38 

1.635 


33 

1.633 

30 

1.578 


17 

1.578 

17 


* Arbitrary unxts indicating the line intensities. 
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of certain reflections for BigO^.MoO^-AlPO^ with those of pure 

BigO^^MoO^ specimens indicated modifications in the relative 

intensities (Pig* 4). Por example j the intensities for the 

0 

reflections at 2»73 and 1.652 A corresponding to 29 values of 

32.75° and 55.65° have increased compared to their doublet peaks 

0 

at 2.69 and 1.632 A corresponding to 2© values of 33*5° and 56.35°. 
Since the reflections at 2.73 and 1.652 A are for 002 and 331 
planes respectively » it nay be inferred in the present case that 
a tendency for preferred orientations of the planes in these 
directions for Bi20^»Mo0^ exists in the presence of AlPO^. This 
appears to be also the case for Bi 20 ^. 2 MoO^ in the presence of 
AlPO^ (discussed below), 

IV. Bi ^O^; «3 Mo 0 ^ -A1P0^ Catalyst . 

The X-ray diffraction pattern for this catalyst is 
shown in Pig. 5 and the data are presented in Table 12. It can 
be seen from the figure that the relative peak shapes are not 
affected significantly. The peaks of the calcined and the 
unoaleined catalysts appeared in the same 2© ranges. Thus the 
presence of AlPO^ or calcination of Bi 20 ^. 3 Mo 0 ^ together with 
AlPO^ did not bring any modifications. 

V. Bi 20 ^.2MoO ^ -AlPO^ Catalyst . 

The X-ray diffraction patterns of the unealeined and 
the calcined BI 2 OJ .2MoO^-AlPO^ catalysts are shown in Pig. 6 and 
the data are presented in Table 13. A comparison of peaks in 
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A 


Table 12. 

X-RAY DIRRRAGTIOU DATA 

Bi20^.3Mo0^/AlP0^ = 

'd' SP AGINGS and PEAK 

OP Bi^O_.3MoO,-AlPO, 
2 3 3 4 

: 0.25 mol/mol. 

INTENSITIES (I). 

CATAIYST. 

Uncalcined i 

. - ... .... ... f 

Calcined 


d(A) 

1 f 

t ^ ! 

1 I ! 

- - 1- t _ 


* 

I 

4.88 

28 

4.89 

26 

3.59 

12 

3.60 

13 

3.18 

143 

3.19 

127 

3.05 

100 

3.06 

100 

2.87 

45 

2.88 

41 

2.481 

28 

2.483 

25 

2.244 

18 

2.245 

21 

2.003 

39 

2.003 

30' 

1 .916 

35 

1.918 

24 

1 .879 

34 

1 .880 

27 

1 .804 

15 

1 .807 

10 • 

1 .762 

11 

1 .764 

9 

1 .714 

21 

1.716 

18 

1 .694 

29 

1.695 

25 

1.592 

18 

1.593 

17 

1.556 

11 

1.556 

9 


* Arbitrary units indicating the line intensities. 









67 


Table 13. X-RAY DIPREAGTION DATA OP BigO^ • 2Mo0^-AlP0_^ CATALYST. 

Bi20^.2Mo0^/AlP0^ = 0.05 mol/mol. 



*d' 

SPACINGS 

AND 

PEAK INTENSITIES 

(I). 

Uncalcined 


“1 

t 

f 

Calcined 


0 

d, A 

! 

i 

♦ 

i 

* 

I 

f 

f 

0 i 

d, A \ 

• 



I 

7.97 


16 


8.19 

20 

4.90 


26 


4.95 

12 

3.55 


95 


3.53 

48 

3.20 


347 


3.18 

375 

3.15 


429 




3.06 


155 


3.06 

59 

2.99 


50 


2.90 

12 

2.87 


53 




2.80 


87 


2.81 

23 

2.75 


113 


2.77 

124 

2.70 


100 


2.75 

100 

2.488 


37 


2.501 

39 

1.991 


42 


2.008 

10 

1.943 


129 


1.949 

88 

1.926 


124 


1.934 

135 

1.884 


55 


1 .890 

27 

1 .687 


76 


1.691 

14 

1.653 


100 


1.659 

118 

1 .635 


96 


1.638 

98 

1.595 


21 


1.600 

9 

1.576 


42 


1.580 

47 


* Arbitrary units indicating the line intensities. 





68 


‘ 2MoO^-A 1PO^ Catalyst does not show any changes in the 
positions of the peaks for the two substances except that some of 
the relative peak intensities have increased slightly* This is 
indicative of the preferred orientation of certain crystal planes 
of Bi20^.2MoO^ in the presence of AlPO^ which has already been 
explained earlier. This feature was observed in the case of 
Bi^O^ .MoO^-AlPO^ catalyst also (see Section 4.5, v). 

In addition, a comparison of the intensities for the 
peaks of calcined Bi^O^. 2Mo0^-AlP0^ catalyst with the correspon- 
ding peak intensities for the unoalcined Bi20^.2MoO^-AlPO^ 
catalyst indicates certain structural changes pertaining to 
Bi20^.2MoO^ compound. The reduction in the relative intensities 

in the peaks of a-bismuth molybdate at 'd' spacings of 3 •53, 

0 

3.06 (most characteristic) and 2.008 A at 20 values of 25.1°, 
29.2° and 45.15° respectively in the calcined catalyst reveals 
the reduced presence of a -bismuth molybdate in the same. It may 
be inferred that the calcination of Bi20^.2MoO^ sample together 
with AlPO^ results in inhibiting its decomposition into a- and 
Y -phases. 

4.6. Infrared Measurements . 

Having established the overall structure of the bismuth 
molybdates and bismuth molybdate -aluminium phosphate catalysts 
through X-ray diffraction studies, their infrared spectra were 
obtained in order to understand the nature of the bond between 
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molybdenim and oxygen atons and the influence of aluminiun 
phosphate on it , The X-ray diffraction pattern was of no use in 
establishing the structure of aluniniun phosphate due to its 
amorphous nature. Accordingly, an attempt has been made to 
determine its structure using infrared technique. 

The infrared spectra were recorded on a Perkin-Elmer 
Model 521 spectrometer in the range 4000-250 cm . The samples 
were prepared in the form of KBr discs. 

Results of Infrared Analyses ; 

1 . Ortho-Aluminium Phosphate (-AlPO^ ). 

The ir spectrum in the range 250-1200 cm~"* for o-AlPO^ 
(Pig. 7d) exhibited two broad and intense bands at 1070 and 500 

.I 

cm . The former band has been ascribed to the triply degenerate 
stretching vibrations mode, Rakamoto [90b] ) of tetrahedral 
(PO^)"^". The band at 500 cm was attributed to the triply 
degenerate bending vibrations mode, Rakamoto [90b] ) of 
(PO^)” tetrahedra. Thus in the prepared aluminium phosphate, the 
phosphorous atom is surrounded by a tetrahedron of oxygen atoms. 
It IS also known that in the ortho-aluminium phosphates, the 
network of aluminium and phosphorous atoms is analogous to a 
Si(SlO^) network, with both aluminium and phosphorous surrounded 
by a tetrahedron of oxygen atoms [127] . In addition, two bands 

— "I 

(not indicated in the figure) around 1614 cm (intense) and 
3410 cm”"' (wide and intense) were ascribed to 0-H in-plane- 
bending vibrations and —OH stretching vibrations respectively. 
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^ Wav« mifbP" ^ * I *« 

® (o) ^McO% im catalyst 

(c> Catalyst CM55.mol Bl203St«03/ »««* o-AlP04,W> 
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The broad band observed in the range of 300-800 cn"’^ also contains 
the absorption band due to 0-H out-of -plane bending vibrations^ 

11 Bismuth Molybdates . 

In the bismuth molybdates, the molybdenum atoms are in 
octahedral configuration with comer-sharing or edge-sharing of 
the oxygen atoms f^S] . Hence, at the surface we can have either 
a terminal oxygen 0^, which is bonded to only one molybdenum atom 
or bridging oxygen that xs bonded to two molybdenum atoms. These 
bonds are represented as Mc^-O^ and Mo-0^ respectively. 

Mitchell and Trifiro [90a] have assigned infrared frequencies to 
the various vibrational modes of the octahedral molybdenum 
identified by Clark and Doyle [37]* The ir spectra of various 
bismuth molybdates prepared in the present study are interpreted 
below. 

The ir spectrum of Bi20^.MoO^ indicated in Pig. 7a 

exhibits bands at 950 (weak), 940 (weak), 910 (weak), 850 

(sharp + shoulder at 870), 800 (sharp), 740 (broad and intense 

with a shoulder at 780), 570 (intense with shoulders at 600 and 
— 1 

520) cm” . The bands observed at 740 and 850 em~ were assigned 
to Mo-0^ and Mo-0^ stretching vibrations. A prominent 
shoulder band of Mc-0^ at 800 cm is indicative of the 
tetrahedral molybdenum. 

The spectra for Bi20^,3Mo0^ and Bi20^.2Mo0^ are 
presented in Pigs. 8a and 8b. The band pattern in the region 
600-900 cm”^ has undergone modification in the present cases 
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compared to the patterns observed for Bi20^.Mo0^ sample. Intense 

bands are observed 850, 725 » 570 and 460 cm”\ The bands 

occurring at 740 and 850 cm are relegated to the stretching 

vibrations of Mo-0^ and Mo-O^ bonds. The presence of 

edge-shared molybdenum octahedra in the two bismuth molybdates is 

-1 

evidenced from the shaip bands observed in the range 900-950 cm , 

-1 

A weak band observed around 995 cm in the case of BigO^.^MoO^ 
was due to the Mo—0^ stretching vibrations of free MoO^. 

The broad bands noticed in all the bismuth molybdates 

-"i 

in the range 400-500 cm are attributed to the absorption by 

Bi20^. 

iii. Bi ^O ^.MoO ^ -AlPO^^ Catalysts . 

Bigs. 7b and 7c show the ir spectra of the two catalysts 
with compositions 0.25 mol BigO^.MoO^/mol AlPO^ and 0.05 
BigO^ .MoO^/mol AlPO^ respectively. The ir spectra for these 
catalysts exhibit the general features of the Bi20^.MoO^ compound 
(Pig. 7a). The bands due to adsorbed water as in the case of 
o-AlPO^ were also observed (not shown in the figure) aroimd 1615 
and 3410 em*"^ in the present ease. 

iv. Bi 20 ->.2Mo0v-AlP0_ ^ Oatalvsts . 

The ir spectra of these catalysts are shown in Pigs. 

8e and 8d, Both spectra exhibit the general features observed 
in ir speetmm of Bi20j.2MoO^ (Pig. 8a). The absorption bands 
due to adsorbed moisture were also indicated at 1615 cm and 
3410 cm"^. 



74 


4*7* Thearmal Characteristics * 

The differential thermal analyses (DTA), the thermo- 
gravimetric analyses (TGA) and the differential thermogravimetric • 
analyses (DTGA) were carried out in order to understand the 
thermal stability and possible phase transitions taking place in 
the individual starting materials used in the preparation of the 
catalysts* These patterns enable 1) the location of the phase 
changes at different temperatures 2) determination of the 
transformation energies involved and 3) the modifications in 
the sample weights during the transformations. 

The DTA| TGA and DTGA of the samples were recorded 
simultaneously using MCM Derivatograph (Hungary) with platinum- 
platinum rhodium thermocouple as the temperature sensing element. 
The samples were heated in an alumina crucible using -AlgO^ as 
the reference. The furnace heating rate was maintained constant 
at 5°/minute. 

Differential thermal analysis has also been carried out 
on DDPOHT 900 DTA unit fitted with 1200°C furnace assembly 
using Pt-Pt Rh thermocouples with heating rate maintained 
at 12°C/minute. This has been conducted to obtain the DTA 
pattern on cooling* recording the curves corresponding to the 
reactions in the transformation region. This operation was not 
possible with MCM derivatograph and enables a clear under- 
standing on the nature of the reversibility of transformations. 
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Results of Thermal AnaXyses : 

I . Ortho-Aluminium Phosphate (o-AlPO^ ). 

The DTA, TG-A and DTGA curves are shown in Pig. 9. The 
DTA curve exhibited an endotherm commencing at 60°C attaining a 
maximum around 120°C. This endotherm was attributed to the 
dehydration of aluminium phosphate. Between 120°G and 810°C, 
phosphate does not show any transitions [129] . The TGA and 
DTGA curves exhibited a weight-loss of 10 per cent in the 
dehydration range. 

II . Bismuth Molybdates . 

The DTAj TGA and DTGA curves for Bi20^.Mo0^ are 
indicated in Pig* 10a. The DTA curve shows well defined 
endotherms with peaks at 650, 705 and 720°C in addition to the 
initial endothermic slope. These endotherms correspond to the 
structural transformation, y Y* where y' is the high 

temperature modification of y-t'ismuth molybdate. It appears 
from the DTA pattern that the y -> Y ' transformation takes 
place in several steps as could be seen from these individual 
endotherms occurring within the temperature range of the 
transformation. 

Heating and cooling curves indicated in Pig. 10b 
were obtained to understand the nature of reversibility of 
these transformations. It is seen from the figure that upon 
first heating, a large peak is observed at 686 °C with a 
shoulder at 680°G undergoes a reduction in its area on first 
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Fig. 9 :TGA,DTG and DTA curves tor 0 -AIPO 4 . Sample 
weight 356 mg, heating rate 5*C/min. 

Reference: oc-Al 203 - 
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Tempera turc,®C 


Fig. 10 a: TGA.DTG and OTA curves for 81203. M0O3. 

Sample weight 1056mg, heating rate SXimm. 
Reference: o6-Ai203. 




1st heating 


1st cooling 


2nd heating 


KJ 2nd cooling 


3 rd heating 

686 


■ig, 10 b : Differential thermal pattern tor 81203. M0O3 
indicoting the forward and backward 
transformations (arrows indicate heating, 
and cooling cycles). 
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cooling, thus, indicating that some of the transitions observed 
in the first heating curve are irreversible. Repeated heating 
and cooling of uhe same sample (second heating - second cooling - 
third heating) showed that the areas under the peaks remained 
unchanged. Hence, it can be inferred that the peak at 610°C 
and the shoulder at 680°C are associated with the irreversible 
y Y' transition, while the peak observed at 686 °C is 
due to the melting of Y' phase. 

The TGA and DTGA curves (Pig. 10a) indicated minimal 
weight-loss (aroimd 0.6%) for the sample when heated upto 800°C . 

The DTA curve for the p -bismuth molybdate (BigO^ . 2MoO^ ) 
(Pig. 11a) exhibits two endotherms. The shallow endotherm 
observed at 550°C is relegated to the solid-state reaction 
(Equation (4.1)) between a- and Y-phases present in the 
Bi20^.2Mo0^ sample. An exothermic hump developing at about 
650°C immediately followed by an intense endothermic peak 
attaining a maximum around 700°C were due to the decomposition 
of p-phase into a- and y' - phases (650°G) and subsequent melting 
of these two phases (700°G). 

The DTA pattern indicated in Pig. 11b clearly brings 
out the reversibility of a P + Y* transformations. In the 
forward run, the intense endothermic peak attained a maximum at 
700®C while the corresponding cooling curve was split into two 
peaks with maxima at 656°C and 630°C . The 656°C peak represents 
solidification of a- and Y'“Phases, the peak at 630°G indicates 
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Fig. 11a ; TGA, DIG and DTA curves for Bi 203 . 2 Mo 03 . 
Sample weight 1130 mg , heating rate 6*C/min. 
Reference. ot-Al 203 . 
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Fig.llb: Differential thermal pattern tor Bi 203 . 2 Mo 03 

indicating the forward and backward 
transformations {arrows indicate heating 
and cooling cycle). 
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the reversible reaction, a + y’ P • Thus it can be 
concluded that the p-phase is not stable beyond 630°C, One 
conclusion is apparent at this stage, i.e., the structural 
transitions observed in the present case are reversible, while 
the same were irreversible in the case of y-bismuth molybdate. 

The T&A and BTGA indicate minimal weight-loss around 
0.1% for the sample upto a tempera tui« of 800°C (Pig. 11a). 

This could be due to the sublimation of molybdic oxide occluded 
in Bi 20 ^. 2 MoO^ during its preparation. 

111 . Bi 20^ .2Mo0^-AlP0^ Catalyst . 

In order to understand the nature of the phase 
transitions undergone by the individual starting materials 
during calcination, the thermograms of uncalc ined and calcined 
Bi20^.2Mo0^-AlP0^ catalysts were obtained (Pigs. 12 and 13). 

The DTA curves for the uncalcined as also the calcined 
catalysts are identical. The endothermic peak first observed 
at 110°C IS due to dehydration. The identical nature of the 
two thermograms with no noticeable peaks beyond the dehydration 
range confirms the absence of any chemical reaction between the 
two starting materials, namely, Bi20^.2Mo0^ and AlPO^. It may 
be mentioned here that the endothermic peaks as observed in the 
pure bismuth molybdate sample (Pig. 11a) axe not evidenced in 
this case since the amount of bismuth molybdate present 
together with AlPO^ is only 4.76 per cent. As such the energies 
involved in the transformations would be correspondingly low 
with the peaks smaller than that could be noticed. 
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0*05 mol Bi203'2Mo03 per mol o~AtP04. 
Sample weight 422 mg, heating rate 5*C/min. 
Reterence : 0C.-AI2O3. 



-^Endotbei 

Weight loss Exotherr 



Temperature,®C 


Fig. 13 : TGA,DTGand DTA curves for calcined 81203. 

2M0O3-AIPO4 cotalyst. Catalyst composition 
0 05 mol Bi203.2Mo03 per mol 0-AIPO4. 
Sample weight 424 mg, heating rate 5®C/min. 
Reference : oC“Al203» 
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The TG-A and DTGA curves indicated a weight -loss of 
about 12% in the case of uncalcined catalyst sample. This loss 
corresponds to the dehydration of the aluminium phosphate 
present in the sample. The patterns for the TGA and DTGA curves 
in the Case of calcined samples, although similar in trend to 
the corresponding uncalcined samples, exhibited relatively lower 
we ight -loss (8% ) . 

4.8. Discussions and Conclusions . 

Aluminium phosphate, bismuth molybdates and the bismuth 
molybdate -aluminium phosphate catalysts have been characterised 
in detail using X-ray diffraction, infrared spectra and thermal 
analysis techniques. The bulk density and the surface area 
measurements have been carried out on the starting materials as 
also the catalysts. 

The X-ray diffraction studies for ortho-aluminium 
phosphate indicated its amorphous nature , making it difficult to 
establish its structure. However, the infrared spectra revealed 
that the phosphorous in the same was tetrahedrally coordinated 
to the surrounding oxygen atoms. Considering the fact that 
aluminium phosphates are iso-structurals with silica [127 ] it 
can be stated that the prepared aluminium phosphate had a 
network structure similar to that of Si(SiO^) network with both 
aluminium and phosphorous surrounded by tetrahedra of oxygen 


atoms. 
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X-ray diffraction studies for Y-bismuth molybdate 
(Bt20^iMo0^) established that its structure was analogous to 
the structure of mineral koechlinite* The infrared spectra of 
the Same indicated the presence of c or ner -shared molybdenum 
octahedra with terminal oxygen bonded to only one molybdenum 
atom (inferred from the band at 850 cm~ ). Such oxygen sites 
are indicative of an active dehydrogenation catalyst [l8, 90a]. 
X-Ray diffraction and infrared studies of BigO^.MoC^-AlPO^ 
catalysts indicated that the calcination of Bi20^.Mo0^ together 
with AlPO^ did not bring about any changes in its structure, 
However, the thermal analyses revealed that the Bi20^.MoO^ 
undergoes a phase transition, Y ^ at temperatures of 630°C 

and above. This phase transition is irreversible and results 
m a loss of catalytic activity [18, 90a]. 

X-Ray diffraction studies for the p-bismuth molybdate 
(Bi20^ . 2MoO^ ) indicated the presence, in restricted quantities, 
of a- and p -bismuth molybdates. The infrared spectra also 
confirmed the presence of a-phase in the prepared Bi20^.2Mo0^ 
from the occurrence of sharp bands in the range of 900-950 cm 
(characteristic of a-bismuth molybdate). The DTA analysis 
indicated that these a- and y - bismuth molybdates recombined at 
a temperature of 550°G. The DTA analysis also indicated that the 
thermal decomposition experienced by p -bismuth molybdate above 
630°G was of reversible nature. 
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The X-ray diffraction studies evidenced that the 
calcination of Big 0^. 2J/Io0j with o-AlPO^ resulted in the inhibition 
of the decomposition of y^ismuth molybdate into c£~ and Y-phases. 
Similar observations were reported by Batist et#al. [16] in the 
case of BigO^ . 2MoO^/silica gel Catalysts. It may be stated that 
the ortho -aluminium phosphates are iso-structurals with one or 
more forms of silica having a direct parallelism between all 
crystalline modifications of AlPO^ and SiOg [127]. The DTA 
analyses of the uncalcined and the calcined Bi20^,2Mo0^-AlP0^ 
indicated no reaction between the two constituents. This was 
also confirmed in the infrared analyses. 

In the case of a-bismuth molybdate (Bi20^,3Mo0^) the 
X-ray diffraction data indicated that the compound has a 
monoclinic structure. In addition, the X-ray results have also 
revealed that the calcination of Bi20^.3Mo0j together with AlPO^ 
did not bring about any changes in its structure. Since the 
monoclinic structure of Bi20^.3Mo0^ results in edge-sharing 
molybdenum octahedra, this compound was not considered as a good 
dehydrogenation catalyst [18]. 

The absence of any thermal transition in the normal 
temperature range (400-500 °C) experienced in the oxidative 
dehydrogenation of butane and butenes suggests that both BigO^.MoO^^- 
AlPO^ are potential catalysts for the oxidative dehydrogenation 


of butane. 
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CHAPTER 5 

STUDIES ON BISMUTH MOLYBDATE-ADUMINIUM PHOSPHATE CATAXYSTS 
USING RESPONSE SUEPACE METHODOLOGY 

5.1. Introduction . 

Meagre information seems to be available in literature 
on the oxidative dehydrogenation of butane. In the present 
investigation, in which the suitability of Bismuth 

Molybdate-Aluminium Phosphate as a catalyst has been studied, 
the chemical process has been treated as a black box and its 
input-output relationships were studied using response surface 
methodology. The inputs to the black box were the experimental 
variables like the temperature of reaction, contact time, 
butane-oxygen ratio, oxygen-nitrogen ratio, composition of 
catalyst, method of catalyst preparation etc. The output of 
the black box could be the selectivity or the yield of any 
particular product and the like. Choice of the response and 
the independent variables has been discussed in Section 5.5. 

An introduction to the response surface methodology is given 
below. 


5.2, Response Surface Methodology . 

Let y be the response of a chemical process dependent 
on the levels of k quantitative factors or process variables 




f ^2 * * * 


, 3^ which can be precisely measiired and controlled. 



The model for the u combination of factor levels is given by, 

~ ***r U * 1, 2, *t‘t H (5*l) 

where N is the number of experiments and 0 is the functional 

relationship. It is further assumed that the response y^ has a 

_ 2 

mean y^ and variance a • 

A geometrical portrait of the response function in 
the factor space is called a response surface. 

The experimental region R, is a bounded subspace of 
the whole factor space, R is bounded because of practical 
limitations. The response surface methodology deals with the 
problem of locating a point (x°, x^, x^) within R at which 

y IS an extremum [3, 29, 31, 32, 38, 40, 59 ]» The response 
surface methodology, originally introduced by Box and Wilson 
[31 ], is essentially applied in two stages in order to find the 
optimum conditions. The two stages are termed the first-order strat- 
egy and the seccnd order strategy respectively. The methodology 
is applied as follows: 

After conducting preliminary experiments using any 
suitable experimental design, a hyper-planar surface is fitted 
to the responds data, and the adequacy of the model is tested. 

The direction of steepest .ascent/descent is established and the 
factor level combinations lying along this direction are 
computed. Experiments will have to be performed at these 
factor levels to verify the predicted value of the response. 

One of these factor level combinations is chosen as a base level 



for the next experimental design and then the experiments are 
perfoimed around these chosen base levels. The whole procedure 
IS repeated till a near-stationary region is reached where the 
hyper-planar approximation of the response surface is no more 
valid (i.e. the first order strategy fails). Thereupon, a 
second order strategy is employed in which the interactions and 
the higher order effects of the various factors (i.e. non- 
linearity) are determined by augmenting the experimental design. 

The traditional one factor-at-a-time method, where 
the experiments are earned out holding all the factors 
constant except the one mder consideration, often fails to 
locate the true optimum [31]. 

The response surface methodology cuts down experimental 
efforts by making use of experimental designs which permit the 
experimenter to assess the strength of the interactions between 
the various factors while varying them simultaneously [38], 

In the following section, the results of the feasibility 
studies of the oxidative dehydrogenation of butane are presented. 

5 . 3 . Molybdic Oxide -Aluminium Phosphate Catalyst — A Peasibilitv 
Study . 

Exploratory woik: of Gaspar and Pasternak [50 ] 
indicated that molybdic oxide -aluminium phosphate catalysts 
could dehydrogenate n-butane to 1, 3-butadiene in the presence 
of oxygen. Since the aim of the present investigation was 
to dehydrogenate n-butane to 1, 3-butadiene in one step, 
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experiments were first conducted using the molybdic oxide- 
aluminium phosphate (MoO^-AlPO^) catalyst prepared in the 
laboratory. This gave an estimate of the molybdic oxide 
content of the bismuth molybdate -aluminium phosphate Oatalysts 
subsequently tested. 

MoO^-AlPO^ catalysts prepared according to the 
procedure of Section 3«4 were used in the preliminary experiments. 
The catalyst con^osition was varied from 0 to 0,3 mol MoO^/ 
mol AlPO^, The experiments were conducted at the following 
levels of the variables: 


Temperature 
Catalyst -feed ratio 
Butane-oxygen ratio 
Diluent 


550, 600 °C 

0*75, 1 • 25 g/ (mgmol/minut e ) 
1:1 
Bil 


Experimental data are presented in Appendix VII, for 
runs numbered 1 through 21 , The results are shown graphically 
in Pigs. 14 through 21. 

The overall conversion of n-butane over varying 
catalyst composition is shown in Pig, 14. Catalyst compositions 
of 0.15 mol MoO^/mol AlPO^ and above seem to have very little 
influence on the overall conversion of butane. The effect of 
the molybdic oxide-aluminium phosphate ratio on the yields of 
1, 3-butadiene, 1 -butene, propylene and ethylene are shown in 
Pigs. 15 through 18. These figures show that the individual 
yields of olefins and diolefins pass through a maximum at a 


f 



Yield of butadiene 



•c 

gmin/gmol 


550 

1-25 

- X 

550 

0.75 

0 

600 

1-25 

• 

600 

0-75 



Mol M0O3/M0I AIPO4 

Fig. 14 : Molybdic oxide -aluminium phosphate catalyst 

overall conversion. 







Mol M0O3/ Mol AIPO4 


Fig. 17 ; Molybdic oxide -aluminium phosphate catalyst' 

Propylene yield. 





\o\Q ot ethylene /100 Mole butane 


®C g min /g mol 

550 V2S 
K 550 0-75 

o 600 1*25 

’• 600 0-75 





0 0-1 0-2 0-3 

Mol M0O3/ Mol AlPO/i 

Fig. 18 ; Molybdic oxide-aluminium phosphate catalyst: 


Ethylene yield. 







Fig. 20 : Molybdk oxide-aluminium phosphate catalyst: 


Selectivity to olefins and dioletins. 



Fig. 2T :Molybdlc oxide- aluminium phosphate catalyst: 

Formation of carbon oxides. 
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catalyst composition of about 0*1 mol MoO^/mol AlPO^, Pigs. 19 
and 20 show that the peak selectivity to total olefins 
diolefins together with their peak yield ocdUrs afound the 
same catalyst composition. 

These results show that the best catalytic activity 
in the case of MoO^-AlPO^ catalyst occurred at a catalyst 
composition of 0.1 mol MoO^/mol AlPO^. Howeveri the activity of 
the catalyst for the formation of undesired products like CO 
and COg (shown in Pig. 1 1 ) is also the highest at the same 
catalyst composition. 

The following sections are devoted to the enhancement 
of (butene + butadiene )yield by experimenting on the constituent 
elements of the catalyst and the process variables affecting 
the oxidative dehydrogenation. 

5.4. Studies on Bismuth Molybdate -Aluminium Phosphate Catalysts . 

Having established the feasibility of the one-step 
oxidative dehydrogenation of n-butane to 1, 3 -butadiene, 
attempts were made to in^irove the gross yield of butenes and 
butadiene by incorporating bismuth in the catalyst. Bismuth 
oxide when present together with molybdic oxide enhances the 
selectivity of butene to butadiene reaction. Hence the 
oxidative dehydrogenation of butane was further investigated 
over bismuth molybdate -aluminium phosphate catalysts. Bismuth 
molybdate catalysts were prepared according to the procedures 
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outlined in Section 3.5. The different factors affecting the 
oxidative dehydrogenation process are analysed in the following 
section. 

5.4.1. Identification of Factors and the Experimental Region . 

Preliminary experiments (not reported herein) showed 
that for the reactor diameter used and the operating conditions 
chosen, catalyst dilution was not necessary. A catalyst 
particle size range of 1-2 mm was used throughout the 
inve st igat ion. 

Preliminary experiments also indicated that air 
enriched in oxygen was needed to get enhanced yield of butadiene. 
Further, in order to reduce the experimental effort, the 
oxygen -nitrogen ratio was kept constant at 0.5. 

Of the numerous factors affecting the process, the 
following four were considered to be the most important: 

1. Reaction temperature. 

2. Catalyst-feed ratio. 

3. Butane -oxygen ratio. 

4. Composition of the catalyst. 

The range of experimentation is chosen based on the 
following considerations: 

(^) Reaction Temperature : The lower limit of the reaction 
temperature was fixed at 400°C , since at temperatures of 350°C 

below, oxygenated products like alcohols are reported to 
be formed in preference to dehydrogenated products [53] - 
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The upper limit was set at 500°C since the hismuth molybdate 
catalysts tend to loose their activity above 500°C [53]. 

( 2 ) Catalyst -Feed Ratio ; The catalyst -feed ratio (W/p) was 
chosen in the range of 0 , 3 - 2. 6 g catalyst/{mgmol/minute ), This 
was equivalent to a contact time of about 0.25-2 seconds, a 
range normally used in practice. 

(3) Butane Oxygen Ratio ; The choice of butane-oxygen ratio 
was governed by the explosion range of the butane-air mixture. 
The explosition range of butane is 1.9-8. 5 vol % in air [ 86 ], 
Therefore, any butane-air mixture having a butane to oxygen 
ratio between 0.092 and O .445 is explosive and hence cannot be 
used in the experiments. 

At butane to oxygen ratios of 0.09 and less, oxygenated 
compounds tend to dominate the product stream [53, 61], The 
upper limit of butane to oxygen ratio was fixed at 2.0 since 
at very low oxygen concentrations, bismuth molybdate tends to 
deactivate rapidly due to the loss of lattice oxygen. 

( 4 ) Composition of the Catalyst ; The bismuth molybdate- 
aluminium phosphate catalysts were prepared with different 
proportions of bismuth molybdate to aluminium phosphate. The 
composition of bismuth molybdate was also treated as a variable, 
since bismuth molybdates of different phase compositions show 
good activity for the oxidation dehydrogenation of butanes 

[53, 61, 104, 115], 
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Two phases of bismuth molybdates, namely^ P- and Y- 
with compositions Bi 20 ^. 2 Mo 0 ^ and Bi^O^.MoO^ respectively were 
used in the present experiments. Structural studies of the 
bismuth molybdates discussed in Chapter 4 concluded that ct- 
bismuth molybdate (BigO^ ,3MoOj ) with its edge shared MoOg 
octahedra structure was likely to be less active compared to 
the p- and yP^^s^ses. Moreover, the Catalyst grannules of 
Bi20^ .3Mo0^-AlP0^ catalyst were very friable. Hence this 
catalyst was not used in the present study. 

The p- and .Y -phases of bismuth molybdates were 
treated as qualitative variables in the present response 
surface strategy. The nature of the qualitative variable 
demands that the strategy has to be applied in parallel for 
these two phases. Therefore, experiments were conducted 
side-by-side on these two phases in the linear region of the 
response surface in order to locate the best possible values 
of the optimisation factor (to be defined in the following 
section) obtainable on them. Once the approximate extrema were 
located, a decision was taken, as to which of the qualitative 
factors should be used in the second order strategy. 

The ranges of experimentation for the various 
factors used in the present study are summarised in Table 14. 
The choice of the optimisation parameter is discussed in the 
following section. 
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Table 14* EXPEEIMEUTAL REGION OF FACTORS USED IN THE STUDY OF 
OXIDATIVE DEHYDROGENATION OF BUTANE 


S.No. 

”7 

! 

} Factor 

r 

1 

! [ 

ICodei 
INamel 

I 1 

t f 

Experimental Region 

1. 

Temperature 

^1 

400-5 00 oc 

2. 

Catalyst-feed ratio, 

W/F 

^2 

0. 5-2.6 g-minut e/mol 

5. 

Butane -oxygen ratio 


0. 5-2.0 

4. 

Catalyst composition 
(mol bismuth molybdate 
per 100 mol aliuninium 
phosphate ) 

^4 

0.0 and above 

5. 

Phase of bismuth 
molybdate 


p -bismuth molybdate 
(BigO^ . 2MoO^ } 

T -bismuth molybdate 
(BigO^ .MoO^ ) 
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Oxygen consumption was an important variable as it 
would determine the total volume of gases to be handled in the 
up-stream and down-stream units of an actual chemical plants 
In the present study, it was observed that the oxygen consumption 
was always almost complete and hence it could not be considered 
as an optimisation parameter. 

Based on the above considerations, the total yield of 
(butenes + butadiene) was chosen as the optimisation parameter 
(response) in the present investigation. 

5 . 4 . 3 . Selection of an Bxnenmental Design . 
k-1 

A 2 fractional factorial design was used in the 
present study. According to this design, 8 experiments had to 
be performed for each qualitative variable. The fractional 
factorial design reduces tho amount of experimentation without 
loosing at the same time the advantages of symmetry, ortho- 
gonality, normalisation and rotatability required of a good 
experimental design [38]. The details regarding the fractional 
factorial design are given in the following sections. 

5 , 5 . First Order Response Surface Strategy: Annlication of 
Steepest Ascent to Move to a Near-Stationary Regio n. 

Having defined the factors governing the process and 
the process response to be optimised, a base level had to be 
chosen within the experimental region of the factor space, R. 

This experimental region has already been given in Table 14. 



The base level should be the combination of factor levels 
corresponding to the best conditions determined from an 
analysis of the a-prion information. 

itirther^ the variation interval for each factor had to 
be chosen. The variation interval of a factor when added or 
subtracted from the base level gives the upper or the lower 
level of the factor respectively. 

To simplify the recording conditions of an experiment, 
processing of the experimental data the values of the 
were coded as follows: 

X = -JlZ 10 3 = 1,2, ..., k (5.2) 

1 

the coded value of the 3 -th factor, 
the natural value of the 3 -th factor, 
the natural value of the base level, 
the variation interval of the 3 -th factor. 

Conventionally, the coded value of the upper level 
corresponds to + 1 , the lower level to -1 and the base level to 0 . 

The first order approximation of the response surface 
resulted in the simplification of the functional form of 0 
in Equation (5.1) to the linear fonn, 

y = bo b,x, + + bjXj + b^x^ (5.3) 

where b^ , bg, b^ and b^ are the coefficients associated with 
the factors x^ , X 2 , and x^ respectively. The above 


and the 
factors 


where 


30 


IS 

IS 

IS 

IS 
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coefficxents together with were evaluated experimentally. 

The experiments were performed according to a fractional 

factorial design [38] shown in Table 15. 

In this design, the experimental levels of factor x^ 
were generated using the relation, 

^4 x^x^x^ (5.4) 

Such a design gave the product x^X 2 X^x^ = +1. This product is 
called the defining contrast. Which effect is confounded with 
whioh interaction, was calculated from the defining contrast [38] . 
For the above design the confounding pattern may be seen to be, 


1 

= 234 

2 

= 134 

1 

= 124 

i 

= m 

12 

= 14 

11 

= 

21 

= li 


This means that the main effect of variable x^ was confoimded 
With the interaction of the variables 2, 3 and 4 (i.e. XgX^x^) 
and so on. The estimates b^ , b 2 > and b^ will be the 
estimates of their true values provided the third order 
interactions of the effects are negligible. 

The coefficients b^, b^ , b 2 , and b^ can be 
estimated in the least squares sense as, 



Table 15. DESIGI MATRIX FOR 2 ^ ^ FRACTIONAL FACTORIAL 
EXPERIMENT 


r 

Design ! 
Trial i 
No. 1 

f 

* 

^0 

^1 

^2 

i 

X3 I 

T 

1 

! 

f 

f 

? 

^4 1 

(= X^XgX^) 1 
! 

Response 

y 

1 

4-1 

-1 

-1 

-1 

-1 

yi 

2 

+1 

+1 

-1 

-1 

+1 

yp 

3 

+1 

-1 

+ 1 

-1 

+1 

^3 

4 

+1 

+1 

+1 

-1 

-1 

^4 

5 

+1 

-1 

-1 

+1 

+1 

^5 

6 

+1 

+1 

-1 

+1 

-1 

^6 

7 

+1 

-1 

+ 1 

+1 

-1 

^7 

8 

+1 

+1 

+1 

+1 

+1 

^8 


* Xq IS a dummy variable bearing a constant value of +1 . 
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^ , T ^-1 T 

b = (x x) Z 


(5.6) 


where x is the design matrix, £ is the vector of observed 
A 

yields and b is the vector of coefficients. In the case of 
present experimental design, can be alternatively calculated 
from the equation, 




b 


3 


N 

■C” 

/- 

1=1 



1=1 


^I3=^x 



3 — ^ > 2, ...,1c 


(5.7) 


The first-order strategy was applied till such time the 
Equation (5.3) described the yield-response surface within 
experimental errors. This was verified by testing the adequacy 
of the model and the significance of coefficients describing 
the model. 


Adequacy of the Model ; The experimental variance of the 

2 

optimisation factor, s was compared with the variance of 

*/ 

adequacy of the model, (defined below) by constructing an 

E test [60], The test is described below. 

The experimental error variance of the optimisation 
parameter was calculated using the formula, 



n , 

^ (Zq - 

9 ^ 1 

n - 1 


(5.8) 


where, 
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o 

s is the variance of the optimisation parameter j* n is the 
number of repeat observations ^ y is the mean value of the 

2 

repeat obseivation. The variance of adequacy of the model s^^ 
was calculated as, 


r: (1 - y, )‘ 
.2 _ 1=1 ^ ^ 
’ad ” N - (k + l) 


(5.9) 


where, 

s^, IS the variance of adequacy of the model 
ad 

y^ IS the predicted value of the optimisation parameter. 

Equation (5»3) was used to calculate the predicted values of 
the total yield of (butene + butadiene). The adequacy of the 
model represented by Equation (5.3) was tested by constructing 
the Eisher ratio (E test) as, 


E 


calc 



( 5 . 10 ) 


The model was deemed inadequate if , 

Wo> I'o.os'” 


where , 

(N - (k+1 ) , n - 1 ) is the tabulated value of E at a 

0.05 

significance level of 0.05 with E - (k+1) degrees of freedom 
(d.f.) in the numerator and n-1 degrees of freedom in the 
denominator. The tabulated values of E are found in most 
books on statistics [60] . 
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Significance of the Coefficionts : Significance of each 

coefficient was tested "by constructing the confidence intervals 

2 

[3, 129]. First the variance of the regression coefficient 

5 

b was calculated as, 

J 


fL 

N 


(5.11) 


The confidence interval, was then calculated ast 

2 


= ± t 


0.05 b. 


( 5 . 12 ) 


where , 

IS the standard deviation of the regression coefficient. 

tn^/MT IS the tabulated value of the "student's t" at 

u. ub 

5 % significance level and n-1 degrees of freedom [60], 

V/hen some coefficients of a model were not 
significant, the variation intervals of the factors associated 
with these coefficients are increased in order to test whether 
these coefficients are really not significant. If they are 
still found to be not significant, then the terms aasociated 
with these coefficients are deleted from the model and again 
the adequacy of the model is tested. 

Once the model fitted to Equation (5.5) was found to 
be adequate the steepest ascent technique was applied to 
locate the ccmbination of factor levels which improve the 
yield. 
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Steepest Ascent Along the Response Surface : The direction of 
steepest ascent was determined using the relation, 




c>x 


1 


u. + "Uo + 

1 2 


4 - 


■k 




( 5 . 13 ) 


where, 

0 IS the function describing the response surface, 
V0 IS the gradient of the response function. 


S»0 


dx 


IS the partial derivative of the function with respect 


3 "t" V\ 

to the 2 factor. 


u,| , U2, Uj^ are the unit vectors in the direction of the 

coordinate axes. 

It can be easily verified that the components of the 

gradient ^^0 i.e., ,3 = 1 , 2 , ...,k are the same as the 

D 

regression coefficients in Equation (5.3). The details 
regarding the numerical evaluation of the direction of the steepest 
ascent are found in Sections 5.5.1 and 5.5.2. 

The steepest ascent calculations were used to locate 
the possible combination of factor levels along the steepest 
ascent direction which predicted increasing yields of 
(butenes + butadiene). The predicted yields of (butenes + 
butadiene) were not experimentally verified, since the 
experimental range in the present investigation was narrow. 
Moreover, the step sizes chosen were small, so that the 
predicted responses were bound to be observed experimentally. 



One of the above combination of factor levels was 
used as the base level for the next experimental design. This 
above procedure was repeated. This iteration process was 
discontinued when the response surface was no more describable 
by Equation (5.3). This indicated that the factor levels were 
now close to the stationary region of the response surface. 

The strategy was applied in parallel to the two 
qualitative factors namely, p- and Y-phases of bismuth molybdat 
as described in Section 5.4.1. Once the best possible yields of 
(butenes + butadiene) attainable with these two factors were 
known, one of these factors was eliminated based on the values 
of the best yields obtained and the process conditions at 
which these best yields were realised. The second order 
response surface strategy was then applied only to this 
qualitative factor. 

A base level of quantitative factors will have to 
be chosen to begin the search for the optimum yield. A 
discussion on the choice of initial base levels is presented 
in the following section. The first order strategy as applied 
to the two qualitative factors is described in Sections 5.5.2 
and 5 . 5 . 3 . The two qualitative factors are discriminated 
against each other in Section 5*6 where the second order 
response surface strategy adopted has also been discussed. 
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5.5.1. Choice of a Base level . 

Having chosen the factors and the experimental region 
of the factor space , the optimisation parameter and the 
experimental design (Section 5.4), experiments were conducted 
to search for an optimum yield of (butenes + butadiene). 

The search was initiated by choosing a base level of 

4-1 

factors, and setting up a 2 fractional factorial design 
around this base level. This base level of factors was chosen 
on the basis of a-priori knowledge. The levels of temperature, 
the butane to oxygen ratio and the catalyst to feed ratio, 
were chosen close to the best values reported in literature 
for the butene oxidation dehydrogenation process [55, 61, 

104, 115] . 

The proportion of bismuth molybdate to aluminium 
phosphate was chosen as follows: Experiments on MoO^-AlPO^ 
catalysts had showed an optimum ratio of molybdic oxide to alumincum 
phosphate to be in the range of 0.1-0.15 mol per mol. The - 
bismuth molybdate -aluminium phosphate catalyst should contain 
about 0.1-0.15 mol of bismuth molybdate per mol of aluminium 
phosphate. The composition of the bismuth molybdate was, 
therefore, fixed at 20.0 mol bismuth molybdate/lOO mol 
aluminium phosphate. 

The variation interval of temperature was chosen to 
be 25 °C, large enough to take care of the non-isothermal 
effects in the catalyst bed. The variation intervals of other 
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factors were chosen rather arbitrarily. 

Table 16 lists the base levels and the variation 
intervals of the various factors used for the first move on 
the response surface. A fractional factorial design was set 
up using these values and the response surface methodology was 
initiated using the experimental data obtained. Experiments 
conducted on 312©^ ,2MoO^-AlPO^ and BigO^ .MoO^-AlPO^ are found 
in Sections 5.5*2 and 5.5.3 respectively. 

5.5.2. Optimisation of Butenes + Butadiene Yield Using 
Bi^O ? . 2MoO^ -AlPO^ C at aly st . 

In the present section, experiments performed on 
BigO^ . 2Mo0j-AlP0^ catalysts are dealt with. The response 
surface methodology was initiated by choosing the base lovel 
of the factors and their variation intervals at the values shown 
in Table 16. 

Fractional factorial design used in the experiments 
IS shown in Table 17. The experimental trials were run in a 
random manner to avoid any possible bias on the part of the 
experimenter. The randomisation sequence is shown in column 3 
of the table. The experimental results are given in Appendix 
YII for runs numbered 22 through 31. The yield of (butenes + 
butadiene) was calculated as the moles of total butenes and 
butadiene formed per 100 moles of butane fed. 

Equation (5.3) was fitted to the experimental data, 
and the coefficients b^, b^, bg, b^ and b^ were estimated 
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Table 16. FIRST ORDER RESPONSE SURFACE STRATEGY. 
INITIAD BASE LEVEIS AND VARIATION INTERVADS OF FACTORS. 


Factor 

^ 

|C ode 

t 

1 

IBese 

1 

i 

{Variation 

— 

1 

IPactor Level 


1 

1 

« 

! 

|Ievel {Interval 

f t 

t ! - 

-1 

1 

+1 

Temperature °C 

Xi 

450 

25 

425 

475 

W/F g-min/mg mol 

Xp 

1 .3 

0.3 

0.7 

1 .6 

Butane/oxygen ratio 

X3 

1 .0 

0.25 

0.75 

1.25 

Moles of bismuth 


20 

5 

15 

25 

molybdate/l 00 mol 
aluminium phosphate 

4 






Type of bismuth molybdate Qualitative Factor; P- and Y -phases. 
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Table 17.' FIRST ORDER RESPONSE SURFACE STRATEGY. 

FIRST MOVE ON THE RESPONSE SURFACE 

BigO^-SMoO^-AlPO^ CATALYST. DESIGN MATRIX, EXPERIMENTAL 
RESULTS AND TEST FOR MODEL ADEQUACY. 


t f t 

Run {Design {Rand- { Factor Levels 

No. {Trial {omised { * ** 

{No. {Sequence { x^ x.^ x^ 

t f 

^4 i 

_ t 

Yield 

Expt . 

y 

Pred. 

y 

22 

1 

10 +1 

-1 

-1 

-1 

-1 

5.751 

5.9C3 

23 

2 

5 +1 

+1 

-1 

-1 

+1 

0.503 

1 .0^2 

24 

3 

1 +1 

-1 

+1 

-1 

+1 

2.029 

2.719 

25 

4 

7 +1 

+1 

+1 

-1 ■ 

-1 

1.152 

0.656 

26 

5 

4 +1 

-1 

-1 

H-1 

+1 

4.346 

6.712 

27 

6 

8 +1 

+1 

-1 

+1 

-1 

1.746 

2.649 

28 

7 

9 +1 

-1 

+1 

+1 

-1 

8.649 

6.326 

29 

8 

2 +1 

+1 

+1 

+1 

+1 

4.171 

1.463 

Reoeat Trials 








30 

9 

3 +1 

-1 

+1 

-1 

+1 

2.707 


31 

10 

6 +1 

-1 

+1 

-1 

+1 

3.298 


_2 

N 

IZ (Yi 

i=1 ^ 


20.4461 

= 6.815, 

d.f. = 3 


^ad “ N - (k+1 ) 

8-(4+T7 



= 0.4065, 

d.f. = 2 







^calc " 0*.4065 ■" 

^0.05^5’ 2) = 19.164 

^0 

db 

.05 

0 ^ 

= 2.92, 

- ■*^0.05 

at 

v = 

d.f. = 2 

= 0.6565 



^calc^ ^0.05^^’ 


Model is adequate. 


* Xq IS a dunaiiy variable having a constant value of +1 . 

** mol (butenes +■ butadiene) formed per 100 mol. of butane fed. 
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using Equation (5.7). The yields of the repeat trials 3, 9 
and 10 were averaged before making the calculations. The yield 
of (butenes + butadiene) was represented by the following 
—Equation- 

y = 3.6245 - 1.7315 - 0.893 x^ + 1.035 x^ - 0.700' 

(5.14) 

The above model was tested for the adequacy of fit and the 

significance of the regression coefficients. These calculations 

described below are summarised in Table 17. 

The yield values reported in trials 3, 9 and 10 were 

used to calculate the error variance of the yield of (butenes + 

2 

butadiene), s , employing Equation (5.8). The error variance 
«/ 

was found to be 0.4065 at a degrees of freedom of 2. 

The adequacy of the model represented by Equation 

( 5 . 14 ) was tested next. The predicted values of the yields 

2 

were calculated using Equation (5.14) and s^^, the variance of 
adequacy of the said model, estimated using Equation (5. 9), was 
6. 81 5 at a degrees of freedom of 3. 

The P ratio, calculated using Equation (5.10), was 
found to be 16.56. The tabulated value of P [60] at a 
significance level of 0.05 and degrees of freedom (3, 2) is 
19.164. This value is denoted as Pq 05 ^^’ Table 17. 

Since less than P^ q^( 3, 2), Equation (5.14) represents 

the response surface of yield adequately at 5% significance 
level. 



116 


Having established the adequacy of the model, the 
significance of the estimated coefficients was tested. The ''t" 
value at a significance level of 0.05 and degrees of freedom 
of 2 was found to be equal to 2^92 [60]. The confidence 
interval, /Jib. of the regression coefficients was calculated to 

vl 

be + 0,6565, using Equation (5.12), All the coefficients of 
Equation (5,14) were fo\md significant since their absolute 
values were larger than their confidence intervals. 

The adequacy of the linear model indicated that the 
subregion in which the response surface was explored, was away 
from its near-stationery region. Hence the steepest ascent 
technique described in the previous section was applied to the 
yield surface represented by Equation (5.14) to locate possible 
combination of factor levels which enhance the yield. The 
direction of the steepest ascent was located using Equation 
(5.13), in which the components of the gradient (the partial 
derivatives of response function 0 represented by Equation (5-14)) 
are the regression coefficients in Equation (5.14). Thus by 
changing the independent variables x^ , x^ and x^ in 
proportion to the values of their regression coefficients, the 
movement along the steepest path may be realised. From the 
fitted model, it will be now interesting to compute these 
factor level combinations which predict an increase in the 
yield of (butenes + butadiene). These calculations, summarised 
in Table 18, are described below. 
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Table 18* CALCULATION OP STEEPEST ASCENT* FIRST MOVE ON THE 
RESPONSE SURFACE. 

Bi oO, * aioO, -AlPO . C AT ALYST . 

2 :) D 4 



h 

I 

♦ 

t 


X3 

mol 

mol 

\ 

i "^4 

1 

! mol 
ilOO mol 

ICalculatu(? 

'Yield 

i ^ 

1 . Base level 

450 

1.3 

1.0 

20 


2. Variation interval (unit) 

25 

0.3 

0.25 

5 


5. Estimated slope (change 

in yield per unit change 
in x) 

-1 .732 

-0.893 

1.104 

-0.700 


4. Unit 

-43.3 

-0.268 

0.276 

-3.5 


5. Step when x-j changes by 

22 

-22 

-0.134 

0.14 

-1.8 


6 . 5Jhe direction of 1 

428 

1.17 

1.14 

18.2 

6.42 

steepest ascent 2 

406 

1.03 

1 .38 

16.4 

9.22 

represented by a 5 

428* 

1.03 

1.38 

16.4 

8.70 

series of 4 

428* 

0.90 

1.42 

14.6 

8.95 

possible trials 5 

on it 

428* 

0.76 

1.56 

12.8 

10.24 


* Temperature lerel was held at this level so that it may be used 
as basic level in the next step. 
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The base levels and the variation intervals of the 
factors are shown in the first two rows respectively of the 
table. In the third row, the regression coefficients of the 
response surface model (Equation (5.14-)) are given. These 
coefficients also represent the gradient of the yield in their 
respective directions. 

In order to follow the path of the steepest ascent, 
the factor levels were changed in proportion to their slopes 
(i.e. the respective regression coefficients). The fourth row 
was obtained by multiplying the variation interval of each 
factor with its regression coefficient. Physically these 
values represent the coiiq)onents of the steepest ascent vector 
along the coordinate directions of the factor space. The 
choice of the size of this vector is a matter of experimental 
convenience. In the present case, a step corresponding to 
22°0 in the direction of factor x^ was chosen and the corres- 
ponding changes in Xg, x^ and x^ are shown m the fifth rcw. 
These values were subsequently roimded off depending upon the 
accuracy with which the factor was measured in the laboratory. 

The sixth row was obtained by algebraically adding 
the calculated values, in row 5, of changes in the levels of 
the variables (equivalent to unit step in that direction) to 
the original base levels in row 1. In other words, rows 1 and 
5 were added together element by element giving the first set 
of values in rov/ 6. This process was iterated giving 
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different sets of values in row 6, which lie in the direction 

of steepest ascents 
"till 

The 6^ column in Table 18 represents the expected 
yields at the combination of factor levels indicated by the 
corresponding row. These yields were calculated using Equation 
(5. H). after coding the factors in the trials using Equation 
(5.2). 

It may bo seen that, in the second trial of sixth 
row, the level of the factor x^ approaches the lower bo\md of 
its experimental region given in Table 14. The movement in 

the direction of steepest ascent below this value was 

experimentally not feasible. Moreover, any future experimental 
investigation of the yield along this direction required that 

a new base level be chosen for the experimental design. The 

new base level chosen should be such that the upper and the 
lower levels of the factors used in the next experimental 
design should be within the region of experimentation. The 
factor x^ (temperature) was constrained at a level of 428°C 
in trial 3 , 4 and 5 , so that with a variation inteiwal of 
25 °C, the lower level of the same was within the experimental 
region. It is seen from column 6 that the yields of 
(butenes + butadiene) calculated along the steepest ascent 
direction show an increasing trend. 

The response surface strategy requires that the 
predicted yields along the steepest ascent direction be 
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experimentally verified. ThiSi in the present case would mean 
that the catalysts of varying compositions (trials 1, 2, 3, ... 
in the steepest ascent direction) be prepared to verify the 
experimental yields. The experimental effort in preparing 
and testing a single catalyst was considerable, while the same 
in conducting the run was much less. Hence, one of the factor 
level combinations in the steepest ascent direction was chosen 
and a set of experiments were conducted around these base levels. 
The experimental design for this situation is described below. 

Second Move Towards the Hear-Stationarv Region : The new base 
level was chosen on the basis of the calculations in the last 
row of Table 18. The factor level combinations which predicted 
a yield of 10,24 percent was chosen as the base level for the 
present experimental design. The base levels and the variation 
intervals of the factors so chosen are shown in Table 19. 

Here, the variation intervals of factors Xg and x^ were 
increased over the previous values (shown in Table 16) in 
order to improve the significance of the coefficients associated 
with these factors. The variation intervals of factors x^ 
and x^ were not changed. 

The experimental design matrix, the randomised 
sequence of trials and the corresponding experimental yields 
of (butenes + butadiene) are presented in Table 20. The 
experimental data are given in Appendix VII under runs numbered 
32 through 41 . A glance at the reported yields in 
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Table 19. FIRST ORDER RESPONSE SURFACE STRATEGY. 

SECOND MOVE ON THE RESPONSE SURFACE. 

Bi20^.2Mo0^-AlP0^ CATALYST. 

BASE LEVELS AND VARIATION INTERVALS OF FACTORS. 


Factor 

{ Code 

! 

f 

f 

f 

1 

» 

{ Base 

J Level 

1 

t 

f 

1 

{Variation 

{Interval 

f 

Factor Level 
-1 +1 

Temperature, 

^1 

425 

25 

400 

450 

W/F, g min/mgmol 

^2 

0.8 

0.4 

0.4 

1 .2 

Butane/oxygen ratio 


1.5 

0.25 

1.25 

1.75 

Mol of bismuth/lOO mol 
aluminium phosphate 

^4 

13 

8 

5 

21 
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Table 20 < FIRST ORDER RESPONSE SURFACE STRATEGY. SECOND 
MOVE ON THE RESPONSE SURFACE. 

Bi20^.2MoO^-AlPO^ CATALYST. DESIGN MATRIX, EXPERIMENTAL 
RESULTS AND TEST FOR MODEL ADEQUACY. 


r 1 r 

Run {Design {Rand- | 

No. {Trial {omised { 

{No. {Sequence} 

j 1 _i_ 

^o 

Factor Levels 

*4 

1 ** 

{ Yield 

{ Expt , Pred. 

t 

1 y y 

32 1 

10 

+ 1 

-1 

-1 

-1 

-1 

12.146® 

12.511 

33 2 

6 

+1 

+ 1 

-1 

-1 

+1 

5.499 

5.693 

34 3 

4 

+1 

-1 

+ 1 

-1 

+1 

6.039 

4.089 

35 4 

3 

+1 

+1 

+ 1 

-1 

-1 

11.538 

12.960 

36 5 

9 

-h1 

-1 

-1 

+1 

+1 

2.278 

3.870 

37 6 

5 

+1 

+1 

-1 

4*1 

-1 

14.810 

12.815 

38 7 

2 

+ 1 

-1 

+ 1 

+1 

-1 

10.902 

11.131 

39 8 

1 

+1 

+1 

+ 1 

+1 

+1 

4.132 

4.312 

Repeat Trials 









40 9 

8 

+ 1 

-1 

-1 

-1 

-1 

12.532 


41 10 

7 

+1 

-1 

-1 

-1 

-1 

12.105 


N 

2 1=1 

A 

1 - ^1 

)2 

13 

i.44 

4.48 


d.f. = 

3 

®ad “ N - 

(k+1 ) 



3 

5 

s^ = 0.1099, 

d.f 

• — 

2 






p 4.48 

^calc " 0.1099 

= 40.' 

773 


"'^0.05 " 

2.92 


at d.f. 

= 2 

^0.05'5> 2) = 

19.164 


1 = 

i ^0. 

,05 = 

0.3425 

^oalo^" ^0.05 

(3, 2) 



. * . Model IS 

inadequate. 



* Xq 13 a dummy variable having a constant value of +1 . 

** mol of (butenes + butadiene) formed/100 mol butane fed. 

© use average of trials 1, 9 and 10 for calculation of regression 
coefficients. 



Tables 1? and 20 shows that the response surface methodology 
has been effective in achieving improved yields. 

A linear model represented by Equation (5.3) was 
again fitted to the present set of experimental data. The 
coefficients were evaluated using Equation (5.7). The 
subregion of the response surface was therefore described by 
the equation, 


y = 8.4118 + 0.5331 - 0.3087 Xg - 0.5184 x ^ - 3.9423 x^ 

(5.1 


Comparing this equation with Equation (5.14) of the 
previous subregion, it may be seen that the sign of the 
coefficient b^ has changed from negative to positive and that 
of b^ from positive to negative. A sign reversal implied 
that the optimal values of these factors lie somewhere in the 
region of the hyper-plane of these factors covered by the two 
experimental designs. This also indicated the proximity of 
the factors to the noar-stationary region. 

2 

The experimental error variance s reported in 

</ 

Table 20 has been calculated using Equation (5.8). The 
confidence interval of the regression coefficient b was 
calculated from Equation (5.12) as, 


A»b. 


2.92 sJO. 1099/8 


0.3425 
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Table 20 indicates that the model represented by 
Equation (5.15) is grossly inadequate, implying that the 
first order response surface assumption is no longer valid. 
Hence a second order strategy had to be adopted in the further 
experimentation (see Section 5.6). 

At this [juncture, it may be recalled that the phase 
composition of the bismuth molybdate* , namely, p- and Y-phase 
was used as a qualitative variable (as mentioned in Sections 
5.4,1 and 5.5.1). During the course of investigation, the 
linear portion of the (yield) response surface of 
BigO^ . 2MoO^-AlPO^ and Bi^O^.MoO^-AlPO^ catalysts were 
investigated side-by-side in order to compare the effect of 
the qualitative factor on the yield. Hence, the experiments 
on BigO^ .MoO^-AlPO^ catalysts which were performed in parallel 
with Bi^O^ . 2MoO^-AlPO^ are discussed in the following section. 
The second order response surface of (butenes+ butadiene) 
yield was explored only after the information regarding the 
qualitative variable was available. 

5.5.3. Ontimisation of (Buten^+ Butadiene) Yield Using 
Bi^O ^.MoQ^-AlPO^ Catalyst . 

The search for the optimum was commenced by choosing 
a base level for the experimental design and following the 
procedure described in Section 5.5.1. The base levels of 
factors for the present catalyst were chosen to be the same 
as those reported in Table 16. The initial design matrix and 
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the corresponding response values ite* the yields of (butenes + 
butadiene) are shown in Table 21 < The pertinent experimental 
data are tabulated in Appendix VII for runs numbered 42 through 
49* The experiments were conducted in a randomised sequence. 

The coefficients of the linear response surface model 
were estimated using Equation (5*7) and the corresponding 
equation describing the response was, 

^ = 4.056 - 1.7716 x^ - 0.376 Xg + 0.8919 x^ - 1.1515 x^ 

(5.16) 

The experimental error variance for the present set 

of experiments was assumed to be the same as that of the 

experiments on Bi^O^. 2Mo0^-AlP0^ catalyst reported in Table 17 

under Section 5.5*1. Then by Equations (5.11) and (5.12) the 

A 

confidence interval, estimated coefficient, b^ was 

found to be + 0.6565. Since in Equation (5.16), the absolute 
value of the coefficient bg (0.376) was less than its confidence 
interval, its variation interval was purposely increased in the 
next set of experiments in order to check whether it was 
really not significant. The calculations shown in Table 21 
indicated that the model of the yield response surface represented 
by Equation (5*16) was adequate. This indicated that the 
subregion of the explored response surface was away from the 
near-stationary region, and therefore a decision was taken to 
conduct another set of experiments using the first order response 
surface strategy. 
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a?able 21. PIEST OBDER EESPOITSB SURPJiCE STRATEGY. 

PIRST MOVE ON THE RESPONSE SURE ACE. 

BXgO^.MoO^-AlPO^ CATALYST. DESIGN MATRIX, EXPERIMENTAL 
RESULTS AND TEST POR MOIEL ADEQUACY. 


t I 

1 I 

Run {Design {Rand- 
No. {Trial {omised 
{No. {Sequence 

f t 

f « 

* 

^0 

Factor Levels 

X 2 x^ 

^4 

f ¥r¥r 

{ Yield 

{ Expt. Pred. 

i y y 

42 

1 4 

+1 

-1 

-1 

-1 

-1 

6.228 

6.443 

45 

2 8 

+1 

+1 

-1 

-1 

+1 

2.678 

0.597 

44 

3 6 

+1 

-1 

+1 

-1 

+1 

5.012 

0.388 

45 

4 2 

+1 

+ 1 

+1 

-1 

-1 

0.608 

2.148 

46 

5 7 

+1 

-1 

-1 

+ 1 

+1 

4.522 

5.924 

47 

6 1 

+1 

+ 1 

-1 

+1 

-1 

3.949 

4.684 

48 

7 3 

+1 

-1 

+1 

+1 

-1 

9.965 

7.475 

49 

8 5 

+1 

+1 

+ 1 

+1 

+1 

1.325 

1.629 


N 2 

^(y - y) 

N - ( k+lT 

15.62 c 

- 3 - 

.207, 

d. 

f. = 3 


4 = 

0.513, d. 

f. = 

2 






^calc 

5.207 

" 0.513 " 

10.15 


^ 0.05 “ 

2.92 

at d . f . 

= 2 

^0.05 

(3, 2) = 19. 

164 


Ab. 

= ± 

"^0.05 

= 

0.6565 

P < 

®'o.05'3. 2 ) 



• 

* Ik 

Model 

IS adequate 



* Xq IS a dummy variable having a constant value of +1 . 

** mol - of (butenes + butadiene) formed per 100 mol of butane 
fed. 
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The linear nature of the response surface permits the 
application of steepest ascent technique to locate the next 
factor leVel combinations that piroduce a maximum increase in 
the yield on the surface. The steepest ascent direction was 
located, and the combination of factor levels lying in this 
direction were calculated according to the procedure already 
described in Section 5.5.2. These calculations are stunmarised 
in Table 22. It may be seen that the predicted values of the 
yield show an upward trend. These predicted yields were not 
experimentally verified for reasons described earlier in 
Section 5.5.2, 

The first order strategy was continued by choosing a 
new base level that gave a higher yield of (butenes + butadiene ) . 
The next set of experiments performed are described below. 

Second Mcnre Towards the Near-Stationary Begion ; In the model 
represented by Equation (5.16), the coefficient \>2 was found 
to be not significant and hence the variation interval of the 
factor ^2 was increased to 0.6 con^jared to its original value 
of 0,5 in the previous design. The variation intervals of the 
factors x^ and were also increased in order to improve the 
absolute numerical values of regression coefficients associated 
with them. 

The factor levels of the experiments in design trial 
numbered 7 of Table 21 were chosen as the new base levels, 
since the observed yield at this point was more than the 
predicted values of the yield along the response surface. 
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Table 22 i CALCULATION OP STEEPEST ASCENT. FIRST MOVE ON 
THE EESPONSB SURFACE. 

BI 2 O 5 . MoO^ -AlPO^ C AT ALYST . 




< 

I *1 

{ °G 

I 

f 

X 2 

gmin 

X 3 

mol 

*4 

mol 

! 

Calculated 

yield 



mg mol 

mol 

100 mol 

7 

1 . Basic level 


450 

1.3 

1 .0 

20 


2. Variation interval 


25 

0.3 

0.25 

5 


(unit ) 







3. Estimated slope 


-1 .772 

-0.576 

0.892 

- 1.515 


(change in yield 
per imit change in 







4 . ^ X variation interval 

0 

-44.3 

-0.128 

0.225 

-5.756 


5. Step when x. changes 
by 22 


-22 

-0.064 

0.112 

-2.9 


Roimding off 


-22 

-0.06 

0,11 

3.0 


6. The direction of 

1 

428 

1.24 

1.11 

17 

6.76 

steepest ascent 
represented by a 

2 

406 

1.17 

1.23 

14 

9.53 

series of 

3 

428* 

1.17 

1.23 

14 

7.97 

possible trials 
on it 

4 

428* 

1.05 

1.34 

11 

9.20 


* 


Held stationary to accoBsnodate variation interval for the next 
set. 
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The base levels and the variation intervals used in 
the present experimental design are shown in Table 23. The 
relevant experimental data are reported in the Appendix VTI 
denoted by runs numbered 50 through 57. The effectiveness of 
the response surface methodology was established since the 
yields reported in Table 23 were superior to the yields 
reported in Table 21 . 

A linear response was fitted to the yield data, and 
the regression coefficients were estimated using the Equation 
(5.7), The fitted equation was, 

y* = 7.94 - 2.644 x^ + 0.656 + 1.956 x^ + 1.5028 x^ 

(5.17) 

Comparing this equation with Equation (5.16), it is found that 
the coefficients and b^ have reversed their sign indicating 
that the optimum values of the factors X2 and x^ lie somewhere 
in between the corresponding factor levels used in the previous 
design (Table 16) and the present design (Table 23). 

The experimental error variance was assumed to be the 
same as that obtained in Table 20. The confidence intervals 
of the coefficients were given by Ab = + 0.3425. Thus all 

ci 

the coefficients in Equation (5.17) were found to be significant. 
The model adequacy test in Table 23 showed that the linear 
approximation of the response surface (Equation (5.17)) was 
no more valid. At this stage, it was decided to terminate 
the first order strategy. 
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Table 23. flRST OEDEE RESPONSE SURFACE STRATEGY. 

SECOND MOVE ON THE EESPONSE SURFACE. 

CATALYST: Bi^O^ .HoO^-AlPO^ BASE LEVEL AND VARIATION INTERVALS 
OF FACTORS, DESIGN MATRIX, EXPERIMENTAL RESULTS AND TEST FOR 
MODEL ADEQUACY. 



Basic level 
Variation interval 
Upper level 
Lower level 


Coded value of factors 


Run 

Random- 

Design 

No. 

ised 

Trial 


Sequence 

No. 

50 

6 

1 

51 

5 

2 

52 

7 

3 

53 

2 

4 

54 

3 

5 

55 

8 

6 

56 

4 

7 

57 

1 

8 


N 


^2 

_ tii' 

- y2) 

®ad 

N - i 

[FfTT “ 


425 

1.6 

1.25 

15 

25 

0.6 

0.4 

8 

450 

2.2 

1 .65 

23 

400 

1.0 

0.85 

7 


^2 

X3 : 

^4 


+1 

-1 

-1 

-1 

-1 

+1 

+1 

-1 

-1 

+1 

+1 

-1 

+1 

-1 

+ 1 

+1 

+1 

+1 

-1 

-1 

+1 

-1 

-1 

+1 

+1 

+1 

+1 

-1 

+1 

-1 

+1 

-1 

+1 

+1 

-1 

+1 

+1 

+1 

+1 

+ 1 

f .4137 

D .1099 


40.01 

y 


7.259 

6.469 

5.829 

4.187 

9.183 

10.787 

1.594 

2.493 

12.550 

13.387 

3.551 

5.093 

13.396 

11.693 

10.209 

9.411 


d.f. = 3 


= 0.1099, d.f- = 2 


calc 


4.4137 


0.1099 


= 40.01 


^0.05^5^ 2) = 19,164 


^ 0 . 05 ^^’ 


to 05 = 2.92 at d.f. = 2 
± -0.05 

. * ** . Model IS inadequat e . 


* X IS a dummy variable having a constant value of +1, 

** mol of (butenes + butadiene) formed per 100 mol of butane fed. 





Having established the best possible yields of 
(butenes + butadiene) in the near-stationary regions of the two 
response surfaces, namely, the response surfaces associated 
with p- and y-pbases of bismuth molybdates, the two phases were 
evaluated for their suitability as dehydrogenation catalyst. 

The second order response surface strategy was applied only to 
one of the above two catalysts found superior. The following 
section deals with the choice between the two qualitative 
variables and the second order methodology ^plied subsequently. 

5.6. Second Order Response Surface Strategy. Exploration of 
the Response at the Hear-Stationary Region . 

The first order response surface strategy applied to 
BigO^ .2MoO^-AlPO^ and BigO^.MoO^-AlPO^ led to a situation 
where the linear approximation of the response surface function 
0 was no more valid. This suggested a curvilinear nature of 
the subregion of the response surface with a possibility of max- 
imum being located close to this subregion. Therefore^ the 
second order strategy was resoirted to and a nonlinear model 
was employed to describe the response surface. 

In the first order study, the phase composition of 
the bismuth molybdate was chosen as a qualitative factor. It 
was considered appropriate to eliminate one of the qualitative 
factors at this stage of experimentation, so that the better 
amongst the two may be used for further investigation. The 
choice between the phases of bismuth molybdates is discussed below. 
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5.6;1. Screenirijs: of the Qualitative Eactors. Choice Between 
Bi^0 v.2Mo0 ^ and Bi ^O^ .MoO ^, 

Both the BigO^'SMoO^ and Bx^O^.MoO^ phases used as 
qualitative factors in the first order response surface 
methodology were found to give comparable yields of (butenes + 
butadiene) at the near-stationary regions of the two response 
surfaces (see Tables 20 and 23). The yields were slightly 
higher in the case of BigO^ .2Mo0^-AlP0^ catalyst when compared 
with the yields obtained over Bi^O^.MoO^-AlPO^ catalyst. 
Therefore, it was not possible to choose between the two 
qualitative factors on the basis of the yields obtained in the 
subregion of the two response surfaces. 

A choice between the two qualitative factors was 
made by analysing process conditions at which the near-optimum 
yields were located. A comparison of the models describing the 
two response surfaces (Equations (5.14) to (5.17)) was necessary 
for this purpose. The values of the regression coefficients 
associated with the two response surfaces are shown in Table 24. 
The influence of the process variables on the yields of 
(butenes + butadiene) were scrutinised with the help of the 
above table and subsequently relative merits and demerits of 
the two catalysts were evaluated. The influence of the process 
variables on the yields of (butenes + butadiene) are as follows; 

1 , T empe rature : (Refer Table 24). 

The temperature coefficient b^ for Bi^O^.BMoO^-AlPO^ 
catalyst changed its sign from minus to plus, when the region 
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Table 24/ BEaEBSSION COEFFICIENTS OF THE MODELS FROM THE 
FIRST ORDER RESPONSE SURFACE STRATEGY* . 


Qualitative Factor 

“T r“ 

t e 

{Average JRegression Coefficients of Factors 
{Yield { , 

1 \ ^ 0 

{ ^0 { ^2 ^3 ^4 

i « , . .... - . - . 

1. Bi20^.2Mo0^-AlP0^ 


catalyst 


First move 

3.625 - 1.7315 -0.8930 1.T055 -0.7000 

( 450 ) ( 1 . 3 ) (1.0) (20) 

Second move 

8.412 0.5331 -0.3087 -0.3814 -0.9425 

(425) ( 0 . 8 ) (1.5) ( 13 ) 


2 . BI2O3 • MoO^ -AlPO^ 


catalyst 







First move 


4.036 

-1.7716 

( 450 ) 

-0.3760 

(1.3) 

0.8919 

(1.0) 

- 1.1515 

(20) 

Second move 


7.940 

-2.6440 

(425) 

0.6560 

(1.6) 

1.9560 

(1.25) 

1 . 5028 
( 15 ) 


* Base levels are shown in parenthesis. 
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of the response surface was shifted from a base level of 450°C 
to 425°C. The sign reversal indicated that the optimum value 
of the temperature lies in between the two base levels of 
450°C and 425°C- The temperature coefficient for Bi202.MoO^- 
AlPO^ catalyst was negative in both moves, and thus the optimum 
temperature was below 425°C. 

ii. Catalyst-Peed Ratio : (Refer Table 24). 

The coefficient associated with the catalyst feed 
ratio was b2. The Bi20^.2Mo0^ catalyst required a catalyst- 
feed ratio less than 0.8 g/(mg mol/min) as indicated by the 
negative values of the coefficient. However, the optimum 
contact time for Bi20^.Mo0^ catalyst was existed somewhere in 
between a catalyst-feed ratios of 1.3 to 1.6 g/(mg mol/min). 

lii. Butane -Oxygen Ratio ; (Refer Table 24)* 

The coefficient of hydrocarbon-oxygen ratio for b^ 
indicated that the optimum butane-oxygen ratio was located 
between 1.0 and 1.5 for Bi20^ . 2Mo0^-AlP0^ catalyst. In the 
case of Bi20^ .MoO^-AlPO^ a butane-oxygen ratios in excess of 
1,5 were found beneficial. 

IV. Bismuth Molybdate -Aluminium Phosphate Ratio : (Refer Table 24) 
Lower proportions of bismuth molybdate to aliiminiiun 
phosphates are favoured in the case of Bi20^.2Mo0^. The 
regression coefficient b^ (Equation (5.15)), associated with 
Bi20^.2MoO^, had a large negative value showing that the 
optimum proportion of Bi20^.2MoO^ to AlPO^ lies far below the 





value of 13 mol/lOO mol. Howeveriin the case of Bi20^»MoO^- 
AlPO^ catalyst, it may be concluded that the optimum proportion 
lies between 15 and 20 mol/100 moli 

Prom Table 20 it may be seen that the Bi202.2Mo0^- 
AlPO^ catalyst with a BigO^ , 2Mo0^/AlP0^ ratio of 5 mol/lOO mol 
has the highest activity for the yield of (butenes + butadiene), 
whilst the optimum proportion for Bi20^.Mo0^-AlP0^ catalyst 
lies somewhere in the range of 15-20 mol/lOO mol. This 
suggested that the Bi20^.2Mo0^ was a better catalyst, since a 
catalyst with a lower proportion of bismuth molybdate to 
aluminium phosphate would result in reduced catalyst costs. 

The method of catalyst preparation also favoured the 
choice of Bi20^.2Mo0^. The catalyst preparation method in 
the case of Bi20^.2MoO^ was a sin5)le precipitation reaction 
followed by filtration of the precipitate, while the preparation 
of Bi20^.Mo0^ involved very intricate solid state reaction 
between (see Section 3.2). 

The expected value of the optimum temperature in 
the case of Bi20^.2Mo0^ (425-450°C ) further supported its 
choice over Bi20^.MoO^ since the optimum value of the temperature 
in the latter ease was below 425°C5, which could lead to a 
region where the reactants were partially converted to 
oxygenated products. 

The effects of contact time as well as butane-oxygen 
ratios favoured the choice of BigO^.aMoO^ which required a 
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lower contact time (0,8) and a butane to oaygen ratio of 1.0 
to 1#5< However/ these factors were minor compared to the 
above variables ^ The surface areas of the two catalysts 
reported in Chapter ^ deserve consideration at this stage. It 
may be seen from Table 7 (Section 4.4-) that the surface area 
of BigO^.MoO^-AlPO^ catalyst is higher than the surface area of 
Bi^O^ ,2 Mo 02-A1P04 catalyst. One would be tempted to favour 
Bi20^ .MoO^-AlPO^ on the basis of surface area data. Higher 
surface area would be a disadvantage in the present case since 
it would result in large heat generation (exothermic reaction) 
per mit volume of the catalyst bed. The control of the 
exothermic oxidative dehydrogenation reaction would be 
difficult in such a situation. 

The above arguments strongly favour the choice of 
Bi20^.2MoO^ as an active component when compared to its 
competitor, Bi20^.Mo0^, and hence it was chosen for further 
study of the yield-response surface. The second order strategy 
was applied to Bi20^. 2Mo0^-AlP0^ catalyst with a view to 
locate the factor level combinations within the experimental 
region that gaVe highest yield of (butenes + butadiene). The 
second ordeh response methodology as implemented in the 
present investigation is described below. 

5.6.2. Quadratic Model of the Response Surface . 

In the near-stationary region, the curvilinear effects 
dominate, and the functional form of the response function 0 
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(Equation (5.’1)) will have to be suitably chosen to explain 
the above nonlinearity. In the present investigation, the 
response function was approximated by a quadratic (second 
degree) polynomial suggested by Box and Wilson [3l]. Tbe 
polynomial is given byj 

y = '’o * ’^ 2^2 ’’3^3 '’11^1 '’22*2 * 

>>33^3 + + 

bgjXgXj + bg^xjx^ + bj^xjx^ (5.18) 

This model contains the linear terms in x.j , x^, x^, the 

2 2 2 2 

squared terns in x^, x^, x^, x^ and the cross product terms in 
Alternatively, the above equation can be written in 

(5.19) 

b3 b^ ] 

b-13 

^23 ^24 

^33 ^34 

^34 ^44 


the matrix notation as, 


T 


Z = bo + .felS s b, ,x 


Where 

T 


^1 


-11 


[ X., 

tb, 


'11 


'12 


'12 


'22 


'13 23 

h4 ^24 


T denotes the transpose 
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The central c 0115)031 te design developed by Box and 
Wilson [31 ] was employed to collect the experim#fital data. 
These data were used to estimate the coefficients in Equation 
(5.18). The central coD5)osite design is constructed by adding 
further factor level combinations to those obtained from the 
2^ factorial design. These extra factor level combinations 
are as follows: 


( 0 , 0 , 0 , 0 ); 


(- a, 0, 0, 0) ; 
( 0,-a, 0, 0); 
( 0, 0,-a, 0); 
( 0, 0, 0,-a); 


(+ a, 0, 0, 0); 

( 0, + a, 0, 0); 

( 0, 0,+ a, 0); 

( 0, 0, 0,+ a). 


Thus the new design contained (2^ + k + 1) factor 
level combinations with k = 4 in the present case. The value 
of a was set at 2 to make the design rotatable and the 
regression coefficients orthogonal to one another. The 
experimental design is presented in Table 25 under the heading 
"Design Matrix". The central point (or the base level) of 
the design represented by the factor level combination 
(0, 0, 0, 0) was repeated seven times in order to get an 
estimate of the experimental error variance of the optimisation 
parameter y as recommended by Cochran and Cox [38]. 

The design trial number 6 of Table 20 had the best 
combination of the factor levels with an experimental yield 
of (butenes + butadiene) equal to 54.8% » Therefore, this was 
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Table 25. SECOND OEDER RESPONSE SURFACE STRATEGY. 

CENTRAL COBIPOSITE DESIGN MATRIX AND EXPERIMENTAL RESULTS. 

. 2MoO,-AlPO. catalyst . 

2 3 3 4 



Temp 

Cat/Feed 

Butane/ 

Bi^O- . 2MoO^/AlPO. 


oc 

g/(mg mol/ 
min) 

Oxygen 

mol/mol 

^ j j 4 

mol/l00 mol 



^2 


^4 

Base level 

440 

0.7 

1.4 

5 

Variation interval 

20 

0,2 

0.2 

2 


{Design 

{Trial 

{No. 


{Random- 

{ised 

{Sequence 


Design Matrix 
Factor Levels 


Yield 

Expt. 

7 


Pred. 

7 



Continued. . . 







Table 25 (continued) 


1 ! 

2 ; 

3 

: 4 

5 

6 

7 

_ 8 _ L 

9 

10 

66 

9 

5 

+1 

-1 

-1 

-1 

+1 

9.335 

9.728 

67 

10 

4 

+1 

+ 1 

-1 

-1 

+1 

8.521 

8.486 

68 

11 

8 

+1 

-1 

4-1 

-1 

+1 

10.453 

10.073 

69 

12 

3 

+1 

+1 

+1 

-1 

+1 

9.641 

8.831 

70 

13 

6 

+1 

-1 

-1 

+1 

+1 

12.702 

12.763 

71 

14 

7 

+1 

4-1 

-1 

+1 

" t -1 

11.961 

12.323 

72 

15 

2 

+1 

-1 

+1 

+1 

+1 

12.556 

12.768 

73 

16 

1 

+1 

4-1 

+1 

+1 

4-1 

12.376 

12.328 

74 

17 

3 

+1 

-2 

0 

0 

0 

13.022 

12.444 

75 

18 

8 

+1 

+2 

0 

.0 

0 

9.854 

10.093 

76 

19 

5 

+1 

0 

-2 

0 

0 

12.528 

11.456 

77 

20 

7 

+1 

0 

+2 

0 

0 

10.079 

10.812 

78 

21 

9 

+1 

0 

0 

-2 

0 

6.415 

6.955 

79 

22 

10 

+1 

0 

0 

+2 

0 

12.789 

11.911 

80 

23 

- 

+1 

0 

0 

0 

-2 

10.995 

10.291 

81 

24 

~ 

+1 

0 

:o 

0 

+2 

11.090 

11.306 

82 

25 

-6 

+1 

0 

0 

0 

0 

10.576 

10.799 

83 

26 

12 

+1 

0 

0 

0 

0 

11.383 

10.799 

84 

27 

2 

+1 

0 

0 

0 

0 

10.892 

10.799 

85 

28 

13 

+1 

0 

0 

0 

0 

10.958 

10.799 

86 

29 

11 

+1 

0 

0 

0 

0 

10.114 

10.799 

87 

30 

1 

+1 

0 

0 

0 

0 

11.401 

10.799 

88 

31 

4 

+1 

0 

0 

0 

0 

10.116 

10.799 


* Xq IS a dummy variable having a constant value of +1 . 
** mol (butenes + butadiene) per 100 mol butane fed. 
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the best choice for the baso level of central composite design. 
Some numerical ad^justments were made in the above base levels 
of factors to accommodate their upper and Icwer levels within 
the region of experimentation. The base levels and the 
variation intervals of the factors used in the present study 
are shown in Table 25. 

The experimental data obtained were fitted to the 
quadratic polynomial given in Equation (5.18) and the 
coefficients b^, b^ , were estimated. Unnecessary 

terms in Equation (5.18) were omitted through significance 
testing. The significance of any coefficient was evaluated 
as follows [58]: 

The variance -covariance matrix of the regression coefficients, 
s (b) was evaluated first by 


s^(b) 


T 

X X 


/at t 
- b X X 


U - p 


(is)''' 


( 5 . 20 ) 


where, 

X IS the matrix of the independent variables, 
y IS the vector of experimentally observed yields. 

is the vector of estimated coefficients. 

N is total number of experiments in the design, 
p IS the number of coefficients (parameters). 

The "Student's t" value is calculated using the relation. 
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'calc 




( 5 . 21 ) 


where ,* 


^calc calculated value of "Student's t". 

IS the estimated value of the coefficient. 

J 

s(b ) IS the variance of the coefficient obtained from 

J 


Variance -covariance matrix s^(b). 


If 


'calcj ^ ^(1 - a/2, N - p)' coefficient was 

considered as not significant. Subsequently, the coefficients 
of Equation (5.18) were re-estimated after setting all the 
non-significant coefficients to zero. 

Adequacy of the model was tested as before by the 
E-test. The multiple correlation coefficient, R [60], was 
defined by the equation, 


R = 





i=1 


^ (y. 
1=1 ^ 


y)^ 


( 5 . 22 ) 


where y is given by the relation 

z: y:, 

i=1 ^ 


(5.23) 


The multiple correlation coefficient R is a measure of the 
linear association between the dependent variable y and the 
set of independent variables, 
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, ...» x^, . . . , x^Xg, . . . , X^x^ . R has a value close to 1 
for a good fxt of the data by the model. 

Experimental data reported in Table 25 were used to 
develop a model for the yield-response surface. The coefficients 
of the quadratic model represented by Equation (5.18) were 
estimated using Equation (5*6). The variance of the experimental 
error was estimated from the design trials numbered 25 through 
31. The various steps in the development of a correct model 
are as follows: 

A. All the Coefficients in Equation (5.18) were Assximed to be 
Significant . 

The coefficients were estimated using Equations (5.6) 
and (5.18). In the matrix notation these are as follows: 

= 10.777 

o 


[ -0.5876 

-0.1611 

0.1239 

0.2539] 

~ 0.11 97 

-0,0457 

0.2005 

0.1672 

-0.0457 

0.0861 

-0.0851 

0,2487 

0.2005 

-0.0851 

-0.3592 

0.3938 

0.1672 

0.2487 

0.3938 

0.0209 



Model adequacy : 

2 

Variance of experimental error s^ = 0.2866, 

d.f. = 6 
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VarxancG of adequacy = 0.5151? 

d.f. =51-15 


= 16 


calc 


¥ = = 1 -^ 974 ' 


Fo.05^''^’ 6) = 3.9254 

^calc"^ ^0 05^^^^ hence the model xs adequate. 

Multxple correlatxon coeffxcient using Equatxon (5.22) 
R = 0.9341 (x.e. the fxt is good). 


Significance of the coefficients : 

The computed "t" values of the various parameters in 
the matrix notation were as follows: 

t^ = 146.44 


t, = 

"^1 

[-23.38 

-7.50 

57.73 

11.83] 



6.65 

-1 .42 

6.23 

5.19 



-1 .42 

4.78 

-2.64 

7.75 

t, 

^11 


6.23 

-2.64 

-18.83 

12.23 



J 5.19 

7.75 

12.23 

1.16 


The tabulated value of "t" at 5% significance level 

and 16 degrees of freedom is 2.12. It was found that the 

A /\ 

calculated "t" values of the coefficients and b^^ 

( = 1.16 and 1.42) were less than the tabulated "t" 
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value of 2 i 12 , and hence the coefficients b ^2 ^44 were 

not significant 4 Their inclusion in the model may lead to 
undue regression of the experimental error, and hence it was 
considered proper to drop them from the model* 

B. Elimination of Non-Significant Coefficients . 

The remaining coefficients in Equation ( 5418 ) were 
re-estimated, and the statistical tests were applied to this 
new model. 

The re-estimated values of the coefficients were as 


follows* 






II 

10.799 




bi = 

C- 

-0.5876 

-0.1611 

1.2391 

0.2539] 



0.1175 

0 

0.2005 

0.1672 

II 


0 

0.2005 

0.0839 

-0.0851 

-0.0851 

-0.3414 

0.2487 

0.3938 



0.1672 

0.2487 

0.3938 

0 


The model adequacy was tested again. 


Model adequacy: 

Variance 

of experimental 

2 

error s = 

«/ 

d.f. 

0.2866 , 

= 6 


2 

Variance of adequacy s^^ = 0.4604, 

d.f. =31 - 13 = 18 



146 



The tabulated value of "t" at 5% significance level and 18 degrees 
of freedom is 2.10. All the coefficients in the present model 
were found to be significant. Hence the model of the response 
surface is given by, 
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y = 10.799 - 0.5876x^ - O.ISIIXg + 1.2391x5 + 0.2539x^ + 

0.1175x^ + 0.0859x| - 0.3414x5 + 0.2005x^X5 + 0.1672x^x^ 

- 0.085 1X2X5 + 0.2487x 2X^ + 0.3938x5X^ 

(5.25) 

The model fitted to the subregion of the response 
surface was examined further to locate the factor level 
combinations that result in a maxiinum yield of (butenes + 
butadiene). This is discussed in the following section. 


5.6.3. Analysis of the Quadratic Surface. Location of the 
Ontimum Factor Level Combination . 

The maximum of y was located by differentiating 
Equation (5.25) with respect to each x in turn and setting 
them equal to zero flO], The resulting set of simultaneous 
algebraic equations were solved to obtain the optimum value of 
X that maximised y. 

Upon difforontiation with respect to x^, Equation 
(5.18) yields. 


^ = 



= 0; i=1, 2,3,4. 


( 5 . 26 ) 


Or, in the matrix notation, 



2b,1 


X = 


0 


(5.27) 
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The solution of the above equation gave the vector x , that 

” —max 

maximised the yield y: 


-max 


1 ^-I^T 
“ 2 - 11-1 


^e 


(5.28) 


Substitution of Xq m Equation (5.20) would give the optimum 


value of yield as, 


y. 


max 


y^ = + b. X, + b.. x„ 

•' e o — 1 — e — e — 1 1 — e 


(5.29) 


where , T represents the transpose of the vector. Thus, the 
optimum (coded) values of factor levels calculated using 
Equation (5. 28), and tho corresponding natural values are shown 
in Table 26. The maximum value of yield of (butenes + butadiene) 
was calculated using Equation (5.29). 

Presence of any saddle points or ridges in the fitted 
response surface (Equation (5.25)) would result in the location 
of a false optimum. It was therefore considered worthwhile to 
analyse the fitted response for the presence of any saddle 
points or ridges within the experimental region. This analysis 
is as follows I 


5.6.4. Simplification of the Quadratic Response Surface — 
Canonical Analysis t29. 601. 

By suitable transformation of the coordinates, 
Equation (5.19) may be reduced to the form, 


(5.30) 
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Table 26. ‘‘OPTITOI' VALUES OP FACTOR LEVELS AND THE 
YIELD OP (BUTENES + BUTADIENE). 

CATALYST; Bi20^ . 2Mo02-AlP0^. 


Factor 

t 

! Code 

J Name 
< 

-? 

j Coded 
! Value 

f 

f 



J Natural Value 

f 

! 

f 

-...f - - . -- 

Temperature 

^1e 

-3.306 

374 “o'" 

0 at aly st/feed ratio 

^2e 

-1.595 

0.38 mol/ (mg mol/min) 

Butane/oxygen ratio 

^3e 

2.089 

1.82 mol/mol 

Bi 20 ^. 2 MoO^/AlPO^ ratio 

^4e 

1.576 

8,15 mol/ 100 mol 

'Optimum' yield = 13. 

39 mol 

100 

(butenes + butadiene)/ 

mol butane fed. 


* Lies outside the expenmental region. 
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called xts canonical fona. This was achieved by translation 
and rotation of the coordinate axes as follows! 

The first order teims, i.e, b-jX, in Equation (5.19) 
or ( 5 . 26 ) were eliminated by shifting the origin to the value 
x„ predicted by Equation (5-29). Thus, the new coordinate 
system (x^ , x^, x^ , x|) was defined by, 


xj 


or 


X - X 

1 le 


X - ^e 


(5.31) 


Introducing this in Equation (5.19), 

y = -1 -e^ -e^^ + x^) 

Subtracting Equation (5.29) from Equation (5.32) one gets, 

y -Je = in 2 e +4 -^ii 2') +2''' 2' 


(5.32) 


(5.33) 


The expression in the brackets of Equation (5.33) is identical 

T 

to zero because p remult iplic at ion of Equation (5.28) by ,x' b^^ 

results in the following equation: 


X' ■b-,1 

— —11 — e 


(x’^ ) 


while the equality, 

T ^ ,T ,T ^ 1 / 

— e — 1 1 — “ — — 1 1 — e 2 — —1 
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IS obtained by the multiplication of the elements of the 

T T 

respective matrices. Also observing b^x' =2' 
simplify Equation (5.33) as, 




(5.34) 


The translation of the coordinates of the present 
response surface represented by Equation (5.25) to the was 
carried out to eliminate the linear terms namely, x^, x^, 
and x^ in it. Hence Equation (5.25) was converted into the 
matrix notation using Equations (5.19) and (5.24), from which, 
the translation of the coordinates to x was affected. Hence, 
Equation (5.34) for the present case was 


0.1175 

0 

0.2005 

0.1672 

0 

0.0839 

-0.0851 

0.2487 

0.2005 

-0.0851 

-0.5414 

0.3938 

jo . 1672 

0.2487 

0.3938 

0 


(5.35) 


where , 


X' 


x^ + 3.31 

Xg + 1.60 

x^ - 2.09 

x^ - 1.58 


(5.56) 


The expansion of Equations (5.54 ) , or (5.35) proved 
that the cross product terms like x^x^, xJjx^ etc. were still 
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present in the equation* These were eliminated from the 
equation by rotating the coordinate axes defined by (xJJ, 
x^, x|) to align it with the principal axes of the quadratic. 

This method, known as "reduction by orthogonal transformation" 

[29, 31, 40, 60], was adopted to affect such a change. In the 
orthogonal transformation, the new coordinate axes (X^ , X^, X^, 
X^) are related to the old coordinates (x^ , x^, x^, x^) through 
the transformation, 

X = U"'* X’ = (x - Xg) (5.37) 

-1 T 

where U is called the unitary matrix. The property U = U 
of the unitary matrix was used in Equation (5.37). The 
component vectors of this transformation matrix were the same 
as the eigenvectors of the matrix b^ ^ . The equation of the 
response surface will then be reduced to the form represented 
by Equation (5.30), where >» 2 » A 3 and are the 

eigenvalues of b^ ^ . 

Thus, the eigenvalues and the eigenvectors of the 
matrix b^^ given in Equation (5.24) were evaluated. The 
determinant of this matrix was evaluated first and was found 
to be equal to +0.006458. It may be observed that the matrix 
b^-j is real and symmetric, with positive value of the 
determinant. The method of Jacobi was therefore used to 
calculate the eigenvalues and the eigenvectors of b^ ^ matrix 
[ 33 ]. A FORTRAN programme of this algorithm is given in 
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reference [33] as SUBEOUTEHE JCOBI. The eigenvalues and the 
eigenvectors of ^ in Equation (5.25) computed by this 

subroutine were 4 

= 0.49066 

Ag = 0.15536 

(5.38) 

X 5 = -0,66528 
= -0.12068 

and 


0.5222 

0.5228 

-0.0852 

-0.6684 

0.3254 

-O. 8 O 7 I 

0.2709 

- 0.4115 

0.4109 

0.2271 

0.7880 

0.3982 

0.6727 

-0.1541 

-0.5463 

0.4747 


respectively. Thus, the canonical form of Equation (5.25) is 
given by, 

y - 13.39 = 0.49066 + 0.15536 - 0.66528 X^ - 

0.12068 X^ (5.40) 

where the coordinate system (X^ , Xg, X^) is related to 
the original coordinates (x^ , x^r x^) by the following 
relations 
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**“ — 

1 ) 


■ 

0.5222 0.3254 0.4109 0.6727^ 


X^ + 3.31 ; 

X 2 


0.5228 -0.8071 0.2271 -0.1541 


Xg + 1.60 

1 h 


1-0.0852 0.2709 0.7880 -0.5463 

i 

i ^ 

i 

Xj - 2 . 09 , 

\ 


-0.6684 -0.4115 0.3982 0.4747 

1 1 

^4 - ''•58| 


(5.41) 


It IS seen that the eigenvalues and /\ 2 were both positive 

while X 4 were both negative. It is evident that the 

yield of (butenes + butadiene) will show an increasing trend if 
one changes the factor levels along the positive or negative 
directions of and X 2 axes while any movement along X^ and X^ 
directions would decrease the yield. This established that the 
optimum point x (Table 26) computed using Equation (5.29) was 
not a true maximum, but a saddle point. The absolute values 
of the eigenvalues were comparable in magnitude, Hence , the 
response surface was an elliptical hyper-saddle surface. Factor 
level combinations, which predict yields higher than that 
predicted at the saddle point x^were calculated as follows: 

Since a movement along X^ and X^ directions resulted in 
a decreased yield, it was decided to constrain such a movement 
by setting the values of X^ and X^ in Equation (5.40) to zero. 
Equation (5.40) now becomes, 

y - y^ = 0.49066 X^ + 0.15556 x| 


(5.42) 
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The above equation represents a plane parallel to the X^X2 
coordinate plane passing through the saddle point Xg* 
constraints X3 = 0 and = 0 result in relations. 

= -0.0852(x^ - x^g) + 0.2709 (x2 - + 

0,7880(x^ - - 0 . 5463 (x 4 - x^^) 

= 0 

X^ = -0.6684(x^ - x^g) - 0.4115 (x 2 - Xg^) + 

0.3982 (x^ “ ^ 0.4747(x^ - x^^) 

= 0 


(5.43) 


(5.44) 


In Equation (5.42), X^ and X2 are given by the following 
relations : 


X^ = 0 . 5222 (x^ - x^^) + 0.3254 (x 2 - Xg^) + 

^4 ■ -*^4e' 


0.4109(x 3 - x^g) + 0 . 6727 (x^ - x^g) 


(5.45) 


Xg = 0.5228(x^ - x^g) - 0 . 807 l(x 2 - X2g) + 

0.2271(^:3 - ^3e) - 0 - 1541 (x 4 - ^^43) 

Equations (5.43) to (5.46) were obtained frcm Equation (5.41). 
Solving Equations (5.43) and (5.44) for x, and x^, one gets, 

(X, - X, ) = 0.6856(x^ - x,^) . 0.1627 (x 2 - XjJ (5.47) 

j j 6 

(x^ - x^g) = 0.8329(x., - x.,g) + 0.7305 (x2 - X 2g) (5.48) 



Substituting Equations (5.45) to (5.48) in Equation (5.42) and 
simplifying, the following relationship for the yield of 
(butenes + butadiene) was obtained, 

y' = 0.9602 (x^ - x^^)^ + 0.5042(x2 - Xg^)^ 

+ 1.0321 (x^ - x^g)(x2 - X2g) (5.49) 

The response values calculated using the above 
relationship within the experimental region are presented in 
Table 27. The factor levels were decoded using the data given 
in Table 25. It can be seen from the table that the hyper-plane 
described by Equation (5.49) lies outside the experimental- 
range of butane to oxygen mtio. It was, therefore, decided to 
map (numerically) the entire response surface for various 
factor level combinations so that the best feasible combination 
of factor levels may be chosen. 

The yields were calculated at different factor level 
combinations using Equation ( 5 * 25 )* The results of such 
calculations are presented in Appendix V. 

It is apparent from these tables (in Appendix V) that 
an increase in the Bi 20 ^. 2 Mo 0 ^ content of the catalyst results 
in an increase in the total yield of (butenes + butadiene). 

For example, the yield of (butenes + butadiene) predicted at a 
temperature of 460 °C, W/F = T.2 g catalyst/mg mol/min of feed, 
butane to oxygen ratio = 1 .75 and bismuth molybdate to 
aiiiTmTinTm phosphate ratio of 0.1 is 16.06 percent. This value 
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Table 27. CANONICAL ANALYSIS OP THE EESPONSE SURE ACE - LOCATION 
OP THE OPTIM PACTOR LEVEL COMBINATION FOR 
BigO^ . 2MoO^ -AlPO^ C AT ALYST . 


f — 

t 

Temperature } 

°C jg/( 

1 

W/P 

mg mol/min) 

1 7^ — 

P 4 H 10 /O 2 jmol 

I mol/mol {per 

1 » 

r 

BigO^ , 2 MoO^ } 
100 mol AlPO^ j 

Yield* 

400 

0.3 

1.98 

9.7 

14.57 

400 

0.4 

2.00 

10.5 

15. 16 

400 

0.5 

2.02 

11.2 

16.01 

400 

0.6 

2.03 

11.9 

17.11 

400 

0.7 

2.05 

12.7 

18.46 

400 

0,8 

2.06 

13.4 

20.07 

400 

0.9 

2.08 

14.1 

21.92 

400 

1.0 

2.10 

14.8 

24.03 

410 

0.3 

2.05 

10.6 

15.85 

410 

0.4 

2.07 

11.3 

16.71 

410 

0.5 

2.08 

12.0 

17.81 

410 

0.6 

2.10 

12.8 

19.17 

410 

0.7 

2.12 

13.5 f 

20.78 

410 

0.8 

2.13 

14.2 

22.64 

410 

0.9 

2.15 

15.0 

24.76 

410 

1 .0 

2.17 

15.7 

27.12 

420 

0.3 

2.12 

11.4 

17.62 

420 

0.4 

2.14 

12.1 

18.73 

420 

0.5 

2.15 

12.9 

20.09 

420 

0.6 

2.17 

15.6 

21 .71 

420 

0.7 

2.19 

14.3 

23.58 

420 

0.8 

2.20 

15.1 

25.70 

420 

0.9 

2.22 

15.8 

28.07 

420 

1.0 

2.23 

16.5 

30.69 


* Yield 


mol (butenes + butadiene) formed per 100 mol of 
butane fed. 
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needs to be expemnentally verified so as to check its validity 
since these factor level combinations lie on the boundaries of 
the range of experimental variables covered in the second 
order design (see Table 25). 

It can be seen from the design trial 8 of both Tables 
17 and 19 that similar combinations of factor levels resulted 
in poor experimental yields of (butenes + butadiene) of 4.17 
and 4.13 percent respectively. Therefore, 16.06 percent yield 
predicted by Equation (5.25) at the above-mentioned factor 
level combinations is erroneous and hence it was not considered 
to be the best value. 

Next the predicted yield at catalyst compositions 
of less than 7.5 mol Bi^O^ . 2MoO^/l 00 mol AlPO^ were examined. 

At this stage, it may be pointed out that a very high butane 
to oxygen ratio (1.5 or larger) is not desirable since the long 
tenn activity of the bismuth molybdate will be reduced 
drastically. The catalyst has the best long term activity when 
it IS operated at a butane to oxygen ratio of 1.0 to 1.5. 

Hence, a catalyst composition which has reasonably uniform 
yield over a wide range of butane to oxygen ratio would be the 
best from an industrial point of view. Erom a comparison of 
the two catalysts at composition levels of 5.0 and 7.5 mol 
Bi20^.2Mo0^/l00 mol AlPO^, it can be seen that at a temperature 
of 420°C and a butane to oxygen ratio of 1.0 or less the 
yield for the former is better. Prom a comparison of the 
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catalysts at composition levels of 2.5 and 5.0 mol Bi 20 ^. 2 Mo 0 ^/ 
100 mol AlPO^, it IS seen that the latter catalyst gives hotter 
yield of (butenes + butadiene) at butane to oxygen ratios of 
1.25 and higher. Brom these arguments, the best catalyst 
composition was assumed to be 5.0 mol Bi20^ . 2MoO^/l 00 mol AlPO^. 

Effect of temperature on the yield of (butenes + 
butadiene) on the 5.0 mol Bi 20 ^. 2 MoO^/lOO mol AlPO^ catalyst 
was examined next. It can be seen that for the butane to 
oxygen ratios between 0.75 to 1-75, the yield of (butenes + 
butadiene) decreases with an increase in the temperature. A 
temperature of 420°C was chosen to be the best value in 
preference to 400°C in order to avoid the formation of 
oxygenated compounds during the oxidative dehydrogenation 
reaction. 

Having chosen a temperature of 420°C, and a catalyst 
composition of 5.0 mol Bi20^.2Mo0^/l00 mol AlPO^ experiments 
were conducted on this catalyst at various catalyst to feed 
ratios and butane to oxygon ratios in order to verify experi- 
mentally the predicted values of (butenes + butadiene) yields. 
These experiments (discussed in the following chapter) are 
reported in Appendix VII as runs numbered 89 through 106. A 
comparison of these experimental yields with the yields 
predicted by the second order model (Appendix VI) supports 
the conclusions drawn in the above analysis. 
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5 . 7 . Conclusions . 

The response surface methodology has been successful 
in establishing the best conditions for the oxidative dehydro- 
genation of butane to butenes and butadiene over bismuth 
molybdate — aluminium phosphate catalysts. The first order 
approximation of the response surface, and the subsequent 
application of the steepest ascent technique to the fitted 
response surface gradually increased tho yield of (butenes + 
butadiene). The first order strategy also led the process 
conditions like temperature, catalyst-feed ratio, butane- 
oxygen ratio, catalyst composition to their near-optimum 
conditions. 

The lack of fit of the linear model to the response 
surface indicated that the curvature effects could no more be 
neglected. This indicated the nearness of the stationary 
region of the response surface. 

At this stage, a decision was taken regarding the 
qualitative variable, i.c. the type of bismuth molybdate phase 
used in the preparation of tho catalysts. The two phases, namely, 
Bi20^.2MoO^ and Bi^O^.MoO^ were found to be equally active in 
the conversion of butane to butenes and butadiene. However, 
Bi20^.2MoO^ was preferred to Bi20^.Mo0^ since the former was 
stable at higher temperatures, and also its proportion in the 
catalyst (compared to AIBO^) was less. Therefore, the near- 
stationary region of only Bi 20 ^. 2 Mo 0 ^-AlB 0 ^ catalyst was 
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studied by making use of central composite design and 
subsequently, a quadratic model was fitted to the experimental 
data. A canonical analysis of this quadratic polynomial 
indicated that the response surface was not a multidimensional 
hill, instead, it was an elliptical hyper-saddle surface. 

Exploration for the optimum factor level combination 
along the rising section of the hyper-saddle proved futile, 
since better yields were predicted only at butane to oxygen 
ratios above 2.0. Such an optimum was not feasible, since, the 
bismuth molybdate catalyst would deactivate rapidly at high 
butane to oxygen ratios. 

As an alternate step, the yield of (butenes + butadiene) 
was evaluated at all possible combinations of factor levels in 
the experimental range. Prom an observation of these values, 
a temperature of 420°C and a catalyst composition of 5.0 mol 
BigO^ .2MoO^/lOO mol AlPO^wero chosen to be the best values 
since at these factor level combinations, the yields of 
(butenes + butadiene) were uniformly better at various butane 
to oxygen ratios. 
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CHAPTER 6 

MODEILING OP OXIDATIVE DEHYDROGEN ATI OF REACTION 

6.1. Introduction . 

In the yield studies using response surface methodology, 
the oxidative dehydrogenation of butane was investigated without 
any prior knowledge of the mathematical model describing the 
dehydrogenation process. A black-box approach was used and the 
localised regions of the response surface were approximated using 
simple algebraic expressions which do not give any insight into 
the physical nature of the dehydrogenation reaction. In the 
present chapter, an attempt has been made to arrive at a most 
plausible mechanism based on the physical nature of the system 
and the experimental observations. An understanding, in the 
microscopic level, of the possible steps involved in the 
oxidative dehydrogenation reaction aids in the developme>"t of a 
mathematical model describing the oxidative dehydrogenation 
process. 

In heterogeneous reactions, one or more of the reactants 
are chemisorbed on the solid surface of the catalyst prior to 
their conversion to the final products. Some possible modes of 
butane adsorption on the solid catalysts are discussed b low. 
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6.2* Mechanism of Butane Adsorption . 

Studies of exchange reactions of alkanes with deuterium 
(Pafshall et.al* [92]) indicated that the adsorption of alkanes 
on the metal oxides is dissociative, i.e., 

ECH^ H 

1 ; 

BCH^ + 2S — ^ S + S (6.1) 

where S is an adsorption site on the catalyst and R is an alkyl 
group. An alternate formulation suggested by the same authors 
involved the adsorption of alkane on a site, given by, 

R H 

\ / 

RH + S S (6.2) 

followed by migration of hydrogen atom. 

Mobile lattice oxygen present in bismuth molybdates 
suggest a redox mechanism (proposed by Mars and Van Krevelen 
[85]) for oxidative dehydrogenation of butane. In this scheme, 
the interaction of the butane molecule with lattice oxygen in 
the catalyst leads to the formation of butene and water which 
on further reaction with lattice oxygen forms butadiene. 

On the basis of reaction mechanisms discussed above, 
various plausible models have been developed in the following 


section. 
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6.3. Reaction Rate Models . 

In the derivations given below^ 
r IS the rate of reaction, 
p is the partial pressure of reactant, and 
K is the adsorption equilibrium constant, 

with subscripts A, B, C, 0 and W representing butane, butene, 
butadiene, oxygen and water respectively. 

6.3.1. Redox Model . 

The oxidative dehydrogenation of butane to butadiene is 
represented by the following reaction scheme; 


A + Gat (ox) 

■I 

B + Gat (red) 

(6.3) 


B + C at ( ox ) 


C + Cat (red) 

(6.4) 


Cat (red) + 


C at ( ox ) 

(6.5) 


where Cat (ox) represents the oxidised state of the catalyst and 
Cat (red) denotes the reduced state of the catalyst. 

let, 9 be the fraction of oxidised sites in the 
catalyst. Then, by using the above reaction scheme, one can 
write the reaction rates as, 
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rg = 

"•ox = 

By a balance for oxygen, i.e., = ^Gat(red)’ 

the following equation is obtained, 





Hence, the reaction rates are now given by, 

” “"iPa *'2Pb ’^3Po 

^^iPa “ ^2Pb^^3Po 
' ICiPa + Vs + ^3Po 

k,p| + kgP^ + k^p^ 

The rate of consumption of oxygen and the rate of generation of 
water can be calculated from the equations, 

( 6 . 9 ) 

^0 “ 2 
% = 


( 6 . 6 ) 


( 6 . 7 ) 


(6 . 8 ) 




( 6 . 10 ) 
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6.3.2. Langmuir-H ins he Iw o od Model . 

Assuming irreversible surface reactions, viz., 

k. 


* * 
A +0 


* * 
B + 0 


‘1 


* * 
B + H 2 O 


* * 
C + HgO 


( 6 . 11 ) 

(6.12) 


(* represents adsorbed species), 


the following expressions may be derived when surface reactions 
are rate controlling. 




(1 + + KjgPg + K^Pq +\/KqPq + K^Py/) 


(6.13) 


(kiKAPA - kgKgPg) i/KqPq 


'B 


^A^A ^B^B ^wPw^ 


(6.14) 


^2 ^bPb^^oPo 


^aPa ^'^bPb ^cPc ■*'V^oPo ^Pw^ 


(6.15) 


The oxygen consumption and water generation rates can 
be calculated using Equations (6.9) and (6. 10). 


6.3,3. Rideal-Elev Mechanism . 

It IS possible that only one of the two reactants 
(either butane or oxygen) is adsorbed, and the other reacts from 
the gas phase. In a situation, where butane is adsorbed, and 
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oxygen is In the gas phase, the surface reactions are given by 


* 1^ ^1 * 

A + ^2 ® + ^ 2 ^ 


(6.16) 

* 1 * 

B + ^02 — ^ C + H 2 O 

(6.17) 


Assuming irreversible surface reactions to be rate controlling, 
the following rate expressions can be obtained^ 




(6.18) 


(k,K^p^ - »^2 K^Pb)Po 

^B 1 + K^p^ + KgPg + KqP(. 


(6.19) 




^ 1 + K^P^ + KgPg + KqPq 


( 6 . 20 ) 


Rates of oxygen consumption and water generation are again 
given by Equations (6.9) and (6.10). 

Similar expressions can be derived for the case where 
the reaction takes place between adsorbed oxygen species and 
butane in gas phase. Such a mechanism involving butane in gas 
phase is very unlikely. 

Experimental data obtained in the laboratory, were 
used to test if one or more of the above models are plausible. 
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6(4. Experimental Observations . 

The analysis of the response surface of (butenes + 
butadiene) yield reported in Section 5.8 indicated a catalyst 
composition of 5 mol Bi20^ . 2MoO^/ 1 00 mol AlPO^ to be the most 
suitable one. Experiments v/ere conducted on a randomised order 
on this catalyst at 420°C and at various levels of catalyst to 
feed ratio and butane to oxygen ratio. These data are reported 
in Appendix VII as runs numbered 89 through 106. 

The partial pressures of butene and butadiene 
measured at different catalyst -feed ratios and butane to oxygen 
ratios are shown in Pigs. 22 and 23. It can be visualised 
from Pig. 22 that the partial pressure of butene has a possible 
maximum in the range of catalyst to feed ratios employed. This 
is because butene is an intermediate product of two consecutive 
reactions. The partial pressure of butadiene (Pig. 23) reached 
a plateau at high values of catalyst to feed ratios, which is 
possibly due to cracking. 

Point selectivity, defined as the ratio of the rates 
of formation of any two species [99] , was used to identify the 
plausible models describing the oxidative dehydrogenation 
process. Results and discussion of these aspects are given in 
the following sections. 
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Fig. 22 : Partial pressure of butenes vs. W/F ratio tor 81203. 2M0O3- 

AIPO4 catalyst. 


Temperature = 420®C • Butane/ oxygen =0-75 
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Fig. 23 : Partial pressure of butadiene vs. W/F ratio for 61203. 2Mo03“ 

AIPO4 catalyst . 
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6.5. Identification of Plausible Models — Method of Point 
Selectivity. 


Point selectivity was defined as the ratio of the rate 
of formation of butene to that of butadiene. These ratios for 
the various models proposed in Section 6.3 are as follows: 


Redox Model 


^3 1 

kg ^A^B " ^ 

langmuir-Hinshelwood 

"“b 

^I^A 

Pa 

Model 


■ 

Pb 

Eideal-Eley 

^B 


U _ 

Mechanism 


^2^3 ' 

Pb 

These equations can be 

reduced 

to the 

form 



a -b 

PaPb 


1 


( 6 . 21 ) 


( 6 . 22 ) 


(6.23) 


(6.24) 


where k' is a constant and a = b = 1 for both Langmuir- 
Hinshelwood and Eideal-Eley models. Thus all the models 
describe the response (point selectivity) for the situation 
when a = b = 1 . However, discrimination is possible when the 
exponents are different from unity. 

The rates of formation of butene and butadiene were 
obtained by graphical differentiation of partial pressure Vs. 
catalyst to feed ratio curves which are shown in Pigs. 22 and 
23. The parameter k' in Equation (6.24) was evaluated at 
different values of a (= 0.5 to 2) and b (= 0.5 to 1) by linear 
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regression and the confidence limits [60] of the same are 
presented in Table 28. It can be inferred from the table that, 
since the parameter k* is positive for all cases it is not 
possible to reject any models. The combination (a = 0.5, 
b = 0.5) has the least estimate of error variance. Taking this 
as the best model (in the least squares sense), the plausibility 
of most of the combinations of a and b cannot be ruled out. 

Hence, in order to achieve better understanding about the 
models, a multiresponse technique has been proposed as a scope 
for future study. 

6.6, Pa-pameter Estimation — Multiresponse Situation . 

Models constructed in Section 6.3 on the basis of 
redox, Langmuir-Hmshelwood and Rideal-Eley mechanisms, in each 
of the oases were a set of nonlinear ordinary differential 
equations. The two step nature of the oxidative dehydrogenation 
reaction required that all the reliable responses be considered 
in evaluating the parameters [62], In the present investigation, 
butene and butadiene measurements were found to be most reliable 
and hence they were used in the parameter estimation. 

Bard [15] analysed the problem of parameter estimation 
in systems of nonlinear ordinary differential equations and 
concluded that the use of surface minimisation methods (such as 
Davidon-Pletcher-Powell method given in subroutine DEMP and 
Pletcher-Eeeves conjugate gradient method given as subroutine 
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Table 28. POINT SILEOTIYITY IVIBTHOD OP MODEL DISCEIMINATION. 
PAEMBTER ESTIPIATBS MD THEIR CONFIDENCE INTERVAL. 


S.No. 

Exponents 

a b 

Parameter k' 

+ Confidence 
Interval 

Estimate 

of 

Variance 

Total 

Residi 

1 

0.5 

0.5 

0.195 + 0.059 

0.559 

4.84 

2 

0.5 

1 .0 

0.194 + 0.102 

1.155 

10.43 

5 

1 .0 

0.5 

0.016 + 0.006 

0.827 

7.32 

4 

1 .0 

1.0 

0.014 + 0.009 

1.365 

12.37 

5 

1.5 

0.5 

0.0012 + 0.0007 

1.069 

9.57 

6 

1.5 

1.0 

0.0011 + 0.0008 

1.534 

15.93 

7 

2.0 

0.5 

0.00009 + 0.00005 

4 .466 

11.49 

8 

2.0 

1.0 

0.00008 + 0.00006 

1 .669 

15.17 
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EMCG in IM Scientific Subroutine Package [117] ) would lead to 
poor convergence of the parameters. This was confirmed in the 
present investigation when using the determinant minimisation 
criterion suggested by Box and Draper [30, 62] « Hence a simpler 
criterion suggested by lapidus and Peterson [76, 95, 107] given 
as, 

_ H ^ /s, 2 H 2 

2—. (Pb " Pb ^ ^Pc “ Pc ^ ~ minimum (6.21) 

i=1 1 1 1=1 i i 


was used instead. 

The algorithm used in the parameter estimation is 
presented in the form of a flow diagram in Pig. 24. Marquardt's 
nonlinear least squares algorithm [84], as implemented by Ball 
and Groenweghe [l3, 14], was used for the minimisation of sum 
of squares. A listing of this subroutine is found in Kuester 
and Mize [72] as SUBROUTINE BSOLVE. The system of ordinary 
differential equations describing the model was integrated 
using the method suggested by Gil [67b]. A computer programme 
for the same has been developed by Krishnamurthy and Sen [67a]. 

The derivative matrix, /\ required in the Marquardt's 
algorithm may be computed as the solution of the sensitivity 
equations [101, 107], with initial conditions as given below: 


>13 = 


n 


Is 


X 




"dk 


( 6 . 22 ) 


STAKE 



Assume a set of starting 
values for parameters* 


Reset the 
values of the 
parameters to 
the improved 
values 


Integrate the set of 
differential equations 
(model) using Gil’s method 
xo obtain the predicted 
values of the responses. 


Calculate the corrections 
using Marquardt's least 
squares minimisation 
algorithm and evaluate the 
error sum of squares, 0 , 


{ 0-0 

'^min 


STOP 


Pig. 24. Plow Diagram of the Algorithm Used in 
Multiresponse Parameter Estimation. 
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with initial conditions, 


^,,( 0 ) = 0 


1 0 » t f 

J = u 2i ..., P 


(6.23) 


where , 


^ f derivative with respect to time, t, 

IS the derivative of the response with respect to the 

* d P* 

various parameters given as 

0^ is the function defining the i-th response (model) 
p^ IS response, p^, p^, etc. 
k- IS the parameters to be estimated 

J 

n IS the number of responses in the model (Five in the 
present study, viz., the partial pressure of butane, 
butene, butadiene, oxygen and water) 

P is the total number of parameters to be estimated. 

In the matrix form Eq.uation (6.22) may be represented as. 


A = aa +1 


where , 


A = -Tr y 

r-'’A 

t— t— 

II II 

H -73 

B - 

1 Jk* 

1 = 1 

1 


(6.23) 


2, ..., n 

(6.24) 

2, n 


2 , . . . , n 

(6i25) 

2 , . • . , P 
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The details i^egarding the construction of these sensitivity- 
equations hawe been described in Appendix VI. Por n responses 
and P parameters to be estimated, n(P+1 ) sensitivity equations 
(a set of SDiaultaneous differential equations) have to be 
solved in order to evaluate the derivative matrix analytically. 
The effort required in the numerical evaluation of the 
derivative matrix is even more. 

The algorithm described above and shovm by a flow 
diagram in Pig. 24 can be used for obtaining the parameters of 
the multiresponse model. An attempt has been made in this 
direction. However duo to the limitations of the computer 
speed (IBM 7044) and time, it could not be pursued at the 
present time . 

6.7. Conclusions . 

Using point selectivity ratio, defined as the ratio of 
the rate of formation of butene to that of butadiene, an attempt 
has been made to screen the various reaction models represented 
by sets of ordinary differential equations. In view of the 
acceptable nature of the parameter k’ (that it should be 
positive in the confidence region) and a small ratio of the 
error variances between the best (in the least squares sense) 
and the poor model, barring the model represented by S.Uo. 7 of 
Table 28 for which the error Variance was usually large, it was 
deemed not possible to discriminate amongst the various models. 



A suggestion for future work has been included in Section 6.6 
for estimating the parameters for multiresponse situation which 
hopefully gives better discrimination amongst the various models. 
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CHAPTER 7 

SUI, 7 .IARY AND CONCLUSIONS 

7 . 1 . Summary and Conclusions . 

The present v/oik was started with a feasibility 
investigation of the oxidative dehydrogenation of n-butane on 
a molybdic oxide catalyst. Though the yield of (butenes + 
butadiene) on this catalyst was low, the feasibility of one 
step dehydrogenation of butane to butadiene was established. 

In order to improve the yield of (butenes + butadiene), 
bismuth, which has a promoter effect in the oxidative dehydro- 
genation of butenes, was incorporated in the catalyst in the 
form of bismuth molybdate. Three phases of bismuth molybdates, 
namely, a, p and y with chemical compositions Bi20^.3MoO^, 
Bi20^.2Mo0^ and Bi20^.Mo0^ were used in the present investigation. 
The dehydrogenation catalysts were prepared by mixing them 
physically with ortho-alxominium phosphate in varying proportions. 

The physical and chemical characteristics of the 
prepared catalysts v/ere investigate^ first in order to understand 
their structure and thermal s-jjabilities . The X-ray and 
infrared ap-alyses of the bismuth molybdates established that 
the a -phase was likely to be inactive since it had a monoclinic 
structure with edge shared octahedra. The presence of corner- 
shared octahedral molybdenum in p- and Y“l 3 ismuth molybdates 
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(as inferred from infrared and X-ray analyses) indicated that 
they were promising active materials for the dehydrogenation 
catalyst* The calcination of these p- and y- forms with 
o-AlPO^ did not indicate any structural Changes in them. In 
the case of p-bismuth molybdate, calcination with o-AlPO^ 
inhibited, partially, the decomposition of the former, while 
the calcination of y-bismuth molybdate with o-AlPO^ indicated a 
tendency towards preferred orientations of certain crystal 
planes , 

An investigation of the thermal stabilities of p- 
and Y“^is®u-th molybdates indicated that the p -bismuth molybdate 
underwent a reversible decomposition at temperatures close to 
its melting point (about 700°C), while the y-form underwent 
an irreversible transition y — y’ (at temperatures above 
600°G ), y^ being a poor catalyst active material. Both these 
transformations occurred beyond the range of temperatures 
400-500°C used in the oxidative dehydrogenation reaction. 

Hence the p- and y -bismuth molybdates catalysts on o-AlPO^, 
prepared by the above procedure, were used. 

Response surface methodology was applied to optimise 
the yield of (butenes + butadiene). fractional 

factorial experiments were used to explore the linear region 
of the yield -response surface using temperature, catalyst to 
feed ratio, butane to oxygen ratio and proportion of the 
bismuth molybdate in aluminnun phosphate support, as variables. 
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In addition, the phase of bismuth molybdate was used as a 
qualitative variable. Simultaneous exploration of the response 
surface with two qualitaxive variables indicated that they had 
similar activity for oxidative dehydrogenation of butane, the 
maximum experimental yields being 12.6^ and 14.8% for y- and 
p-bismuth molybdate aluminium phosphate catalysts respectively. 

The near-stationary region of the yield response 
surface was explored for p-bismuth molybdate -aluminium 
phosphate catalyst only. This catalyst was preferred to y- 
bismuth molybdate -aluminium phosphate, since all the thermal 
(phase) transitions of p -bismuth molybdate were reversible 
while those of Y-bismuth molybdate were irreversible. 

Simplicity of preparation of p -bismuth molybdate also favoured 
its choice over Y -bismth molybdate . 

The near -stationary region of yield-response surface 
using p-bismuth molybdate -alviminium phosphate catalyst was 
explored by employing a central composite design. The experi- 
mental data, fitted to a quadratic (second degree) polynomial, 
was found to represent the data adequately. A canonical 
analysis of the fitted response surface indicated that the 
surface was an elliptical hyper-saddle with the rising section 
of saddle passing through the principal planes (of the 
quadratic) lying outside the region of experimentation. Hence, 
the response surface was mapped numerically and a region 
where the predicted yields of (butenes + butadiene) v/ero 
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uniformly good (about 10-125^) in the commonly employed range 
of butane to oxygen ratios (0.5-1 -5) was chosen for modelling 
the oxidative dehydrogenation reaction. 

Models based on point selectivity (defined as the 
ratio of rate of formation of butenes to that of butadiene) 
indicated that the redox, Langmuir-Hinshelwood and Rideal-Eley 
models were all equally plausible. 

7.2. Suggestions for Future Work . 

On the basis of the present investigation, the 
following aspects are recommended for future work: 

1. Catalysts prepared by coprecipitation of BigO^, MoO^ and 
AlPO^ could bo studied. 

2. In the preseit study it was observed that yl^^smuth 
molybdate exhibits preferred orientation in crystal planes 
and the decomposition cf p -bismuth molybdate is partially 
inhibited, in the presence of aluminium phosphate. Hence, 
the thermal stability of the bismuth molybdates in the 
presence of aluminium phosphate provides an interesting 
area for research. 

3. Discrimination amongst various models based on redox 
langmuir-Hinsholwood and Eideal-Eley mechanisms was not 
achieved in the pi’esent study using point selectivity 
ratios. Estimation of parameters in these models using 
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raultiresponso techniques and discrimination of these models 
based on the estimates of parameters are suggested for 
future study. 

4. The canonical analysis of the response surface of 
(butenes + butadiene) yield indicated that better yields 
are obtained at higher butane to oxygen ratios (>2). The 
stability of the catalyst under such high butane -oxygen 
ratios needs to be investigated. 

5. The bismuth molybdates undergo structural transformations 
(participation of lattice oxygen) during oxidative 
dehydrogenation, which may lead to a slow change in the 
structure of a freshly prepared catalyst. Hence, a study 
of the aging of the bismuth molybdate -aluminium phosphates 
will be important from an industrial viewpoint. 
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iiPPENDIX I 

CHEMICIE MAI.YSIS OP CATAI.YSTS 


A. Analysis of Bismuth Nitrate Bi(NO, ), .5H„0 and Bismuthyl 
Nitrate BiONO .:H^O [1261" ^3^5-^2 

EPTA solution ; 0.05 M. Dissolve 18.613 g of disodium-dihydrogen- 
ethylenediaminetetra-acetate dihydrate (BDH Analar) in 1 liter 
of water* 

Bismuth ion solution : 0.05 M. About 6.06 g of bismuth nitrat 
was dissolved in about 25 ml concentrated nitric acid and made 
to 250 ml with water. 

Xylenol orange indicator ; 0.05 g in 10 ml of water. 

Procedure ; 25 ml of bismuth ion solution were diluted to about 
150 ml with Water. pH was adjusted to 2 to 3 by cautious 
addition of dilute ajmnonia. 3-4 drops of xylenol orange 
indicator solution were added, and the solution was titrated 
against standard EDTA solution. Titrant was introduced slowly 
when the r^d colour of the original solution started fading. A: 
the end point the colour changed sharply, from red to orange 
yellow. 

1 ml EDIA s 10.4500 mg Bi 

Purity of bismuth compounds estimated as per above procedure 
are given below; 
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Bismuth nitrate (BDH Analar) - 98*97% 

Bismuthyl nitrate 81*09% Bi20j 

B. Analysis of Molybdic Acid [68, 126] « 

Oxine solution (acid) : 4 g oxine (8, hydroxyquinoline , Ward 
Blenkinsop Sc Co., London, A.R. Grade) were dissolved in 8.5 ml 
of warm glacial acotic acid. It was poured into 80 ml of wate 
and diluted to 100 ml. 

Oxine solution (ale oholic ) ; 2 g of A.R. grade oxine were 

dissolved in 100 ml of ethyl alcohol. 

Procedure : 0.1 g of molybdic acid were dissolved in dilute 
ammonia. A few drops of methyl red indicator were added to get 
yellow coluration. The solution was acidified with 2 U sulphuric 
acid to methyl red colour (neutral). A few drops of 2 N 
sulphuric acid weire then added in excess. Further, 5 ml of 
ammoniima acetate solution were added and diluted with water to 
50-100 ml. Solution was heated to gentle boiling and acetic 
acid solution of oxine was added until supematfjit liquid becomes 
perceptibly yellow. The suspension was simmered for 5 minutes 
while stirring. The precipitate was filtered through a G3 
sintered glass crucible and was washed with hot water to free t 
from the oxine reagent. It was dried in an oven at 130-140°C lO 
constant weight, and weighed as MoOgCOglig 011 ) 2 . 


Molybdic oxide content of 
molybdic acid 
(Reachim U.S.S.R.) 


86.934 wt % 

(88.874% as per chemical 
foimila). 
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C . MoO^ Content of P- and Y-Bismuth Molybdates [43 , 97 3 . 


The bismuth molybdate to be analysed was dissolved in 
nitric acid, and a few drops of methyl red indicator were added. Mix 
about 15 ml of 2 wt % BDTA solution, and 5 ml of buffer solution 
(120 ml of 50 wt % ammoniTom acetate solution with 80 ml 

acetic acid and 80 ml water), then dilute to 100 ml. The 
solution was brought to boiling, oxine solution (in acetic acid as 
prepaired in Section B) was added in slight excess, and the 
suspension was simmered for 2-5 minutes. The suspension was 
filtered through a G3 sintered glass crucible ; precipitate was 
washed with hot water until the washings were colourless. The 
precipitate was dried at 130-140°C to constant weight, and 
weighed as Mo 02 (CgHg 0 N) 2 * 


wt % M 0 O 3 found in the 
P -bismuth molybdate 
prepared (Bi 20 ^. 2 MoO^) 


(38.19 as per the chemical 
formula) . 


wt % M 0 O 3 found in the 
Y -bismuth molybdate 
prepared (BigO^ .MoO^ ) 


23.89^ 

( 23 . 6 % as per the chemical 
formula) . 
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APPENEEX II 

CAIIBEIIION OP THEHfiOG OUPLES 

Chromel-alumel thermocouples were fabricated from 20 
and 24 BS gage theirmocouple wires. Thermocouple beads were 
made by twisting the wires and fusing them over a oxyacetylene 
flame with borax as the flux. The length of the thermocouples 
was about 50 cm. Thermocouples made from 24 gage wir^were 
placed in stainless steel sheath of about 4 mm I.D. and were 
aged by heating them in a tubular furnace at 550-600°C for a 
period of two days. 

Thermocouples were calibrated by measuring the e.m.f . 
generated across the thermocouple at the boiling points of water 
naphthalene, mercury and sulphur. Boiling tubes for determining 
boiling points of naphthalene, mercury and sulphur were 
constructed according toN.B.S. specifications [102], The 
results of e.m.f. measur(anent at different boiling points are 
given below; 

Naphthalene : Pure, white, sublimed grade from Riedel, Germany 

Mercury : Tripple distilled quality from City Chemical 

Corporation, New York. 


Sulphur 


Resublimed grade 
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Table 29 * CALIBRATION OP THERMOCOUPLES. 


....II 

C ompound 

1 

oC 

— 

}C hr ome l-Alum 1 1 
{Standard e.m.f. 
[128] 

i mY 

t 

1 

E.m.f, Generated by 
Thermocouples - mV 

2 3 4 5 

6 

Water 

100 

4.10 

4.05 

3.85 

4.05 

4.00 

3.95 

4.05 

Naphth- 

alene 

217.96 

8.86 

8.90 

8.50 

8.55 

8.75 

8.75 

8.80 

Mercury 

356.8 

14.59 

14.50 

14.35 

14.55 

14.40 

14.35 

14.35 

Sulphur 

444.6 

18.28 

18.20 

18.18 

18.25 

18.19 

18.10 

18.10 


The e.m.f.s were measured using a Honeywell Model 2720 Portable 
Test Set potentiometer. 

Thermocouples 1, 3 and 6 were used for the measurement 
of temperatures in the reactor. The e.m.f. generated by these 
thermocouples were close to the standard values and hence the 
standard tables [128] were used for e.m.f .-temperature conversion. 

Thermocouple (l) was used for contr*olling the 
temperature of reactor furnace and thermocouples (3) and (6) 

Were used to measure the catalyst bed temperatures 
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APPENDIX III 
ORSAI analysis 

Oreat analysis was employed to analyse carbon dioxide, 
oxygen and carbon monoxide in the gas mixture. Reagents 
employed in different cases are given be low i 

Carbon dioxide * 

A 25-30 wt % aqueous solution of Analar potassium 
hydroxide was used [126] . 

Osz^en. 

Sodium dithionite solution was used. This reagent 
was preferred to alkaline pyragallol reagent for the following 
reasons [126]; The absorption of oxygen though slow, takes 
place at lower temperatures also. No carbon monoxide is evolved, 
and the reagent is unaffected by the unsaturated hydrocarbons. 

The absorption reaction is, 

^^2®2^4 ^2 ^2^ NaHSO^ + NaHSO^ 

The reagent is prepared by dissolving 80 g L.R. grade sodium 

5 

dithionite in 400 cm of distilled water and 48 g Analar 

5 

potassium hydroxide in 64 cm distilled water. The two 
solutions are mixed and stored in a tightly stoppered bottle. 
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0 &rb 9 n monoxide < 

Acid cuprous chloride solution was used to avoid the 
absorption of ethylene [88, 126 ]. The solution was prepared by 
dissolving 100 g of L.R. grade cuprous chloride in 500 cm^ of 
concentrated grade hydrochloric acid, and diluting it to 

1 liter. The solution was stored in a stoppered bottle 
containing copper tubes* extending through the entire depth of 
the solution. The reaction is 

OuCl + 2 H 2 O + CO — GuCl. 2 H 2 O.CO 

The absorption is never entirely complete. For this reason, 
the gas was always contacted with a fresh solution to remove 
traces of unabsorbed carbon monoxide. The gas mixture after 
the absorption of GO was passed over KOH solution to remove 
HCl vapours. 

The gas absorption pipette in which the solution was 
used, also contained copper strips to avoid the formation of 
cupric ions. 

Confining liquid . 

Saturated aqueous sodium sulphate solution containing 

a little of sulphuric acid was used. The solution was prepared 

3 

by dissolving 200 g anhydrous sodiim sulphate in 800 cm*^ water 
and 40 cm^ of concentrated sulphuric acid [88]. This reduces 
the solubility of the gases in the confining liquid. 
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APPENDIX IV 

PROCESS GALCUIATIONS 


Carbon, hydrogen, oxygen and nitrogen atomic balances 
were used to calculate the yield frcm the experimental obser- 
vations recorded by the gas chromatograms and Orsat analysis. 
The calculation procedure is as follows; 

Basis of calculation : A gas mixture containing hydrocarbons 
in the same proportions as the proportions obtained in the gas 
chromatogram i.e. 

moles of ethane, 
moles of ethylene, 
moles of propane, 
moles of propylene, 
moles of n-Butane, 

Cg moles of 1 -butene, 

Orj moles of trsns -2-butene, 

Cg moles of cis-2 -butene , and 
Cg moles of 1, 3-butadiene. 

C^ : C^ : ... ; Cg are chosen in the same proportion as in the 
chromatogram. 
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The peak areas calculated from the chromatograms were 
multiplied by their respective response factors (Chapter 4, 

Table 5 ) , and the resulting value or a suitable multiple thereof 
were chosen as the values of C^, C^, Gg. 

Calculation procedure . 

Let, 

Cq moles of methane, 

C^Q moles of carbon dioxide, 

C^^ moles of carbon monoxide, 

0^2 iQoles of oxygen, 

N moles of nitrogen, and 
W moles of water vapour 

be present in the total moles of geis chosen in the basis of 
calculation. 

A carbon, hydrogen and oxygen atomic balances lead 
to the following relations: 

CABB = 2(G^ + Cg) + 3(C^ + C^) + 4(Cg + Cg + + Gg + Gg) 

+ C^Q + C^^ + Gq (A.IV-1) 

HTL = 6(C^ + C^ + Cg) + 4 C 2 + 8(C^ + Cg + C^ + Cg) 

+ lOCg + 2W + 4 Cq (A.IV-2) 

= 2C^q + C^^ + 20^2 + ^ 


OXY 


(A.IV-3 ) 
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The carbon to hydrogen ratio is fixed by the type of hydrocarbon 
used/ Therefore I for butane gas, 

HYD = 10 (CABB) (A.IV-4) 

The butane to oxygen ratio (HYOX) is known, since its value is 
fixed during the experiments. Hence, for butane gas used, 

4 (CAEB) = 2 (OXY) . HYOX (A.IV-5) 

The nitrogen to oxygen ratio (RNITROX) is also fixed dunng the 
exp eriments . The re f ore , 

I = (OXY) . (MITOX) (A. IV -6) 

Thus we have six equations and nine unknowns (GARB, HYDG, 

OXY, Cq, ^12’ ^ three more relations 

required, can be derived from the results of the Orsat analysis. 

let 0RS1 , 0RS2 and 0RS3 be the volumes of carbon 
dioxide, oxygen and carbon monoxide gases measured per 100 
volumes of gas drawn into the Orsat apparatus. With the water 
vapour pressure (at temperature of analysis) represented as VP, 
the following relationships may be obtained. 


mol COg 

mol unabsorbable gas 


0ES1 

100 - (0RS1 + 0RS2 + 0RS3 + ^ . VP) 

(A.IV-7) 


mol 0^ 


0BS2 


mol unabsorbable gas _ (qrsi + 0RS2 + 0RS3 + - '^R) 


(A.IV-8) 
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mol CO 

mol unabeorbable gas 


0RS3 

100 ~ (0RS1 + 0RS2 + ORS5 + ^ . VP) 

(A.IV-9) 


In the basis of calculation chosen, the total moles of gases 
that are not absorbed by the Orsat reagents are given by. 


moles of 

unabsorbable = 0^ + Cg + ^4 ^6 +0^ + Cg + Cg 

+ Cq + N 

= TOTMOL + Cq + N (A.IV-IO) 

Substituting the above equation in Equations (A.IV-7) to 
(A.IV-9)» one gets# 


(ORS1) . (TOTMOL + Gq + N) 

Cio = (ORSKYC ) 

(ORS2) . (TOTLIOL + Cq + N) 

C^2 = ( ORSHYC ) 

(0RS3) . (TOTIIOL + Cq + N) 

C^^ = (ORSHYC) 

where , 


(A.IV-11) 


(A.IV-12) 


(A.IV-13) 


(TOTMOL) = C^ + Cg + + C^ + C^ + Cg + C^ + Cg + Cg 

(A.IV-14) 

and 

(ORSHYC) = 100 - (ORS1 + ORS2 + 0RS3 + ^ . VP) (A.IV-15) 
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Equations (A«IV-1 ) through (A*IV-6) and (A.IV-11) through 
(A. IV -15) now consist of eleven equations v/ith GARB, HYD, OXY, 
Cq? ^10’ ^11' ^12* * TOTMOL^ ORSHYG as the unknowns, and 


hence they may bo solved. 


Algebraic manipulation of the above equation give 


the follcwing equations: 

1 4 . 0RS1 + 0RS2 (p , (0RS1 + 0RS3) . (RIJITOX) 

_ ORShYC J ^0 4 ( ORSHYG) i (HYOX) 


1 . (c;j?B) 


= - (GARBOH) - j * 


(A.IV-16) 


4.Gq - 10 (GARB) + 2.W 

2(0RS1 -f ORS2) + 0RS5 
(ORSHYG) 


- (HYDGEN) 


0 


0 


(k.IV-17) 


2(ORS1 + 0RS2) + (ORS3).(RHITOX) 1 

4 (HYOX). (ORSHYG ) " 2 (HYOX ) 


(GARB) + W 


= - • (2(0RS1 + 0RS2) + 0ES5) (A.IV-18) 

where , 

GARBON = 2(0^ + G^) + 3(G^ + C^) + 4(G^ + Cg + G^ + Gg + 0^) 

(A.IV-19) 

and 

HYDGEN = 6(G^ + + 0^) + 4 G 2 + 8(0^ + Gg + G^ + Gg ) + lOC^ 


(A.IV-20) 
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The computational algorithm is as follows: 

1. Obtain the values 0^ through 0^ from the chromatogram. 

2. Obtain the values ORS1, 0RS2 and ORS3 from Orsat analysis. 
Find the vapour pressure of water, VP, at the temperature 
of Orsat analysis. 

3. Calculate (TOTMOl) from Equation (A.IV-U). 

4. Calculate (ORSHYC) from Equation (A. IV -15). 

5. Calculate (CARBON) from Equation (A.IV-19). 

7. Calculate (HYRGEN) from Equation (A.IV-20). 

8. RNITOX and HYOX are the nitrogen to oxygen and butane to 
oxygen ratios respectively, used in the experiments. 

9. Solve Equations (A.IV-16) through (A.IV-18) for (GARB), 

W and Cq. 

10. Compute (OXY) from Equation (A.IV-5). 

11. Compute N from Equation (A.IV-6). 

12. Compute C^q, and C^^ from Equations (A.IV-ll) through 
(A.IV-13). 

13. Compute (HYD) from Equation (A.IV-4). 

The yields of the various products were calculated 

from the following definition. 

mol products formed 
% yield = -jQQ ^^tane fed 

The above algorithm was implemented on a IBM 7044 


computer. 
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APPENDIX V 

MAP OP RESPONSE SURFACE 

Table 30. LOCATION OP OPTIMTO FACTOR LEVEL COMBINATION. 
SECOND ORDER PBSPONSE SURFACE. 

Bi20^ . 2MoO^ -AlPO^ CATALYST . 
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APPENDIX VI 

SENSITIVITY EQUATIONS Of ORDINARY DlfPEREI'TIAL EQUATIONS 



The differential 

equations representing the 

reaction 

scheme can be written as, 



dPA 

at " 

£(t), k) 

c+ 

II 

o 

II 

> 

o 

(A.VI-1 ) 

dpj 

dt " 

D(t), k) 

= 0) = 0 

(A. VI -2) 

dPc 

dt ~ 

k) 

p^Ct =0) =0 

(A. VI -3) 

dpo 
dt ~ 

k) 

Po(t = 0) = p^ 

(A.VI-4) 

dF" ^ 

)> k) 

1 1 = 0 ) = C 

(A. VI-5) 


The number of responses in the present case are, n = 5. These 
are the partial pressures of butane, butene, butadiene, oxygen 
and water. 

Let the number of parameters to be estimated be P and 
let the data be obtained at r discrete values of time. The 
sensitivity equations for the above set of diffearential 
equations are given by the relations : 
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A.jCo) = 0 

1 = 1, 2, . . . , n 

D = 1, 2, •••; P 

where, 

a Pi 

/^ij “ dk 

J 

The above set of equations can be represented in the 
as, 

^ ( 0 ) = 0 

i = 1,2,»««,n, ” 1,2,#*«,n 

1 = 1 ,2, . . .,n, a = 1 ,2, . . . ,P 

Hence the sensitivity equations in the present case 
as follows : 

d (£Za, . ^ 

dt 'Okj' - aPA ■ ‘>>^3 

^ ,*P0 3 

3Po ■ ^ 

3=1, 2, 




+ 


dpg 

3 0« 


^Pw 


Opo ■ 

a%. 


-"dk. 


c>Pc, 

^Dk ' 



• n ^ 0 


C)f!, 


3-3 ^a k. 


(A.VI-6 ) 


(A. VI -7) 

matrix form 


(A.VI-S; 


(A. VI -9) 

(A.VI-10) 

can be written 


(A.Vl-11 ) 
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a_ f-Il) 

dt 


“ *'A) + 
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d 0r 
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00^ dv-n 00' 
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D0 
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DPv/ * ^k^ Dk 


= 1» 2, » • * t B 


(A. VI-12) 


(^uT’) = 


^ Pa 


ppa (?£S) . ^ . {^) 

(.■Tir"/ + -v_ * vw ■' .1 -n- '■?ik 


■ 3Pb * ‘■'kj i)PC ' ^^ 3 ' 


'SPq 


,7)Po, /^iPw^ >""0 
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W 1 OPb 


'■ak, aPo 
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3Po 
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dPo 
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The partial derivative in the above equations can be evaluated 
from the Equations (AiVI-l) to (A.VI-5 ) defining the model. 

The. Equations (AiVI-ll) to (A. VI-16) v/ill have to be solved 
together with the Equations (A. VI-1 ) to (A, VI -5) in order to 
obtain the elements in the derivative matrix^-^ . Hence, in the 
present case, which has five responses, a total of 5(1+1*) 
number of equations v;ill have to be solvedto obtain the elements 
in the matrices A and B given in Equation (A.VI-8). 



APPENDIX ¥II 


sxperi:3:i:pai data 

In the following tables, experimental data have bee 
reported on dry basis. In these: 

Conversion is defined as moles of butane reacted per 100 mo3 
of butane fed, 

Yield is the moles of product formed per 100 moles of butane 
reacted (bute = butenes, bdn = butadiene). 
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